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I. 


SUMMARY 


The broad objective of this program was to investigate injection, mixing 
and combustion processes using gaseous propellants over a range of operating 
conditions originally specified for the Space Shuttle Auxiliary Propulsion 
System. The end objective was to relate injector and chamber design parameters 
to combustion performance, chamber heat flux and combustion stability in the 
form of a step-by-step design handbook applicable to any selected operating 
condition or gaseous propellant combination. 

The scope of the program was to include analysis, design, fabrication, 
and experimental tasks to develop the injector/ chamber design criteria. The 
principal efforts in this program were to be devoted to evaluating various 
injector element configurations on the basis of single element cold flow and 
hot fire testing. Based on this evaluation and multiple element cold flow 
tests, full scale injectors were to be designed, fabricated and tested to 
verify the design criteria for high performance, low heat flux to the thrust 
chamber, and stability. The program was divided into five tasks: Task I - 
Injector Design Analysis, Task II - Cold Flow and Hot Flow Modeling, Task III - 
Single Element Hot Firing Evaluation, Task IV - Full Scale Injector Evaluation 
and Task V - Data Evaluation. 

Task I was a concept screening and element evaluation effort which 
culminated in the selection of four element types for single element cold flow 
evaluation. The element concepts selected were: (1) the shear coaxial ele- 

ment, (2) the premix element, (3) the external impingement element and (4) 
elements for micro-orifice injectors. Each element concept included design 
geometry variations so that a total of 74 unique element designs were selected 
for detail design during this task. In addition, an analytical study was con- 
ducted to evaluate the effect on chamber heat flux and performance of a low 
mixture ratio barrier at the outer edge of the injector. The results of this 
analysis revealed that chamber heat flux could be lowered considerably by a 
low mixture ratio barrier, but that the performance penalty would be higher 
than if fuel film cooling were used to achieve the same chamber heat flux 
and wall temperature. 

Task II was concerned with the fabrication, cold flow testing and 
analysis of both the single elements and multiple elements designed during 
Task I. The testing consisted of sampling the flow field in the chamber with 
a multi -element probe which was sequenced to measure both local total pressure 
(mass flux) and composition. From these measurements a mixing efficiency could 
be determined at any axial position in the chamber, and radial and circumfer- 
ential gradients were evaluated to obtain chamber compatibility data. An 
evaluation of this data led to the following conclusions: 

1. The shear coaxial element is a relatively low mixing rate concept 
but results in a fuel-rich composition at the outer wall. 
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2. Premix and external impingement elements are high mixing rate con- 
cepts. Nominal mixture ratio composition at the outer wall is characteristic 
of the premix concept, and the external impingement element when the optimum 
injection momentum ratio is selected. 

3. Multi -element cold flow mixing efficiency was not noticeably dif- 
ferent from comparable single element mixing efficiency. 

In addition to the cold flow evaluations, limited combustion testing was con- 
ducted using the same measurement techniques. These experiments were con- 
ducted with the swirl coaxial element and conclusively demonstrated that com- 
bustion retards the mixing rate. 

Task III involved design, fabrication and testing of single element 
thrusters to obtain combustion performance, chamber wall heat flux and sta- 
bility data to compare with the cold flow data obtained during Task II. A 
total of 76 tests were conducted with 7 injector element designs which repre- 
sented variations of the following element concepts: (1) premix, (2) triplet, 

(3) coaxial and (4) swirl coaxial. All of these injectors were evaluated with 
3 different chamber geometries to permit evaluation of chamber length and con- 
traction ratio changes. The results of this single element evaluation corre- 
sponded qualitatively with the Task II single element cold flow testing in 
both combustion performance (mixing efficiency) and chamber heat flux (compo- 
sition at the wall). In addition low frequency instability was encountered 
under certain conditions on all element concepts except the premix. 

Based on the results of Tasks I through III, 2 full scale injectors 
(premix and triplet) were designed during Task IV for maximum cornbustion per- 
formance and a third (triplet) injector was designed with a low mixture ratio 
(barrier cooled) outer row of elements to investigate chamber heat flux effects. 
The experimental data obtained during this task substantiated the design cri- 
teria developed during the earlier tasks: the injectors designed for high com- 

bustion performance achieved energy release efficiencies from 98.4 to 100% with 
a 7.6 (cm) 3.0 in. long chamber. The barrier cooled tests also substantiated 
the results of the Task I analysis, i.e., the performance penalty is excessive 
for the corresponding reduction in chamber wall heat flux. In addition, com- 
bustion was marginally stable during the premix tests and unstable at high 
frequencies during about 25% of the triplet tests. 

The output of Task V, Data Evaluation, was two gaseous injector com- 
bustion models. The first model used the test data and correlated it directly 
with injector/chamber design parameters which are recognized from both theo- 
retical and empirical standpoints as the controlling variables. This empiri- 
cal model has the advantages of (1) inherently being the most accurate pro- 
cedure for gaseous injectors which are to be designed within the operating 
envelopes and using the same propellants as this program, and (2) simplicity 
in the calculation procedure itself. However, it lacks generality since it 
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stressed utilization of the test data and did not concentrate on quantifying 
the mechanistic causal relationships of the mixing/combustion process itself. 
The second modeling approach had the objective of understanding the mixing/ 
cotrbustion process to the maximum extent possible, using both available theo- 
retical knowledge and new techniques suggested and developed from close obser- 
vation of the test data. It is somewhat more complex than the empirical model, 
but it has quantitatively characterized the mixing/combustion process for 
gaseous propellants so that it is general in nature and can handle all gaseous 
propellants and operating conditions. Both of these models have been sum- 
marized into step-by-step design procedures for gaseous injectors with the 
required information displayed in charts, graphs, and tables for clarity of 
presentation. 
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II. 


INTRODUCTION 


In recent years, considerable effort has been devoted to increasing the 
level of knowledge concerning injection and combustion processes in rocket 
engines. The goal of these efforts has been to provide sufficient information 
for the reliable design of injectors having high performance, as well as being 
compatible with the chamber wall. Most of the work in this area has involved 
either liquid-liquid or gas-liquid propellant phases, with many different types 
of injector elements being characterized (e.g., like and unlike impinging 
doublet, triplet, pentad, and coaxial types). Characterization of these ele- 
ments and complete injectors was accomplished by the use of cold flow tests 
with propellant simulant fluids to determine mass and mixture ratio distribu- 
tion profiles. Determination of droplet size distribution was obtained by 
various experimental techniques. Using the experimental cold flow results, 
and correlating such results with hot firing tests, a comprehensive theoreti- 
cal base has been obtained which permits the design of injectors having high 
performance as well as minimum interactions with the chamber wall. 

The work described above has principally been associated with liquid- 
liquid propellant combinations. More recently, some attention has been 
directed to gas liquid combinations. However, to date no comprehensive effort 
has been directed toward developing similar technology which would permit the 
reliable design of gaseous propellant injectors. This program addresses that 
technology void. 

The scope of this program encompassed the analytical screening and 
rating of a multitude of injector elements with the end objective being a com- 
prehensive set of design criteria. These criteria should allow the future 
design of injectors for high performance, stable operation, and maximum com- 
patibility with the thrust chamber, for any selected operating condition or 
gaseous propellant combination. The above was accomplished by single and 
multiple element cold flow tests to investigate the many variables which 
affect the mixing process. As part of this laboratory effort, sampling of 
a combustion gas flow field was accomplished to permit the effect of combus- 
tion on the mixing process to be explored. This was followed by single ele- 
ment hot fire testing to enable a correlation to be made with the previously 
obtained cold flow testing. The multiple element mixing process was then 
investigated on a hot fire basis. The last task was to correlate all the 
data generated on this and other programs using gaseous propellants and 
formulate a design approach for injectors and combustion chambers. 

The significance of the work accomplished on this program can be divided 
into three primary areas. First is the comprehensive set of cold flow data 
which characterized the mixing process for 74 unique element designs over a 
range of operating conditions. Secondly, the influence of combustion on the 
mixing process has been determined for four different element types, using 
gaseous hydrogen, gaseous oxygen propellants. Thirdly, the first two items 
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have provided the experimental bases for the development of gaseous injector 
modeling procedures which can be used to determine optimum design parameters 
for most selected operating conditions and propellant combinations. 

Documentation of the significant work is contained in the following 
sections of this report. In addition a supplemental Design Handbook is 
attached to the end of this report to specify the step-by-step calculations for t 
injector modeling procedures. This Design Handbook contains only the procedural 
steps with most of the required information displayed in charts, graphs, and 
tables for clarity of presentation. In addition, an example problem is included 
so that the reader can be assured that he is using the models correctly. 



III. INJECTOR DESIGN AND ANALYSIS 

The objective of this task was to evaluate and screen potential injec- 
tion elements leading to the selection of four basic element concepts which 
have the potential for being high performing and, at the same time, meet other 
required injector design criteria for the Space Shuttle Attitude Control Propul- 
sion System. The following sections discuss the results of this evaluation 
and selection process, as well as the cold flow hardware designed to incorpo- 
rate the selected elements. Also as part of this task film vs barrier cooling 
studies were conducted and combustion stability characteristics were considered. 

A. CONCEPT SCREENING 

Potential injector design concepts for gaseous propellants evalu- 
ated in the screening process are shown in Table I. The left-hand column, 
labeled "Elements", catalogs basic element configurations. Design deriva- 
tives considered for each element are tabulated in the center column while the 
right hand column lists element geometry variations. The first two columns 
are starred {*) for those design concepts having test information available 
in the literature prior to the starting date of this contract. The applicable 
program has been keyed in the right-hand column. Examination of Table I 
reveals that there is at least some industry experience for seven of the 
twelve elements identified (column 1). The most in-depth experience appears 
to be with the shear elements. 

The screening results are shown in Table II, where estimated ability 
to meet each of the defined evaluation criteria is designated by a number between 
1 and 5 (1 is minimum ability and 5 is maximum ability). The elements are cata- 
loged according to (1) external impingement, (2) premix, (3) shear, (4) micro- 
Qf'-jfTCG snd (5/ rGVGrsG flow concopts. In 8 qtoss sghsGj this rsting includGS 
the orifice geometric derivatives listed in Table I for each design, since the 
particular orifice geometry which was estimated to be best for that element 
was used in determining the rating. For example, the F-O-F triplet has a 
non-circular jet (tri-slot) which was estimated to result in better compati- 
bility, higher potential for decoupling the mixing/reaction processes, and a 
lower design/fabricability rating than a comparable circular-orifice triplet. 

In addition, the rating was based on equal thrust-per-element with the excep- 
tion of the micro-orifice concept. 

Using an equal weighting factor for each criteria, the evaluation 
results are summed in the right-hand column. Four element groupings resulted 
in at least one design derivative which was assessed a rating of 4. Based on 
this rating chart, the four element types and derivatives selected for further 
study were; 
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TABLE I 

GAS/GAS INJECTOR ELEMENT TYPES AND DERIVATIVES 


Elements 

Basic Design 
Derivatives 

Orifice Geometry 
Derivatives 

EXTERNAL IMPINGING 

* Doublets 

* Like Doublets 


Circular 


* Unlike Doublets 


Nonci rcul ar 

Continuous Sheets (1) ( 

* Triplets 

* F-O-F 


Circular (1) 
Noncircular (6) 


O-F-0 


Continuous Sheets 

Quadlets 

O-F-O-F (Trans) 


Circular 


0-0-F-F (CIS) 


Non circular 

Pentads 

F-F-O-F-F 


Circular 


O-O-F-0-0 

• 

Nonci rcul ar 

PREMIX 

* Triplets 

* F-O-F 


Circular 

O-F-0 


Noncircular (7) 

Pentads 

F-F-O-F-F 


Circular 


O-O-F-0-0 


Nonci rcular 

* Swirl 

* 0-Axial , F-Tang 


Circular (8) 


F-Axial, 0-Tang 


Noncircular 

SHEAR 

* Coaxial 

* Fuel Annulus 


Circular (1 ),(6),(7),(1 
Nonci rcular 


Oxidizer Annulus 


Swirler (6) , (7) 
Convolutions (9) 

Showerhead 

Both Propellants 


Circular ( ) 
Nonci rcular 


One Propellant Plus 
Porous Face Plate 

• 

Continuous Sheets ( ) 

MICRO-ORIFICE 

* Showerhead 

* Both Propellants 


Noncircular (3), (4), (5) 

Triplets 

* Like Doublets 
Unlike Doublets 

F-O-F 

O-F-O 


Noncircular (lO) 


KEY: 

* Elements and operation derivatives that have been tested with GH 2 /GO 2 
( ) Program definition for particular orifice geometry derivative 

1. Contract NAS3-14347, "H-0 Catalytic Ign. and Thruster Program" (TRW) 

2. Contract NAS 3-14353, "High Pressure Reverse Flow APS Engine" (Bell) 

3. Contract NAS8-21052, "Advanced Injector Concepts Investigation" (Aerojet) 

4. Contract NAS 8-20672, "Stability Characterization of Advanced Inj." (Aerojet) 

5. Aerojet IR&D Program 

6. Contract NAS 3-14352, "Space Shuttle Auxiliary Propulsion System" (Rocketdyne) 

7. Contract NAS 3-14354, "H-0 APS Engines" (Aerojet) 

8. Aerojet IR&D Program 

9. Aerojet IR&D Program 

10. Contract NAS 8-26188, "Space Shuttle Main Eng. Definition Study, Phase B (Aeroje 

11. Contract NAS9-8285, "Apollo Optimized SPS Injector (OMS)" (Aerojet) 
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1 . Shear Elements 

a. Coaxial 

b. Increased shear area coaxial 

c. Swirl coaxial 

2. Premix Elements 

a. Triplet 

b. Pentad 

3. External Impingement Elements 

a. Triplet (impingement angle tt/4 rad (45°), 
tt/3 rad (60°), and tt/2 rad (90°)) 

4. Micro -Orifice Elements 

a. Parallel sheet 

b. Like Doublet 

B. ELEMENT EVALUATION 

1 . Shear Mixing Elements 

Shear mixing elements mix by turbulent diffusion through the 
shear mixing layer. In order to characterize this element analytically, the 
Aerojet Turbulent Mixing and Chemical Reactions Computer Program was used to 
parametrically investigate various coaxial element characteristics. This pro- 
gram models the interacting effects of mixing and reaction by solving the 
combined energy, streamwise momentum, and species equations using the Prandtl 
boundary layer approximation. A nominal (66.6 N) 15 Ibf thrust coaxial element 
was studied and the fuel annular area was varied from half that of the oxidizer 
area to twice the oxidizer area while mixture ratio was varied from one to 
sixteen. 


The results of this study are shown in Figure la*, where the 
mixing efficiency is plotted as a function of relative injection area and 

mixture ratio for a chamber length to element diameter ratio of 10. This 
analysis draws the conclusions that decreasing fuel annulus is in the direction 
of higher En,. Experimental data pertinent to the coaxial element is shown in 
Figure lb and Ic. The available gas-gas data for full scale coaxial elements 
is shown in Figure lb while the results of an ALRC sponsored single element 
program are shown in Figure Ic. This data indicates that the following design 
parameters are important for shear mixing elements. 


*For reader convenience, the figures appear at the end of each Roman Numeral 
section. 
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a. Relative area ratio 

b. Relative velocity ratio 

c. Absolute size 

d. Shear area between the two propellants 

e. Swirl 


Momentum Mixing Elements 


Since there is no available impinging mixing model which is 
computationally tractable, impinging elements must be evaluated empirically. 
Rupe and his coworkers at JPL as well as other investigators documented the 
results in the 1950 's of an orderly investigation into the impingement mixing 
process for liquid propellants using inert propellant simulants (ref. 1, 2, 3, 
4, and 5). By collecting immiscible fluids with a ganged probe, they were 
able to map mass and mixture ratio distributions for the common impinging 
elements (doublet, triplet, quadlet, and pentad) as a function of design and 
operating conditions. From his investigation, Rupe was able to correlate the 
maximum mixing potential for each of the elements using certain design vari- 
ables. These correlations took the following form: 



K and a are constants which must be determined empirically for each configura- 
tion. Rocketdyne used these correlations directly for gaseous propellants in 
Ref. 6 and concluded that the correlations were in fair agreement with the 
data. Aerojet's impinging premix (triplet) gas-gas cold flow data also was 
examined to determine applicability of the Rupe criterion. The results are 
shown on Figure 2 where it is apparent that the test data did correlate well 
with momentum ratio and geometry parameters. 

Based on this work it was determined that the following 
physical and hydraulic parameters are important. 

a. Relative stream momentum 

b. Relative orifice size 

c. Impingement angle 

It was also inferred that the following parameters are likely to be important: 

a. Physical size (thrust/element) 

b. Spacing 

c. Element orientation 
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C. BARRIER COOLING ANALYSIS 

An evaluation of H2 film cooling vs barrier cooling was conducted 
with the ground rule that the throat heat flux and wall temperature were to be 
maintained at 1960 watt/cm2 (12 Btu/in. 2-sec) and 533°K (500°F) wall, respec- 
tively, based on chamber life considerations. This analysis was conducted 
using the mixing program developed in LeRC Contract NAS 3-14343, (Ref. 7). 
Performance was computed by a mass weighted average Igp based on the barrier 
and free stream at their respective throat mixture ratios. 

Barrier coolant characteristics were calculated using the Aerojet 
Barrier Film Cooling Program (BARFC). The BARFC Program is a film-cooling 
model, where mixing of the film coolant with the mainstream combustion pro- 
ducts is calculated via an empirically adjusted entrainment model. BARFC is 
designed to calculate film temperature and mixture ratio profiles axially 
along the chamber contour. It considers the simultaneous effects of momentum, 
chemistry, and energy transport through the mixing layer. The profiles are 
coupled via empirically derived shape factors. 

The Attitude Control Thruster hardware was modeled with the 19 cm 
(7.5-in.) long chamber composed of a cylindrical section and a 22° half-angle 
cone. Total weight flow was fixed at 1.55 kg/sec (3.45 Ibm/sec) and an over- 
all mixture ratio (0/F) of 4.0 was maintained. Barrier or film coolant was 
injected from the injector face. 

Results are presented in Figure 3; no adjustment has been made 
other than use of the inferred entrainment factor multiplier. Wall temperature 
was taken to be 533°K (500°F) which is typical of throat values required to 
meet the originally specified cycle life (IO6 cycles). Heat fluxes of 1934 
and 1960 watt/cm2 (10 and 12 Btu/in .2-sec) , typical of those accompanying a 
533°K (500°F) wall, were considered. 

Figure 3a shows the heat flux at the throat vs the barrier flow/ 
total flow (injection point) ratio for three barrier injection mixture ratios. 
For a constant heat flux (viz., wall temperature in a regenerati vely cooled 
system), the barrier flow fraction must be increased as the mixture ratio 
increases to offset the higher combustion temperature. The heat transfer 
coefficient remains fairly constant over the range of mixture ratios and 
barrier temperatures predicted at the throat. Dashed lines on these figures 
represent extrapolations, necessitated because of convergence failure of the 
computer program. 

Figure 3b is a cross plot of the preceding data and shows the injec- 
tion barrier mixture ratio vs the barrier flow fraction for the two designated 
heat fluxes. The increased flow requirement with increasing mixture ratio 
is apparent. 
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Figure 3c is a plot of the barrier mixture ratio at the throat vs 
barrier flow fraction at the injection point, and it indicates that throat 
mixture ratio drops rapidly with increasing flow fraction. To achieve a given 
flux level, the barrier mixture ratio remains almost unchanged over a wide 
range of barrier flow fraction, reflecting the temperature sensitivity of flux 
since the heat transfer coefficient is approximately constant. The barrier 
flow fraction at the throat vs the injection barrier flow fraction is shown in 
Figure 3d for the three barrier injection mixture ratios. This figure also 
shows the two heat flux levels. Again, the dependence of heat flux on initial 
barrier mixture ratio is demonstrated. 

Based on this analysis, the decrement in performance was computed 
using a mass weighted two stream tube model using the throat barrier mixture 
ratio and barrier percentage weight flow as defined on Figures 3c and 3d for 
a throat heat flux of 1960 watt/ cm2 (12 Btu/i n. 2-sec ) . The core mixture ratio 
at the throat was then found by difference. The results of this performance 
analysis are shown on Figure 4. This analysis indicates that barrier cooling 
would penalize the engine system more than fuel film cooling [(Barrier 0/F) 

Ini = chamber life requirements and that the performance 

penalty increases as (Barrier 0/F)inj is increased. 

The effects of element characteristics was beyond the scope of 
this barrier cooling study, and therefore the results of this analysis did not 
influence directly the selection of element types for cold flow evaluation. 
However the analysis did provide a theoretical basis upon which to evaluate 
the full scale barrier cooled designs fabricated and tested in Task IV 
(See Section VII). 

D. STABILITY CONSIDERATIONS 

Stability characteristics of gas-gas systems are significantly 
different from liquid/liquid systems investigated to date. Recent experi- 
mental results at ALRC and work at Purdue University (Ref. 8) suggest that 
kinstics pldy the major roll in the stability of a gas- 9 ^^ system rather than 
the generally accepted causal mechanisms for liquid injectors such as pro- 
pellant mixing and vaporization. The influence of kinetics is based upon the 
observed sensitivity in the stability of a system to changes in mixture ratio. 
This effect is predicted with a model such as the Geode model (Ref. 8) which 
uses kinetics as the rate controlling process. The trends summarized in 
Figure 5 reflect a strong influence of mixture ratio and some effect of 
chamber pressure. Such trends in stability are apparent from the model of 
Reference 8 which gives the following expression for this energy release rate: 


Q' = iT 


RT ^ T 


( 2 ) 


where: 

(~j = steady-state quantity 
(') = perturbation quantity 
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From this expression, it can be seen that, increasing temperature 
and pressure, and consequently, density, has a stabilizing effect. This would 
imply that, as the stoichiometric mixture is approached (which generally gives 
a near maximum steady-state temperature), the system becomes more stable. 

Such a trend is shown on Figure 5. 

The significance of such a mechanism being the driver for unstable 
combustion is that the contustion process would not exhibit preferential fre- 
quency. This will make the use of conventional damping devices such as liners 
and baffles difficult because of their limited frequency range of effective- 
ness and could result in a major development hurdle for gas-gas systems. Also, 
with Kinetics playing such a key role in the stability area, it is difficult 
to select elements on the basis of stability characteristics. Therefore the 
primary criteria used for selecting the elements for cold flow evaluation were 
based primarily on the Concept Screening and Element Evaluation analyses dis- 
cussed previously. 

E. ELEMENT SELECTION 

Based on the above analyses the element configurations listed in 
Figure 6 were synthesized and selected as the most promising concepts. The 
element geometry variations selected for cold flow evaluation are also included 
in Figure 6. These geometry variations were influenced by the full scale 
operating ranges specified for the Space Shuttle Attitude Control Thruster, 
which are tabulated below: 

Thrust 

Chamber Pressure 

Mixture Ratio 

Propellant Inlet 

Temperature 

Propel 1 ants 


1790 - 8900 N (500-2000 lb) 
68 - 341 N/cm^ (100-500 psia) 
3-5 

167 - 333 °K (300-600°R) 
O2/H2 


F. COLD FLOW HARDWARE DESIGN 

1 , Singl e Element Injector-Chamber 

The single element hardware assembly is depicted in Figure 7, 
and the selected element designs are shown in Figure 8. Element configura- 
tions were changed by replacing the element plate and the oxidizer orifice 
plugs. Part Nos. 105 and 106, on Figure 7. The element plates shown in 
Figure 8 were removed and changed by unbolting the split ring (Part No. 103) 
and lifting out the remaining assembly, which then allowed access to the 
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element plate. The oxidizer plug access was through Part No. 101. The plug 
was removed and reinstalled by removing Part No. 101 from the assembly, 
extracting the tested plug, inserting a new seal, installing the new plug, 
and reassembling Part No. 101 to the test hardware. Hydrogen and nitrogen 
(the oxygen cold flow simulant) were introduced to the test article through 
Part Nos. 102 and 101, respectively. 

A helium bleed circuit was incorporated into Part No. 104, 
which permitted helium to be injected axially through an annular Rigimesh 
plate that covered the injector face from the element to the chamber wall. 

This helium circuit was designed to suppress face recirculation using the 
criteria described in Reference 9. A recirculation parameter, Cj, is defined 
in Reference ^^h-jch is a function of the source jet velocity, the field 
velocity, and the radius ratio of the jet to the chamber. A jet velocity of 
91.5 m/sec (300 fps) a radius ratio of 0.17, and a field velocity of 7.61 m/sec 
(25 fps) suppressed face recirculation using the Cj criterion defined in 
Reference 9. The addition of a third fluid circuit does not compromise the 
resulting test data as the mass spectrometer used to probe the flow field has 
the capability of determining the mole fraction of any and all species. 

2 . Multiple Element Injector-Chamber 

The design of this hardware was similar to that of the single 
element hardware, with the same consideration given to versatility of element 
configurations and test setup. The assembly of this hardware is shown in 
Figure 9 together with the elements and face geometry. The basic purpose of 
the multiple element test rig was to permit evaluation of inter element mixing 
effects. This was accomplished in two ways. (1) The radial spacing of the 
elements could be varied (three positions). The elements during testing were 
arranged in a square injection pattern as shown in Figure 9. During testing, 
the active elements were located in one of the three hole circles. The remain- 
ing holes were blocked with blank elements. (2) The elements were rotated 
with respect to each other such that the fuel -rich zones of one element were 
directed toward the oxidizer-rich zones of the adjacent element. Based on the 
single element cold flow work it was determined that the coaxial, swirl coaxial 
and premix triplet elements would be investigated. The 67N, (15 Ibf) elements 
are shown in detail in Figure 9. The helium face bleed shown in Figure 9 was 
not operational. It was eliminated because the test facility did not have the 
flow capacity to maintain the required He flow rate. 

3 . Sampl ing Rake and Rake Assembly 

The sampling rake shown in Figure 10 was used to acquire mass 
and pressure data in the flow field downstream of the injector face. It was 
a multiple entry type probe with 20 active positions. Each entry port performed 
two functions: (1) it provided a local gas sample which was drawn from the 

flow field and routed to a mass spectrometer and (2) it served as a total pres- 
sure probe when the test apparatus was sequenced to the pressure scan mode. 
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The rake was located and held by the pintle assently shown in 
Figure 7. The rake and instrumentation leads were routed through Part No. 406 
and sealed against chamber pressure by a Swage-lock fitting. Provision was 
made to traverse the rake both axially and circumferentially. In operation, 
chamber pressure was set by the remotely controlled, belt driven pintle that 
enabled the throat area to be varied until the required pressure was obtained. 


- 15 - 





'I ‘X0U.M3jJj:n ilU|XTW 


17 


Figure 2. Fffect of Momentum Ratio on Premix Impinging Element Mixing Efficiency 
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Figure 6. Element Types Selected for Cold Flow Evaluation 
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Figure 8. Single Element Cold Flow Designs 






IV. 


COLD FLOW MODELING 


A. SINGLE ELEMENT COLD FLOW 

Using the methods described in Appendix B some 223 single element 
cold flow tests were conducted and analyzed. The data from these tests is 
summarized in Table III, Table IV and in Figures 11 through 24. Table III, 
Single Element Cold Flow Test Summary is a tabulation of the injection ele- 
ment design parameters, test operating parameters, injection parameters and 
calculated mixing efficiencies (Em). Table IV is a summary of element design 
and operating effects on mixing efficiency. The figures present both mixing 
efficiencies and compatibility data at the chamber wall as a function of 
thrust/element, (F/E), length, element area ratio (Af/Ao) and mixture ratio. 
Compatibility at the chamber wall is indicated by the measured local mixture 
ratio at the edge of the flow field (0/F)wall divided by the overall nominal 
mixture ratio (0/F)nominal • The open symbols represent the most oxidizer-rich 
composition measurecl at the flow field boundary, the closed symbols represent 
most fuel-rich composition and the half-open symbols represent the average 
boundary composition. 

Figure 11 is a summary graph of E^ for all injectors as a function 
of L/D, F/E, Af/Ao 0/F. The effect of sample position on L/D, where D is 
the oxidizer orifice diameter or equivalent diameter if the element is non- 
circular, is shown in Figure 11 for the area ratio one elements. The char- 
acteristics of the elements are as expected, increasing Em with increasing 
length. The swirl coaxial, premix pentad, and triplet elements group together 
at high mixing efficiencies, while the coaxial element has the lowest element 
Em. The remaining elements group between the extremes of the shear coaxial 
and the swirl coaxial elements. When these elements are compared on a thrust/ 
element basis, the same effects are noted i.e,, the premix designs and the 
swirl coaxial element are high mixing rate elements while the shear coaxial 
is a low mixing rate element. With the exception of the swirl coaxial element, 
the data trends are in the expected direction, reduced mixing with larger ele- 
ments. The swirl coaxial element Em is maximum at 89 N/Element (20 Ib/element). 
It appears that the fuel annulus size and radial oxidizer momentum are related. 
It is interesting to note that the parallel sheet element at low F/E 3 Ibf 
(13 Newtons) approaches the mixing efficiency of the best elements. For small 
thrust/element injectors (such as HIPERTHIN), these data indicate that the 
parallel sheet element is an excellent element concept. 

Area ratio influences for the basic element concepts are illu- 
strated in Figure 11 as well. Both shear mixing elements, the coaxial and 
increased shear element, decrease in performance as fuel to oxidizer orifice 
area ratio is increased. It should be noted that the area ratio was varied 
by changing only the fuel geometry; the oxidizer diameter was held constant. 
Decreasing area ratio is in the direction of both smaller fuel annulus widths 
and higher velocity ratios or delta velocities; effects which analytically are 
predicted to lead toward increased mixing efficiency. 
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TABLE III 

SINGLE ELEMENT GOLD FLOW TEST SUMMARY 



Shear 13.3 (3) 0.49 536 537 11.63 41.32 



SINGLE ELEMENT GOLD FLOW TEST SUMMARY 
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( 289 ) ( 287 ) ( 31 . 90 ) ( 20 . 59 ) 




TABLE III Sheet 3 of 

SINGLE ELEMENT GOLD FLOW TEST SUMMARY 




SINGLE ELEMENT COLD FLOW TEST 



143 X 0*^ I 507 523 10. 18 7.35 .92 . 58 3.64 13 3.95 . 7 50 4.75 95.5 91.6 '^9,88 

144 T 90° G.C.W Jr A (282) (292) (7,02) (5.06) 15 



SINGLE ELEMENT COLD FLOW TEST 



31 



SINGLE ET.RMENT GOLD FLOW TEST SUMMARY 




SINGLE ELEMENT COLD FLOW TEST 



142 210 I I II 484 530 9.56 22.40 2.13 .60 8,10 4.9 3.80 .427 5.22 84.0 75.9 07.70 



SINGLE ELEMENT COLD FLOW TEST 



0) C 
0) 


34 


t x ▼ T ^09 530 36.52 33.91 .93 .93 3.68 



SINGLE ELEMENT COLD FLOW TEST SUMMARY 




OF DESIGN AND OPERATING EFFECTS ON MIXING EFFICIENCY 
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Indicates no influence 




The swirl coaxial and premix designs are insensitive to area ratio 
influences, an effect that was somewhat surprising for the premix case. Cold 
flow testing by ALRC on an element design similar to the premix triplet design 
revealed that mixing efficiency was significantly impacted by orifice geonretry 
changes i see Reference 10 for these data. Although the single-element injec- 
tors of both these series were similar, two important geometry differences are 
evident. The mixing cup length of the element tested on this program was 
significantly longer than that of Reference 10. For this program, cup depths 
ranged from 1.15 to 4.75 L/D, where L is the cup depth and D is the oxidizer 
hole diameter. For the program of Reference 10, the ratio, on the same basis, 
was 0.1. The increased mixing cup depth forces mixing and therefore attenuates 
the effects of area ratio and mixture ratio (momentum ratio) variations. In 
addition this program only investigated rectangular fuel slots while the data 
of Ref. 10 was obtained with a fuel slot geometry that was similar to a capital 
"I" in cross section. 

Referring again to Figure 11, note the very sensitive external 
impinging triplet orifice area ratio influences, where the element optimized 
at an Af/An of 1.25. This element type, in general, is sensitive to most 
geometry and hydraulic variations and a designer must pay particular attention 
to the design parameters in order to obtain an optimum design. 

The effect of mixture ratio on is illustrated in Figure 11 as 
well. Again, the swirl coaxial and premix designs are relatively insensitive 
to this influence, while the shear coaxial element decreases in efficiency as 
a function of increasing mixture ratio (less aV). The remaining elements were 
found to be insensitive to variations in mixture ratio. 

Detailed parametric curves of data obtained from the coaxial family 
of elements are illustrated in Figures 12 through 15. Included in this family 
are the coaxial, increased shear coaxial, and swirl coaxial ele.iients. Figure 12 
compares these elements on a length basis (L/D) (D is the oxidizer tube diam- 
eter), a mixture ratio basis, and a propellant temperature basis. The effect 
of increased shear area and swirl is to increase mixing efficiency. Increased 
surface area does promote increased mixing as illustrated by the increased ^ 
shear element Em data. Radial transport of oxidizer by swirling forces mixing 
when compared to simple coaxial mixing element. When compared on a mixture 
ratio basis, the swirl coaxial shows no influence and the shear and increased 
shear coaxial exhibit decreasing efficiency with increasing mixture ratio. 

This effect is due to a decreased aV influence. The effect of lower propellant 
temperature on Em is shown in Figure 12 as well. With the exception of the 
swirler data, Em increases with decreasing temperature. This trend is contrary 
to hot-fi rG data ganeratod with 1500 Ibf thrusters on Contract NAS 3-14352, 

(Ref. 11), which indicated a slight reduction in performance when tested with 
cold propellants. The contrary temperature effect noted in hot testing may be 
influenced by the contustion influence on the mixing/ reaction process. 
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The effects of area ratio (or fuel annulus width) for these ele- 
ments are illustrated in Figure 13 as a function of thrust/element. In general, 
decreasing fuel to oxidizer orifice area ratio is in the direction of increased 
Em for the three coaxial type elements; again indicating the influence of aV. 

The swirl coaxial data in Figure 13 indicate an optimum in thrust size at 
approximately 67 Newtons/element (15 Ibf/element) . This indicates that the 
radial component of momentum and the physical size are related, or that there 
is a fuel momentum to oxidizer radial momentum ratio which produces optimum 
mixing. The effect of swirl magnitude for various area ratio elements is 
illustrated in Figure 14. These data indicate that large changes in mixing 
efficiency can be obtained with a rather moderate amount of oxidizer radial 
momentum, Vt/Vg = 0.5). 

The compatibility characteristics of the three shear elements are 
indicated in Figure 15, where the coaxial and swirl coaxial elements are com- 
pared in Figure 15a. In the near zone (close to the injector face) the coaxial 
element is fuel rich as compared to the swirler. This effect was noted in 
combustion testing of full-scale hardware in Contract NAS 3-14354 where the 
near zone heat flux using a swirler element was higher than the equivalent 
coaxial element. Ref. 12. Combustion data generated during this contract, dis- 
cussed in a later section, also supports this observation. As a function of 
Af/Ao, the coaxial element becomes more fuel -rich with increasing area ratio 
(see Figure 15c). This fuel rich boundary is reflected in the low Em data for 
this element, which was previously discussed. It is interesting to note how- 
ever that as the Em improves (lower Af/Ao) the compatibility decreases. In 
general, this characteristic is universal, i.e., those elements that exhibit 
high potential compatibility do so at the expense of a low mixing rate. 

An equivalent data set for the premix designs is illustrated in 
Figures 16 through 18. With the exception of length, shown on Figure 16, 
premix elements are relatively insensitive to design and operating variables. 

For all variables investigated, resulting Em was approximately 90% or higher. 
Referring to Figures 16 and 17, note a decrease in Em with increasing thrust/ 
element for both the pentad and triplet; a slight decrease in Em with tempera- 
ture and a small effect of mixture ratio on Em. However, all these influences 
are minor, which indicates that premix elements can run over a wide range of 
engine conditions without suffering a serious performance penalty. As dis- 
cussed earlier, however, the fact that they are insensitive is somewhat sur- 
prising and can be attributed to the fact that the mixing cup depths used on 
this program attenuate the expected response to variations in design and 
operating conditions. 

The premix element compatibility variations, illustrated in 
Figure 18, are not negligible as demonstrated by the large differences in 
maximum and minimum local 0/F. This condition indicated a potential chamber 
wall streaking problem. Note, for the premix triplet. Figure 18b, that 
increasing the element width ratio tended to stratify the flow field. This 
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can be explained by considering the geometry of the injector pattern. For a 
narrow fuel slot width, the available fuel momentum is concentrated on the 
center of the oxidizer stream, penetrating the bulk of the oxidizer resulting 
in a homogenous flow field. As the penetration is reduced by broadening the 
fuel slot, the fuel momentum is insufficient at the center of the oxidizer 
stream. The net effect is to form a flow field edge condition that is alter- 
nately fuel and oxidizer rich at tt/ 2 rad (90°) intervals. The pentad ele- 
ment does not appear to streak as readily because the available momentum of 
any one fuel stream is only half of the equivalent premix triplet element. 

In no case did the pentad completely penetrate the oxidizer stream and, in 
general, always tended to result in unmixed flow rather than completely pene- 
trate and mix with the oxidizer. The streaking characteristics of the 
Af/Ao = 1.0 triplet element can be attenuated by reducing the total fuel 
momentum with larger fuel orifice area, as illustrated in Figure 18c. 

Data similar to those presented earlier for the premix and shear 
designs are found for the like doublet element in Figures 19 and 20. Parti- 
cular attention should be directed at Figure 19 where the element spacing 
effects are parametrically illustrated. These data indicate that an optimum 
is found at a fuel to oxidizer spacing of approximately 0.25 cm (0.1 in.). 
Element spacing lower than the optimum have lower mixing efficiency because 
the primary fans interact strongly forming secondary fans. These fans are 
nonhofnogenous in nature, with fuel on one side of the fan and oxidizer on the 
other. When the doublet spacing is larger than the optimum value, complete 
mixing is prevented because the elements are too far apart. Compatibility 
of this element, as shown in Figure 20 was not good since large circumferential 
variations in 0/F or streaks were measured. This result was caused by testing 
a single pair of oxidizer and fuel elements. Multi-element injectors should 
attenuate some of this unmixed flow. 

Data for the parallel sheet elements are found in Figures 21 and 22. 
Mixing efficiency as a function of length, thrust/element, and propellant tem- 
perature is shown in Figure 21, compatibility in Figure 22. Like the other 
elements investigated, mixing efficiency increases with length and decreases 
with F/E. As a function of 0/F, the data indicates an increase in mixing 
efficiency, contrary to the trend of the similar coaxial element. Compati- 
bility effects are similar to those observed with the doublet element; wide 
0/F variations with (0/F)wall depending on which side the element is probed. 
Again, some attenuation of this effect would be expected in a multiple element 
configuration . 

The external impinging F-O-F triplet data are depicted in Figures 23 
and 24. In general, this element configuration is very sensitive to all oper- 
ating and geometry variables. Figure 23 displays mixing efficiencies for 
various impingement angles and area ratios as a function of thrust/element. 

At the 67 N (15 Ibf) level, the extreme sensitivity of mixing efficiency with 
impingement angle is noted. One also notes that area ratio influences are 
optimized in the 1.0 to 1.25 A^/Aq range. Referring to Figure 23b note that 
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the Af/Ag = 1.0 mixing efficiency increases dramatically with impingement angle. 
Inspection of the detailed flow field reveals that, as the radial momentum com- 
ponent is increased by increasing the impingement angle, increased penetration 
of the axial oxidizer jet occurs, thus yielding a more homogenous flow field. 
However, for the area ratio = 0.50 data, a contrary trend is noted. In this 
case, overpenetration is achieved and increasing the impingement angle only 
amplifies the excess penetration of the fuel. Mixing was also determined to be 
a strong function of mixture ratio coupled with impingement angle as shown in 
Figure 23d. On the fuel-rich side of the optimum, overpenetration occurred; 
on the oxidizer- rich side, under-penetration occurred. 

In general, the compatibility of this element was expected to be 
poor. Figure 24 indicates the degree and magnitude of the potential compati- 
bility problem for the overpenetrated oxidizer jet. The entire flow field 
boundary is higher in O/F than the nominal with a large composition variance 
in the circumferential direction. 

B. MULTIPLE ELEMENT COLD FLOW 

As an integral part of the cold flow program, three prototype 
injection elements were flowed in a multiple element test configuration (four 
elements). The intent of this testing was to generate information that would 
allow evaluation of interelement mixing effects. The three injection elements 
were 67 Newton (15 Ibf) thrust shear coaxial, swirl coaxial, and premix triplet 
configurations. These elements are illustrated in Figure 9. They are arranged 
in a square injection matrix as shown in that figure. Testing was conducted 
as outlined in Appendix B and was essentially identical to the technique used 
in single element testing. The reduced data and test conditions are tabulated 
by test number in Table V. The first sixteen tests were conducted with the 
shear coaxial element, tests 17 to 32 were with the swirl coaxial, and the 
remainder of the tests were with the premix triplet element. Included in this 
table are both the single element and multiple element mixing efficiencies. 

The single element and multiple element data for the shear coaxial 
element are compared in Figure 25. The data in this figure represent identical 
test points for both the single and multiple element testing. Inspection of 
Figure 25b, c, and d indicates that a dramatic increase in mixing efficiency 
occurred when multiple elements were used. However, two observations run 
counter to this conclusion. The first is shown in Figure 25 where the data 
show that element spacing influenced E^ only slightly during the multiple 
element testing. In fact, as the element-to-element distance is increased 
(from Row A to Row C), the apparent multiple element mixing efficiency 
increased. If element proximity tends to increase mixing, then an increase 
in element-to-element distance should decrease the relative mixing efficiency - 
an effect not noted. The second observation was that a large amount of recir- 
culation occurred between the elements. Flow apparently was coming from the 
far downstream region. In this region, the flow is well mixed which, when 
recirculated to the near field and sampled, would produce an apparent increase 
in mixing efficiency. 
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GAS/GAS MULTIPLE ELEMENT COLD FLOW DATA SUMMARY 
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(1) Tangential Fuel Slots 



Some insight into this recirculation mechanism was gained by look- 
ing at the element centerline oxidizer concentration data. This was done by 
inputting, the multiple element design parameters and operating conditions 
into the single element cold flow mixing computer program. It was found that 
the single element model provided an excellent prediction of the multiple 
element fuel and oxidizer element centerline mass fraction composition 
(mixture ratio) and axial mixing distributions. This implied that multiple 
element mixing rates were much closer to the single element mixing rates than 
the higher experimental E 's indicated. 


In addition, it was noticed that the multiple element test inte- 
grated mass factors were significantly higher than unity, indicating gross 
recirculation occurring between the multiple elements. To adjust the data 
for this influence, it was assumed that the recirculation gases flowing between 
the elements would come from downstream and consist of essentially uniformly 
mixed gases (E^ = 100%). Then by the definition of E , the corrected E is: 


E , ME = 100 
m 


n W. r 0/F - 0/F. " 
^ [o/F (1+0/F.), 





0/F. - 0/F 

^Twf- 


(3) 


it can be shown that the corrected 


E , ME = 100 
m’ 


'^int 


[100 - E^, ME, Meas] 


(4) 


Results from this data re-evaluation for the multiple shear 
coaxial elements test are summarized in Table VI. The measured multiple ele- 
ment mixing efficiencies are shown for each of the test conditions evaluated. 
The "mass integration" column in Table VI E^ass* represents the ratio of the 
integrated mass flowrate over the entire flow field divided by the injected 
mass flowrates. This factor is approximately 4 except for two tests. Thus, 
a constant factor of 4 was used as the mass integration factor for all tests 
to determine the corrected multiple element mixing efficiency. Where the data 
is available, the comparable single element mixing efficiency (Ep^, SE) also is 
tabulated in Table VI. A comparison of both raw and corrected multiple ele- 
ment vs single element Em's are shown graphically in Figure 26. The Em, ME 
vs Ejp, SE is the data previously mentioned in Table V and Figure 25. When 
corrected for the mass integration (recirculation) factor the multiple ele- 
ment mixing efficiency is comparable to the single element data. Thus, it was 
concluded that when recirculation effects are compensated for there is no 
significant effect of multiple element interactions in comparison to single 
element mixing. 
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TABLE VI 


MULTIPLE SHEAR CO-AX ELEMENT COLD FLOW SUMMARY 
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The Effect of Length on 
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Figure 22. Compatibility of Parallel Sheet Element 
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Mixing Efficiency, Single Element 


V. 


HOT FLOW MODELING 


It is ALRC's belief that cold flow testing alone does not completely 
describe the mixing process under combustion conditions. Chemical reaction 
does impact adversely the mixing profiles generated in cold flow testing. 

In Reference 10, this impact was described both analytically and empirically 
for a conventional shear coaxial element. In general, the combustion adversely 
biased the mixture ratio profiles, i.e., mixing profiles measured in cold flow 
were displaced downstream axially by cotifcustion. As part of this program, the 
effects of combustion were studied for a high mixing efficiency element, the 
swirl er coaxial element. 

The intent of this portion of the program was not to model the combus- 
tion process, but rather to gain phenomenological information that would 
broaden the insight into the combustion process. In order to accomplish this 
end, a limited test program was designed that was structured to look at the 
differences between "identical" combusting and noncombusting flow fields. 

This testing was conducted with the cold flow swirler element at ambient 
pressure maintaining Mach number similarity in both the injector and chamber. 

The data generated was local mixture ratio distributed axially and 
radially through the chamber. Both hot and cold flow gas flow fields were 
generated with identical hardware and flow conditions. The only difference 
between these tests was that in one case the H 2 /O 2 propellants were burning 
and in the other case they were not. The hot flow hardware was the same hard- 
ware used during the cold flow program. 

★ 

This hot fire testing was conducted in a steam-cooled chamber and 
utilized a double-walled, hot H 2 cooled probe to withdraw the gas sample. The 
sample was routed into the mass spectrometer via copper steam-heated tubing. 

The resulting data conclusively indicated a combustion mixture ratio profile 
substantially different than the cold flow profile. These data are presented 
in Figure 27 as a function of the probe radial and axial distance. This figure 
indicates that the cold flow mixture ratio gradients had essentially dissipated 
at the first axial position, indicative of the desirable mixing qualities of 
the swirl coaxial element. Comparing the hot flow profiles with the cold flow 
profiles at the same length, it is apparent that the mixture ratio profiles 
are steeper (greater deviations from uniformity) and persist to 18.2 cm 
(7.17 in.). This data comparison substantiates that the combustion does 
adversely affect the mixing. These data are plotted in a different format 
in Figure 28 where three mixing zones are postulated. 


♦Particular care must be taken in this type of testing to ensure that the 
combustion products (H 2 O) do not condense anywhere in the system. Thus, 
all components were heated to over 250° F during testing using either steam 
or heated hydrogen. 
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1. The near zone where high mass transfer rates of reactants 
exist and combustion is initiated. 

2. An intermediate zone where the large density gradients due 
to combustion products displace the oxidizer inward and the 
fuel outward. 

3. A region where the continued rotation of the oxidizer-rich 
core promotes additional radial diffusion of the oxidizer- 
rich products into the fuel-rich outer stream tubes. 

Based on this data the scaleability of the cold flow data directly to combus- 
tion conditions must be questioned. In later sections of this report the 
effects of contustion are again investigated and a relationship is determined 
that corrects the cold flow data to compensate for hot fire conditions. 


- 61 - 



lUD ‘U0LSU9U1LQ L6!-P®y 


Ln o Ln 






- 62 - 


Figure 27. The Comparison of Hot and Cold Mixture Ratio Profiles 
for a Swirl Coaxial Element 





VI. SINGLE ELEMENT INJECTOR EVALUATION 
A. TEST ARTICLES 

Based on the cold flow work described previously, four 222N 
(50 Ibf) single element injector configurations were selected as being poten- 
tially high performing or having good compatibility characteristics. These 
injectors were: 

1. Coaxial element 

2. Swirl er coaxial elements 

3. Premix element 

4. F-O-F triplet elements 

The injectors are illustrated in Figure 29 together with the salient design 
parameters. The coaxial and swirl coaxial units were the same hardware used 
during the cold flow program, both high (V^/Vg = 1.0) and low (Vt/Vg = 0.5) 
swirl velocity elements were investigated. The only difference between these 
elements and the cold flow elements was a seal surface added to the hot fire 
injectors. The premix and triplet injectors were also fabricated from the 
same components used during cold flow with the exception that new faces were 
manufactured from OFHC copper rather than the stainless steel used for cold 
flow. 


Three 222N (50-lbf) -thrust chambers were built to interface with 
the cold flow injector assembly shown in Figure 7. These chambers together 
with the pertinent physical dimensions are illustrated in Figure 30. Two 
chambers were 5.1 cm (2.0 in.) long with a 2.5 and 3.5 contraction ratio, 
Ac/At, and referred to as the -1 and -3 configurations, respectively. The 
third chamber was 14 cm (5.5 in.) long with a 2.5 contraction (the -2 con- 
figuration). All the chambers had an exit area ratio, Ae/At = 1.75. They 
were copper heat sink in type and were designed to have a 2-sec test duration 
and provide adequate thermal response to obtain heat flux data. Twelve 
thermocouples were buried at approximately mid-wall depth. The thermocouples 
were located 1.72 cm (0.5 in.) apart axially and in three rows, tt/ 4 rad (45 ) 
apart circumferentially. Provision was made for a spark type igniter to be 
located in the chamber wall. The two Pc taps were located near the injector 
face and at the entrance to the convergent section of the nozzle. 

B. TEST FACILITY AND PROCEDURES 

The uni -element testing was conducted in the Physics Laboratory, 
Test Bay 1; a bay that has been used extensively for H 2 /O 2 ignition and uni- 
element contiustion research. The facilities in Test Bay T include gaseous 
hydrogen and gaseous oxygen feed systems of 0.154 m3 (5.5 ft3) capacity each, 
an alternate hydrogen and oxygen feed system of smaller capacity, and complete 
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instrumentation and process control capabilities. These systems are enclosed 
within insulated cylinders which provided temperature conditioned propel- 
lants. The propellants were conditioned by a LN 2 spray over the vessels for 
subambient tests. Figure 31 is a schematic flow diagram of the propellant 
feed systems. 

For ambient propellant temperature tests, system gas regulators 
were used for setting the supply pressure to critical flow nozzles. For the 
temperature conditioned tests, the feed systems were operated in a blowdown 
mode, and critical flow nozzles were selected to provide the required flow 
at feed system averaae pressure. 

Sy-tcm .niil (lirusi w«'ri‘ iiK'.iMircd ir,iu<i .ix-win' im'.iMirc 

iirmiL ‘.y.lfm iiLili/iiu| '.I tain iMV'.'.mo or loitr I |•.ln^lllu m . Iln* Iran-, 

ducers have 3‘jO Ohm strain gages in a tully <Ktiv«‘ lour arm l)i iili|i‘ ronliiiura 
tion. Typical transducers employed in this system are tlie laber Wodel t’Ui) 
pressure transducer or the BLIi Model U3XXA load cell. 

C. SINGLE ELEMENT TEST RESULTS 

The following paragraphs summarize the resulting performance, heat 
transfer and stability data from the 76 single element tests. Data is pre- 
sented that was obtained using the techniques described in Appendices C and D. 
Performance data is listed in Table VII and VIII in English and SI Units 
respectively, and is plotted in Figure 32 through 35. Heat transfer informa- 
tion is plotted in Figure 36 through 42 and stability data is presented in 
Table IX and plotted in Figure 43. 

1 . Performance 

A summary of the single element hot fire test series per- 
formance is presented in Tables VII and VIII and illustrated in Figure 32 
through 35. Identified for each data point are the applicable test condi- 
tions, the measured and theoretical performance values, specific impulse per- 
formance losses, and calculated performance efficiencies. Included under the 
test conditions are the test numijer, injector and chamber identification, the 
data period, and the chamber pressure, mixture ratio, and the oxygen and 
hydrogen manifold temperature. Some of the ERE data were computed to be oyer 
100% for the triplet element. This inaccuracy is believed to be caused pri- 
marily by the very large boundary layer heat loss calculated for these small 
copper heat sink chambers. Under some test conditions, heat losses were as 
high as 13% of delivered Igp. Also accuracy of measuring thrust values in 
the 222 N (50 Ibf) range is difficult to achieve, although the thrust cell 
was quoted to be calibrated within + 0.5% thrust accuracy. These two effects^ 
combined limit the accuracy of the single element data to be no better than 3% 
on an absolute basis. However, when comparing element performance data for 
different elements at a given test condition, the aERE values should be 
accurate within <1%. 
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TABLE IX 


STABILITY DATA 
FOR SINGLE ELEMENT INJECTORS 


Test 

No. 

Injector 
- Type 

Chamber 

MR 

c» 

psia 

P ( N / cm ) 
c 

C'-P/P . 

c ) ox 

05P/P )p 

c fuel 

Stability 
Freq, Hz 

103 

Swlrler 


.3 

3.83 

264 

182 

0.05 9 

0.119 

US 1050 

104 

(V /V -0.5) 



1,94 

283 

195 

0.026 

0.227 

US 1050 

105B 

c 

a 



5.86 

268 

185 

0.037 

0.040 

s 

106 





7.70 

252 

173 

0.044 

0.023 

s 

107 





0.99 

295 

203 

0.005 

0.438 

US 1050 

109 





3.96 

442 

305 

0.057 

0.110 



110 



> 


4,01 

87 

60 

0.037 

0.098 

US 1050 

112 



-1 

5.90 

275 

189 

0.129 

0.037 

s 

113 





1.99 

283 

195 

0.188 

0.256 

US 1500 

114 





6.03 

89 

61 

o 

o 

0,031 

S 

115 

Coaxial 



3.98 

240 

165 

0.008 

0.098 

s 

116 





2.03 

274 

189 

0.005 

0.200 

MS 750 

117 





5.93 

217 

149 

0.019 

0.062 

S 

118 





6,23 

69 

48 

0,026 

0.056 

s 

119 





1.94 

93 

64 

0.005 

0.208 

MS 750 

120 




t 

1.98 

447 

308 

0.004 

0.207 

MS 750 

121 




r 

3.89 

374 

258 

0.009 

0.103 

S 

122 




■ 2 

3.90 

315 

217 

0.007 

0,057 

MS 600 

123 




1 

2.03 

302 

208 

0.005 

0.199 

US 750 

124 




L 

6.01 

302 

208 

o 

o 

0.032 

S 

125 




Y 

6.00 

96 

66 

0.015 

0,028 

s 

126 



- 

3 

5.99 

207 

143 

0.023 

0.078 

s 

127 





7.73 

194 

134 

0.031 

0.061 

s 

128 





4.05 

73 

50 

0.030 

0.137 

MS 7 50 

129 





3.93 

222 

153 

0.031 

0.137 

MS 900 

130 





1.98 

251 

173 

0.028 

0.318 

US 900 

L3l 





3.91 

369 

254 

0,027 

0.133 

MS 750 

132 





0.99 

300 

207 

0.003 

0.437 

US 800 

134 





3.90 

311 

215 

0.226 

0.288 

S 

135 

Preraix 



6.01 

291 

201 

0.182 

0.208 

S 

138 





7.63 

258 

178 

0,181 

0,194 

s 

139 





2.05 

325 

224 

0.410 

0,581 

s 

140 





3.97 

522 

360 

0.227 

0.291 

s 

141 




f 

4.02 

103 

71 

0.217 

0.278 

s 

143 





0.99 

302 

208 

0.707 

1.159 

s 

144 



- 

2 

3.93 

303 

209 

0.226 

0.291 

s 

145 



1 

1 

2.03 

304 

210 

0.438 

0,625 

s 

146 





6.01 

291 

201 

0.166 

0. 1906 

s 

147 





4.12 

98 

68 

0.221 

0.279 

s 

148 



- 

1 

3.94 

317 

218 

0.221 

0.283 

s 

149 





2.06 

329 

229 

0.398 

0.569 

s 

150 





3.90 

524 

361 

0.218 

0.283 

s 

151 





5.99 

290 

200 

0.176 

0.203 

s 

152 



> 

f 

4. 14 

103 

71 

0.213 

0,270 

S 

153 



- 

3 

6.08 

91 

63 

0.185 

0.208 

s 

154 





1.96 

105 

73 

0.420 

0.600 

s 

155 

High 



3.89 

284 

196 

0.077 

0,077 

MS 700 

156 

Swlrler 


1 

2.02 

299 

206 

0.055 

0.191 

MS 1100 

157 

(V /V -1.0) 

> 


5.93 

275 

190 

0.099 

0.042 

s 

158 



- 

2 

5.93 

271 

187 

0.100 

0.048 

s 

159 



1 


3.92 

304 

210 

0.065 

0.064 

s 

160 



'I 

/ 

2.00 

307 

211 

0.043 

0.161 

MS 1100 

101 

-1 Triplet 

- 

3 

3.86 

262 

180 

0.273 

0.083 

US 3500 

102 





2.00 

312 

215 

0.L22 

0.128 

S 

103 





4.08 

84 

58 

0.314 

0.084 

US 3500 

104 





5.94 

233 

161 

0.395 

0.072 

US 3500 

105 

-2 Triplet 



5.92 

233 

160 

0.414 

0.073 

US 3500 

106 





4.15 

83 

57 

0.354 

0.077 

US 3500 

107 





3.92 

261 

180 

0.271 

0.078 

US 3500 

108 


f 



2.00 

319 

220 

0.116 

0.137 

S 

109 

-3 Triplet 



3.88 

305 

210 

0.197 

0.255 

s 

110 





1.93 

323 

223 

0.098 

0.594 

s 

lU 





4,11 

99 

68 

0.224 

0.254 

s 

112 



> 

f 

5.77 

255 

176 

0.324 

0.154 

s 

113 




2 

5.87 

286 

197 

0,247 

0.141 

S 

114 



1 


4,10 

98 

68 

0.178 

0.232 

s 

115 



1 

/ 

3.92 

300 

207 



S 

116 

Triplet 


2 

1.98 

300 

207 

0.114 

0.644 

s 


(1) US - Unstable 


S - Stable 

MS - Marginally stable 
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Figure 32 shows the variation of energy release efficiency 
with mixture ratio for the five injector elements tested and the two 5.08 cm 
(2.0-in.) long combustion chambers (-1 and -3). These data show that the 
F-O-F triplet element is the highest performing at low mixture ratios, while 
the premix injector is the highest consistent performing element over the 
total mixture ratio range. The high swirler coaxial element, the low velocity 
swirler element and the conventional coaxial element are in order of decreasing 
performance. The effect of mixture ratio and contraction ratio on ERE is also 
shown to be a function of the injector element design. The premix triplet ele- 
ment ERE increases with decreasing mixture ratio over the entire test mixture 
ratio range from 8 to 1 . The swirler coaxial element, on the other hand, 
reaches a minimum ERE near the nominal design mixture ratio of 4 and increases 
from this value for either increasing or decreasing mixture ratio. The coaxial 
element achieves maximum ERE at a mixture ratio of 1 and decreases significantly 
('v-25%) with increasing mixture ratio to a minimum value at a mixture ratio of 6. 
Further increases in mixture ratio result in slightly (-^2%) higher energy 
release efficiency. 

The detailed F-O-F triplet performance is illustrated in 
Figure 33. Energy release efficiencies for the three injectors and two cham- 
bers are represented as a function of mixture ratio. Two characteristics are 
apparent from the data in Figure 33; (1) all injectors are higher performing 

at low mixture ratios; (2) the low area ratio element (Af/Ao = 1.0) is sub- 
stantially higher performing than either of the Af/Ao = 2 elements. Both the 
above trends indicate that this element is sensitive to fuel momentum; i.e., 
there must be sufficient fuel momentum to penetrate the oxidizer jet to effect 
complete propellant mixing. 

A correlation of these data as a function of momentum ratio 
is depicted in Figure 33 as well. The stable data from all the injectors 
correlate well with the mixing parameter, MP. 


MP 


0 0 

w^ sin a 


(4) 


Referring to Figure 33a, note that the mixing parameter maxi- 
mizes at a value of 0.5 in the short chambers and is flat up to a value of 2 
in the long chambers. However, the value of 0.5 is not realistic for practical 
engine systems, since low values of the mixing parameter are in the direction 
of high fuel pressure drop. An O 2 /H 2 engine operating at mixture ratio 4.0 
with an oxidizer pressure drop of 10% of chamber pressure would have a sonic 
fuel injection velocity if a MP of 0.5 were maintained. Depending on the 
particular injector design, the fuel pressure drop would be on the order of 
the chamber pressure, a value too high to be realistically considered for 
practical engine systems. Since the mixing parameter tends to be rather flat 
over the range of 0.5 to 2 for the longer chamber, it is therefore not an 
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adverse compromise to relax the criterion of MP = 0.5 to a value somewhat 
higher. The performance penalty should not be large, and the fuel pressure 
drop is more consistent with practical injector design practice. 


The effect of chamber contraction ratio on ERE is seen in 
Figure 32 to be a function of the injector design. For the premix triplet 
element, ERE increases with increasing contraction ratio. This is caused by 
the increased chamber stay time resulting from the larger L* (6.1 vs 4.8) cm. 
The ERE of the coaxial element is, on the other hand, higher for the lower 
contraction ratio. In this case, it appears that the lower contraction ratio 
enhances the coaxial mixing process by maintaining a higher absolute av for 
mixing so as to overpower the effect of the decreased chamber stay time. 

The chamber axial length (L') has a significant influence on 
the energy release efficiency. The ERE obtained with four of the injector 
elements and the 14 cm (5.5 in.) long chamber is shown in Figure 34a. At this 
chamber length, all four injector elements have ERE's greater than 95%. The 
effect of chamber length is, however, a function of the injector design as 
shown in Figure 34a. The ERE of the premix element, increases only slightly 
from about 96% to 99% as the chamber length increases from 5.08 to 14 cm 
(2 to 5.5 in.). The swirl coaxial element ERE increases from 88 to 98%, the 
coaxial element ERE increases from approximately 67 to 97%, while the F-O-F 
triplet injector increases from 93 to 100% over the same range. The chamber 
characteristic length (L*) has a similar effect on ERE as shown in Figure 34b. 
The L* was varied by contraction ratio or length changes as noted. These data 
also indicate that characteristic length or stay time is the primary variable 
that influences the premix element performance, while physical chamber length 
is the first order influence on the coaxial element performance. 


The effect of chamber pressure and propellant temperature on 
the experimental ERE is presented in Figure 35a for various injector element/ 
chamber configurations and overall mixture ratio. Over the range of chamber 
pressures investigated (45.6 to 341 N/cm2) 70 to 500 psia, little variation in 
ERE with Pc is noted except for the triplet element. Its ERE increased over 
2% when Pq was raised from 68 N/cm^ (100 psia) to 200 N/cm2 (300 psia). 


The propellant temperature shown on Figure 35 is a mass 
weighted average of the fuel and oxidizer manifold temperatures, i.e.. 


T 


P 


1 

0/F + T 


Tf ^ 


0/F 

(ioTTTTT 


T 

0 


(5) 


The fuel and oxidizer temperatures were always within 28°K (50°R) of each other 
and generally were within a few degrees of being the same. In all cases the 
ERE increased with increasing propellant temperature. The premix triplet ele- 
ment exhibited the smallest performance change with propellant temperature, 
increasing 1.0% per 55.5°K (100°R). The coaxial element ERE increased 1.25% 
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per 55.5°K (100°R) while the high velocity swirler element increased 2.7% per 
55.5°K (100°R). Combustion stability variations experienced with the high 
velocity swirler may have influenced the performance variation of this element 
over the tested propellant temperature range. 

2. Heat Transfer 

The detailed calculational techniques are discussed in 
Appendix D. The primary concern of this program was to understand injector 
related heat transfer effects. The data is presented in the form of a heat 
flux ratio, where the local heat flux expr is nominalized by the heat flux 
calculated from a conventional turbulent heat transfer correlation: 

4>nom = 0.026 k/D Rc^-^ PrO.4 (Tq-Tw). Values of one or higher indicate local 
heat fluxes that are higher than nominal and values less than one indicate 
values of relatively low heat flux. 

The design point data for all the injectors is shown in 
Figure 36. This data shows that near zone heat fluxes (in the first 1.27 cm, 
(0.5 in.) of chanfcer length) can vary by an order of magnitude depending on 
injector type. The swirler injector have (|) expr/ci>nom ratios of 0.2 in the near 
zone while the premix injector demonstrated heat flux ratios of over 2. The 
conclusion from this data is that the single biggest variable that influences 
the chamber wall heat flux is the element type. 

The specific results for each element type are discussed in 
the following paragraphs and are illustrated in Figures 37 through 43. 

a. Coaxial Element Axial Heat Transfer 

Characteristics (19 Tests) 

The coaxial element data are shown in Figure 37. These 
data group together reasonably well and a rather gradual increase in heat 
flux with axial distance is indicated. The heat fluxes are generally low and, 
in most cases, are less than one half of the nominal values. Notable effects 
indicated by the data are: 

(1) Contraction Ratio Effect 

Contraction ratio does not appear to exert a large 
influence; however, there is a slight tendency toward higher heat flux ratios 
in the 3.5 contraction ratio chamber (-3 chanter). 

(2) Chamber Length Effects 

A rather sharp increase in heat flux is indicated 
at 4.44 cm (1.75 in.) axial distance in the long 2.5 contraction ratio chamber 
(-2 chamber) but not in the short one (-1 chamber). This behavior is 
unexplained, but it could be an indication that the fiHm cooling effect of 
the fuel annulus deteriorates at about 4.44 cm (1.75-in.). 
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(3) Mixture Ratio Effect 


Relatively high heat flux ratios (up to 0.85) are 
indicated by the 0/F = 1 data obtained at 205 N/cm2 (300 psia) chamber pressure 
in the -3 chamber. This could be largely due to the high mixing rate and thus 
performance of this element at the very low MR operating point. 

(4) Chamber Pressure Effect 

Chamber pressure effects on the heat flux ratio are 
generally very small. The only notable Pc influence exists at 4.44 cm 
(1.75 in.) axial distance in the -2 chamber. 

(5) Propellant Temperature Effect 

The propellant temperature influence appears small 
since the heat fluxes measured at 0/F = 4, Pc = 205 N/cm^ (300 psia) are about 
the same with and without cold propellants. 

b. Swirl Coaxial Element (Vt/Vg = 0.5) Axial 

Heat Transfer Characteristics (7 Tests) 

The low velocity (Vt/Va = 0.5) swirl coaxial element data 
are shown on Figure 38. These data indicate a very low heat flux region exists 
within about 0.6 cm (0.25 in.) from the element face. The heat flux increases 
quite rapidly from this point and is about equal to the nominal value (+25/o 
except for the Pc = 341 N/cm2 (500 psia) test) beyond 3.17 cm (1.25 in.) axial 
distance. 


(1) Contraction Ratio Effects 

Contraction ratio effects appear small. About the 
same results were obtained at 0/F = 5, Pc = 205 N/cm^ (300 psia) with both the 
2.5 and 3.5 contraction ratio chambers. 

(2) Mixture Ratio Effects 

0/F effects are irregular. The 0/F = 2 and 8 data 
obtained at Pr = 205 N/cm2 (300 psia) are about the same, whereas the 0/F = 6, 
Pc = 205 N/cm^ heat flux ratios are about 30% lower. 

(3) Chamber Pressure Effects 

The heat flux ratio is inversely related to Pc 
(decreases Pr)- This is demonstrated by the 0/F = 4, Pc = 68 N/cm2 (100 psia) 
and 341 N/cm2 (500 psia) data as well as the 0/F = 6, Pc = 68 and 205 N/cm2 
data. 


- 74 - 



c. Swirl Coaxial Element (Vt/Vg = 1.0) Axial 

Heat Transfer Characteristics (7 Tests) 

The high velocity swirl (Vt/Vg = 1.0) element heat 
transfer characteristics are similar to those of the Vt/Vg =0.5 swirl ele- 
ment injector as shown on Figure 39. All testing was done at 205 N/cm'- 
(300 psia) chatrter pressure. The only significant difference between these 
data and the Vt/Va = 0.5 swirl element data is the high heat flux ratios 
measured at 0/F = 2 with the long, 2.5 contraction ratio chamber. This 
behavior is similar to that observed with the coaxial element in the long 
chamber. With this injector element, the use of cold propellants reduced 
the heat flux by about 25%. 

d. Premix Element Axial Heat Transfer 

Characteristics (19 Tests) 

Premix element data are shown in Figure 40. A wide 
range of results were obtained with this element as the calculated heat flux 
ratios range from 0.7 to 3.0. The heat transfer environment produced by the 
premix element is fundamentally different from the shear mixing elements and 
is similar to the F-O-F triplet elements, in that the heat flux decreased 
with axial distance except within the 2 cm (0.75-in.) long zone nearest the 
element. This indicates that the premix element produces much more rapid 
combustion than the coaxial or swirl coaxial elements. 

(1) Contraction Ratio Effect 

Two distinct types of heat flux distributions were 
obtained for the 3.5 and 2.5 contraction ratio chambers. The distribution 
drawn through the 205 N/cm2 (300 psia) data in Figure 40a shows that the 3.5 
contraction ratio heat flux ratios are generally higher and that the peak heat 
flux existed further downstream. The 2.5 contraction ratio distribution shows 
a nearly uniform decrease in heat flux with axial distance (Figure 40b and c). 

(2) Chamber Length Effect 

No effect of chamber length is indicated since the 
-1 and -2 chamber data agree reasonably well for the range of 0/F and Pc values 
tested (Figure 40b and c). 

(3) Mixture Ratio Effect 

Mixture ratio effects are very significant with the 
premix element, especially in the 3.5 contraction ratio chamber (Figure 40a). 
The data clearly show a decrease in heat flux ratio with increasing 0/F. This 
characteristic appears related to the momentum of the fuel stream (injected 
perpendicular to the chamber axis in the premix element design) since the fuel 
injection velocity decreases as 0/F increases. This causes a decrease in fuel 
penetration, resulting in more fuel at the outer edge of the jet. 
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(4) Chamber Pressure Effect 


The indicated effect is not great and irregular. 

At 0/F = 4 (3.5 and 2.5 contraction ratio), the heat flux ratio tended to be 
inversely related to Pc; however, the dependency was weak as the variation was 
within + 15%. The opposite effect is indicated for mixture ratios of 2 and 6 
as the 68 N/cm^ (100 psia) data are slightly lower than the 205 N/cm^ (300 psia) 
data (3.5 contraction ratio data). 

(5) Propellant Temperature Effect 

The use of cold propellants at 0/F = 4 and 
Pj, = 205 N/cm2 (300 psia) test conditions reduced heat flux ratios by about 
35%. 

e. F-O-F Triplet Axial Heat Flux Characteristics 

(16 Tests) 

Axial heat flux distributions for the B row (the location 
of this row is shown in Figure 42) of thermocouples are shown in Figure 41 for 
each injector type in terms of the experimental to nominal heat flux ratio. 

The results show that the heat transfer characteristics of the 60° and 90° 
impingement angle injector designs are quite similar, thereby indicating that 
the impingement angle has a small effect on wall heat flux. With both injec- 
tors relatively high heat fluxes were measured at 0/F = 2; where the heat flux 
ratios range from 0.9 to 1.5. This 0/F trend is consistent with the premix 
results discussed previously. The axial heat transfer distribution with the 
60° and 90° impingement angle injectors is considered normal at 0/F = 4 to 6; 
the heat flux increases fairly regularly with axial distance and the heat flux 
ratio is in the 0.2 to 0.8 range. A tendency toward a heat flux peak at 1.0 cm 
(0.75 in.) axial distance is apparent in the 68.6 N/cm^ (100 psia) chamber data. 
A quite different axial characteristic is indicated by the 0/F = 2 data as the 
peak heat flux occurs near the injector face and the heat flux decreases with 
axial distance downstream of this point. The 0/F = 2, Pc = ^05 N/cm2 
(300 psia) data for the 60° impingement angle injector has a characteristic 
almost identical to the previously mentioned premix element; however, the 
external triplet element heat fluxes are about 60% lower than the previous 
data. 


Generally, higher heat fluxes were measured with the 
A^P/Aq = 1 injector element. Heat flux ratios as high as 2.0 were indicated 
as shown in Figure 41. Data were obtained with both the 2.5 and 3.5 con- 
traction ratio chamber designs with this element. The heat flux ratios in 
both chambers are about the same at axial distances of 4.45 cm (1.75 in.) or 
more from the injector face but significantly lower in the 2.5 contraction 
ratio chamber nearer the injector face. 
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Axial and circumferential heat flux distributions are 
given for each injector at 0/F = 4, Pc = 206 N/cm2 (300 psia) operating con- 
ditions in Figure 42. This graph shows that the ^f/^o = 1 injector heat fluxes 
are consistently the highest, but also the most uniform in both the longitudinal 
and circumferential directions. The Af/Ao = 2 injector designs are classical 
"streakers" in that the heat flux along Row C are about 100% higher than the 
heat fluxes along Rows A and B. 

Referring again to Figure 42 it is interesting to note 
that the high heat flux values associated with the Af/A© = 1 injector are con- 
sistent with the high performance of this injector as noted earlier. This cor- 
relation is consistent with a good mixing element, i.e., high heat fluxes and 
a small data spread. The data from the Af/Ao = 2 injectors have a contrary 
characteristic that is again consistent with the performance data, i.e., low 
average heat fluxes with low performance injectors. In addition, these injec- 
tors indicated a chamber hot streak in an area rotated 90" from the plane of 
the fuel holes. These streaks are the result of the poor mixing character- 
istics of the injection elements. The fuel momentum of these injectors is not 
high enough to penetrate the central oxidizer jet. As a result, the propellant 
is incompletely mixed and a near-stoichiometric zone is formed locally at the 
wall. This condition results in the streak noted in Figure 42. It would be 
expected that a full scale injector with Af/A© = 2 elements would behave in a 
similar manner, i.e., lower performance, lower average heat flux with locally 
high values. 

3 . Single Element Stabilit y 

Unstable combustion was encountered during some of the single 
element tests. The test-by-test stability characteristics are listed in 
Table IX. The premix element exhibited stable operation for all test condi- 
tions. The swirler and coaxial elements were, on the other hand, stable only 
at mixture ratios of approximately 6 or greater, while the triplet injector 
exhibited low frequency instabilities at high mixture ratio operating points 
and low fuel side pressure drops. This data is plotted in Figure 43 as a func- 
tion of oxidizer and fuel circuit aP/Pc- At lower mixture ratios, both ele- 
ments were either marginally stable or unstable. The frequency of instability 
varied from 700 to 1100 Hz for the swirler element and 600 to 900 Hz for the 
coaxial element. 


The combustion instability undoubtedly influenced the energy 
release efficiency as noted during Tests 132 (coaxial element) and 162B 
(swirler element). During both these tests, the combustion stability varied 
with attendant performance changes. During Test 132, the combustion would 
shift back and forth from a moderately stable condition to an unstable con- 
dition (800 Hz) and the energy release efficiency would shift from 93.5% during 
stable operation to 89.9% during unstable operation. Similarly, during Test 
162B (a cold propellant test), combustion occurred in two modes: marginally 
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stable and unstable (1000 to 1200 Hz). During the stable mode, the energy 
release efficiency was approximately 83.5% while, during the unstable mode, 
the energy release efficiency was 78.7%. 

The triplet instability at 3500 Hz was evident at mixture 
ratios of four and above for both Af/Ag = 2 injectors. At mixture ratio two, 
these injectors were stable. With the high fuel pressure drop injector 
(Af/Ag = 1h the injector was stable over the entire operating range. The 
triplet data are plotted as a function of (aP/Pc)o vs (AP/Pc)f in Figure 43. 
All injectors were stable at (AP/Pc)f ratios in excess of 10%. No oxidizer 
circuit aP related instability was noted for values of (aP/Pc)q as low as 10%. 
Based on these data, both circuits should be held to aP/Pc ratios higher than 
10% of Pc as a good design criterion. The coaxial and swirler data indicated 
a reverse influence on the fuel side, i.e., high values of aP/Pq were unstable 
while low values of aP/Pq were stable as indicated in Figure 43. 
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FIGURE 29. HOT FIRE SINGLE ELEMENT DESIGN PARAMETERS 




FIGURE 30. SINGLE ELEMENT THRUST CHAMBER 
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FIGURE 32. SINGLE ELEMENT PERFORMANCE 
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FIGURE 34. EFFECT OF CHAMBER VARIABLES ON SINGLE ELEMENT PERFORMANCE 
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FIGURE 36. SINGLE ELEMENT HEAT TRANSFER CHARACTERISTICS FOR VARIOUS ELEMENTS 
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FIGURE 37. COAXIAL ELEMENT HEAT FLUX IN VARIOUS CHAMBERS 
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FIGURE 41. TRIPLET ELEMENT HEAT TRANSFER CHARACTERISTICS 
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VII 


FULL SCALE INJECTOR EVALUATION 


As a part of this prograir roll scc test se''ies was conducted (at 
the 6200 N (1500 1b) thrust leve],’ v.ith irije jrs that were optimized based on 
cold flow and single element testing. The it ent of this testing was to inves- 
tigate performance, heat transfer and stability with full scale prototype hard- 
ware under conditions representative of auxiliary propulsion systems, APS, for 
the Space Shuttle. This testing was accomplished with three basic injectors; 
a deep-cup premix, a F-O-F triplet and two zoned (barrier cooled) F-O-F triplet 
injectors. The results of this testing are briefly summarized below. 

Performance efficiencies approaching 100% ERE were achieved in 7.62 cm 
(3.0 inches) long chambers using both the premix and triplet element injector 
patterns. The zoned injectors were lower performing due to barrier cooling re- 
lated mixture ratio maldistribution effects and due to the core mixture ratio 
shift which caused the core elements to operate off of their uniform mixture ratio 
peak operating point for which they were originally designed. Future barrier cooled 
injectors should be designed so both barrier and core elements are optimized at 
their respective barrier and core mixture ratios to minimize the latter effect. 

High injector-end chamber heat fluxes were recorded for both the triplet and premix 
injectors. Injector face thermal data indicated the need for active cooling of 
the face. All three injectors resulted in at least one case of high frequence 
instability, indicating the need for damping devices such as resonators and/or 
baffles . 


A. DESCRIPTION OF TEST ARTICLES 
1 . Injectors 

a. Deep Cup Premix Injector 

The premix injector designed for this program is con- 
ceptually similar to that injector successfully demonstrated on Contract NAS3- 
14354, Ref. 12. For this program, the basic injector body and manifolding de- 
signs were retained. Only the injection element design was changed. A sketch of 
the element is shown in Figure 44 and it is referred to as the deep cup premix 
element concept. This injector was designed to satisfy the following conditions: 


Number of Elements 
Oxidizer Dia. 

Fuel Slot Geometry 
Pressure Drop 
Injection Velocity 


36 

.60 cm ( .237 in. ) 

.201 X .356 cm (.079 x .140 in.) 
41 N/cm2 (60 psi) 

330 m/sec (1080 fps) 


The oxygen enters the mixing cup through a cylindrically etched passage formed 
of 14 nickel platelets. Fuel is directed radially into the oxygen through a rec- 
tangular passage etched through the same platelets. At the base of the fuel drop- 
through-passage, 6% of the fuel is dive>"ted into the face cooling circuit. A set 
of three platelets redistributes this fuel over the injector face, where it is 
injected axially into the combustion chamber through 384 orifices. 

The deep cup premix concept developed on this program 
utilizes confined mixing of the O 2 /H 2 propellants prior to injection into the 
chamber. This is accomplished by impinging two fuel streams on a central 
oxidizer stream in a mixing cup recessed into the injector face. The basic 
design is similar to the injector that has been demonstrated successfully by 
ALRC on the APS program (Contract NAS3-14354) Ref. 12. This design differs in 
that the mixing cup depth is an order of magnitude deeoer than the APS design 
(0.764 cm deep vs. 0.076 cm deep). 


- 94 - 



This design modification was a result of comparing the 
cold flow and hot firing test performance on this and the APS program. The 
APS single element cold flow data for the premix element with a rectangular 
fuel channel and a shallow cup depth indicated a strong relationship between 
Em and momentum ratio and a relatively low mixing effectiveness. This strong 
functional relationship was not noted in the experimental data from this 
program and the performance level was considerably higher for the deep cup 
design. Also the single element hot fire data obtained on this program 
indicated that the confined mixing imposed by the deep cup increases propellant 
mixing, since the ERE was nearly 100% in a 5.08 cm (2.0 in.) long chamber. 

The injector assembly and platelets are shown in Figure 45. The injector is 
shown at an intermediate step in the fabrication process. The body, platelet 
stack, and individual platelets are depicted. Oxygen enters the injector 
through a 3.8 cm (1.5 in.) diameter stainless steel tube and is distributed 
through the injector body with the aid of sheet metal distribution plates. 

The oxygen enters the face through 36 tubes brazed to the injector body and 
face. Fuel enters the face platelet stack through 36 holes drilled through 
the face plate. The injection platelet stack shown schematically in Figure 44 
is bonded to the face plate, completing the injector assembly. 

The 14 Nickel 200 platelets immediately joined to the 
face plate are the fuel injection platelets, a detail of which is shown in 
Figure 45. The "dumb-bell" shaped passages form the fuel geometry of the 
premix element. Oxygen enters the element through the 36 tubes joined to the 
face plate, and flows into a cylindrical cavity formed by the 14 injection 
platelets. The fuel is injected into the same cavity and intimate mixing of 
the H 2 and O 2 is accomplished. After the propellants are mixed, the gas 
proceeds through the cup plate shown in Figure 45. The function of this plate 
is to provide increased constrained mixing time as well as diverting approxi- 
mately 6% of the fuel into the face cooling circuit. 

The face cooling circuit was comprised of a three- 
platelet stack. The first platelet (the coolant inlet platelet shown in 
Figure 45) was located just below the cup plate, where further H 2 coolant 
metering was accomplished by the 72 0.061 cm (0.024 in.) diameter holes. 

After the coolant passes through this platelet, it was distributed over the 
injector face by the distribution platelet illustrated in Figure 45. The fuel 
leaving the inlet platelet entered the rectangular passages of the distribution 
platelets and flowed over the surface of the injector face, serving to both 
regenerati vely cool the face and feed the face bleed holes. The last platelet 
labeled face in Figure 45 provided the means to inject both the fuel coolant 
and premixed gases into the combustion chamber. 

b. F-O-F Triplet Injector 

Based on cold flow data presented earlier an injector 
composed of an array of external impinging F-O-F triplet elements had the 
potential to be high performing. Also, this element was extremely sensitive 
to design and operating point variables such as relative orifice area and 


- 95 - 



impingement angle. This information was insufficient to design a full-scale 
injector since the effects of combustion were not known but were suspected to 
be important. The experience gained in the single element testing program 
with elements similar to the prototype full scale injection elements, clearly 
indicated that the combustion process did modify the cold flow characteristics 
This data was then used in conjunctions with the cold flow data to complete 
the full scale injector design. An injector was designed based on the design 
criteria discussed in Section VIII, which resulted in the following parameters 


Number of elements 
Oxidizer element diameter 
Fuel element diameter 
Impingement angle 
Oxidizer injection velocity 
Fuel injection velocity 


72 

0.283 cm (0.112 in.) 

0.1325 cm (0.0521 in.) 

7t/ 3 rad. 60° 

86 m/sec (231 ft/sec) 

895 m/sec (2940 ft/sec) 


The triplet injector assembly is shown in Figure 46 
together with details of the injector face. The basic injector body was manu- 
factured from a surplus premix injector from a previous program. Oxygen 
entered the injector via a 3.81 cm p.5 in.) diameter stainless steel tube 
and was distributed throughout the injector body with the aid of sheet metal 
distribution plates. The oxygen was routed to the injection elements through 
72 tubes brazed to the injector body and face plate. Fuel entered the injector 
through a manifold around the oxygen tubes. The hydrogen passed between the 
tubes and entered the elements via holes drilled through the face plate. Up 
to this point the triplet injector and the premix injector were virtually 
identical except that the premix injector was composed to 36 oxidizer tubes 
rather than 72 for the triplet. 

Figure 46 shows the injector element. It was machined 
into a 0.47 cm (0.185 in.) thick copper plate that contains additional fuel 
injection orifices, the oxidizer injection orifices and the face coolant 
metering orifices. This plate was bonded to the injector body such that the 
oxidizer tubes were aligned with the oxidizer injection orifices and the fuel 
galleries were aligned with the holes through the body plate. This plate 
serves the dual function of directing both propellants into the chamber. 
Immediately below the face plate was a set of three platelets that comprised 
the face cooling circuit when bonded to the injection plate. These platelets 
were made of 1/4 hard nickel 200 and were 0.0254 cm (0.010 in.) thick. The 
first platelet, the coolant metering plate, controlled the amount of GH^ 
directed into the coolant circuit (7% of the total fuel flow). Fuel exited 
this plate and was directed into the distribution platelet. Fuel entered 
each shaped passage of this platelet and was distributed throughout the 
injector face. Covering these passages was the last platelet in the stack, 
which injected the coolant into the combustion chamber through 460 0.038 cm 
(0.015 in.) diameter holes. 
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Zoned Injectors 


The third full scale injector was a zoned injector, 
designed to produce an environment of lower temperature barrier gases composed 
of fuel rich combustion products. This '-as achieved by tailoring the triplet 
injector described earlier to produce a low mixture ratio barrier at the 
chamber wall by modifying the element design. The zoned injector was manu- 
factured from the triplet injector by first removing the triplet face and then 
bonding a zoned face and cooling circuit onto the triplet body. Two different 
injector patterns were tested, the -1 and -2 injectors. These injectors 
differed in the values of the barrier mixture ratio and mass flow. The -2 
injector was made from the -1 injector element by increasing the diameter of 
the oxidizer injectioti elements in the outer row. A detail drawing of the 
injection elements and patteni is illustrated in Figure 47. 

The flow distribution in the core and barrier circuits 
of the -1 and -2 injectors were determined empirically, by cold flowing the 
injector with GNo. The individual orifice flow rates were measured with a 
rotameter type flow meter, with the injector flowing GN 2 at the design Mach 
No. The barrier mixture ratios and flow rates were then computed and are 
plotted in Figure 48a. At an overall mixture ratio of 4.0, the -1 injector 
flowed 35. 6X of the total injector flow in the barrier at an 0/F of 2.93 while 
the remaining flow is in the core at an 0/F ratio of 4.9. The -2 injector 
(with increased diameter oxidizer orifices) at an overall 0/F ratio of 4.0 
flowed 39.2% of the toal injector flow through the barrier circuit at a 
mixture ratio of 3.3 while the remaining flow was directed through the core 

at an 0/F of 4.6. 

The zoned injector cooling concept results in a per- 
formance loss due to off design mixture ratio operation of the barrier and 
core circuits. The specific impulse loss was computed by assuming that the 
combustion zone was comprised of two streamtubes operating at the respective 
core and barrier mixture ratios and mass fluxes. The streamtubes were 
assumed to expand one dimensionally through the full engine area ratio without 
mixing. This analysis was completed by calculating an overall specific impulse 
consistent with the sum of the mass weighed specific impulse of the individual 
streamtubes. The results of this analysis are illustrated in Figure 43b. The 
specific impulse loss for these injectors at an overall mixture ratio 0 
is 0.8% and 0.5% of one-dirnensional kinetic specific impulse (Icp qdk) 
the -i and -2 injectors, respectively. Additional performance loss was 
expected because the elements were operating away from the optimum injection 
momentum ratio as discussed in Section VIII. 

2. Chambers 

The injectors were tested in the heat sink copper chambers 
illustrated in Figure 49. These chambers had provision for high frequency 
instrumentation, two static pressure ports and inner wall, C-A thermocouples 
located both between and in line with the i njecti on elements . In addition a 
port was provided for a 2 grain RDX charge for stability evaluation. 
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Two chambers (-1 and -2) and a L* section were fabricated. 

The -2 chamber differed from the -1 chamber in that the throat diameter was 
larger resulting in a smaller contraction ratio. Provision for increasing 
chamber length from 7.62 cm (3.0 in.) to 13.9 cm (5.47 in.) was provided by 
a cylindrical copper section also instrumented with flush-mounted thermocouples 
and both high and low frequency response pressure transducers. The salient 
chamber design parameters are listed in the following table. 


-1 Chamber 


-2 Chamber L* Section 


Length - cm (in.) 7.62 
Throat Diameter cm (in.) 4.87 
Chamber Diameter cm (in.) 8.94 
Contraction Ratio 3.4 
L* cm (in. ) 13.5 
Area Ratio 2.97 


(3.0) 

7.62 

(3.0) 

6.27 (2.47) 

(1.92) 

6.27 

(2.47) 


(3.52) 

8.94 

(3.52) 

8.94 (3.52) 


2.02 



(5.32) 

8.37 

(3.29) 

39 (15.3)0) 


2.97 




B. TEST FACILITIES AND PROCEDURES 

Physics Laboratory Test Bay 6 was selected for the full-scale tests. 
The Bay 6 flow and instrumentation schematic together with the test hardware 
installed in the bay are shown in Figure 50. This bay is equipped with complete 
gaseous hydrogen/gaseous oxygen feed systems with sufficient capacity to run 6600 
N (1500-lbf) engines for 10-sec duration. The facilities include a spherical 
heat sink hydrogen conditioning vessel that was used for the temperature condi- 
tioned propellant tests. Oxygen was temperature conditioned by flowing LN^ 
through a small packed bed. The Op then was chilled by transfer of heat within 
the bed. 

The mass flow rates were measured with critical flow venturis. 

Both ambient and propellant temperature conditioned tests used gas regulators 
for setting the supply pressure to critical flow nozzles. These nozzles 
provide reliable and accurate flow rate data while maintaining a constant 
flow rate to the engine independent of downstream perturbations. System 
pressures and thrust were measured using a six-wire measurement system 
utilizing strain gage pressure or force transducers. The transducers have 
350 Ohm strain gauges in a fully active four-arm bridge configuration. 


C. FULL SCALE TEST RESULTS 


Fifty four full scale injectors tests were completed. The test 
series for the premix, triplet and zoned triplet injectors involved 21, 19 and 
14 tests, respectively. A detailed list of test parameters, test configuration 
and test results is found in Tables X and XI. Table X is in English Units and 
Table XI S.l. units. The following paragraphs summarize the performance, heat 
transfer and stability results from this series of tests. The techniques used 
to yield these results are discussed in Appendices C and D. 


rn With -1 Chamber 
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DEEP CUP PREMIX INJECTOR 
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TEST DATA FOR 
ZONED INJECTOR 
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TEST DATA FOR 
DEEP CUP PREMIX INJECTOR 
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TEST DATA FOR 
I ZONED INJECTOR 


00 O '3' O' 


o csi m 00 


O' -4- 


O' ^ O' O' O' ^ 


"^<M POfOfOtO'COO 
O' O' O' O' O' O' O' O' 


CJ 

M-t 

o 

(M 



vO O r»- 
00 O' 00 
m fo on 


m 

O' 

00 

tn 


00 

(4 

Oi 


O 

0> 

W 


<N <n CM lO 

to O T-I \0 00 

to 'd' tO vO 1 /^ 

CO CO cn CO CO 


CM 

O' 

vO 

CO 


CO'OtOCOtO— lOl^ 
tOvorsirOvOcM'O'O 
'O »-i '£> vO CO r — ■ to \0 
forocofororofoco 


* 

u 


o 

V 


to \D to to 

CM CM »-< sO to 

CM O CM CO CM CO 

CO CO CO CO CO CO 


^ O' 

lO CO to lO 
CM OO CM CM 
CO CM CO CO 


00 to CO CO 
CM <M CM <£> 
O CO CO 
CO CO CO CO 


0 

01 
in 




sO O OO r~~ CO 

to o to to to 'C 


^ ^ ^ ^ CM 


-d-iOCM^t-ICMr-IMO 

iO lO -vO \0 vO \0 vO nO 




tA in to lo ^ to 


<N rsi m oj m Oi 


< 

o 




\0 '■i' O' O CM tO 

o CM o — I 00 


■O’ vO CM CO CO 


tOtOiOiAiOtOtOtO 
O' O' O' O' O' O' O' O' 


CMCMCMCMCMCMCMCM 


•OiOi'.O'Dr'CXJr- 

'O'OO'OCMOO'tO 


cO C- CO MT 'i5 CM CO 



CO 

O 

O 




CO 

O' 

o 

o 




CO 

o 

o 

o 




CM lO CM O' 
^ fo ^ ^ O' r-« 


O 'O O CM CO CO 




tO 03 to to CM 
CM CO CM O '3' r-l 


OO-OcOiOCOCOCOO 

r'OOOOcOOO'D 


CO r-- -ct vO CM co 


— IQOCMCMCMCO.— (00 
^^COCOCMiOOOiO 


cor^Ocj'O'OO'oo 

^cococMco—tr^O 



CJ' O' 00 o 


sD 




H O 


OO 00 o. lO to to 
O' O' ^ ^ O' 
nj CM rsj 04 rsi CM 


tO'OtOCM'OvO'Dr^ 

O'O'OOO'tOCJ'O'iO' 

CMCMCMCMCMCMCMCM 


^ CO 00 00 O O 

^ oO CM CM O O 
O' 00 O' O' O' O' 
CM CM CM CM CM CM 


t£> O 




tr> v£> O 00 f" CO 
vD C'- \0 f'- r-- 
CM CM CM CM CM CM CM 


CM -Ct 


Mt CO O' 


OO O' \0 C 


CM CT' . O' •d’ CO cx: 

o o o o o o 

CO CM CM CO CO CM 


r'OOLOO'cM'O*— <■— ( 

ocoeoooocj'oo 

CMCMCMCMCMCMfOCO 


> 

00 

> 

CIS 

b 


Z 


r- 00 to CM CM 


•4- to UO ^ 

■O CO CO O CO o 

to Mf sO C" O 00 

iTl lO to m CM O' 


^iOCT'OOlO'OiO 


r-Cr'.Vj'vOO'O'CM^ 

OtOoor-vO'J'OO 

00'00c»<t0r-r'« 

fM<J’UOuOiA'OiOO' 


B 

o 


CM 00 CM 

vO . vO 


uO lO ^ O <— I O' 
CO ^ to ^ vS- ^ 

O CJt O r-l O- 'O 
CM ^ CM CM CO 


cj'MtrotncM-J’O'io 


cM^tOtOOCM-d-r^ 
r'»Q0r-C,-(CM<MOvO 
^ CM CM CM CM CO 


M3 CM to CO O' CO 


CO to tO —( -O CO 

O O' O r-C '£> 

CM ^ CM CM CO 


iO^O'iOC 300‘00 


cocMr^p'^ Ocor^.O' 
OhOO^^OCMOsO 
r-l CM CM CM CM CM CO 


e -■c r- to 

cj m CO 

•JC CO tO OO 

J —I CO 


jG 

o 


DO OO CO 

CO CO O 
CO CO CM 


CO « 
(A 
3 
Q 


o ^ ^ xO o o 

o. o CO cO 00 


o -C -t O ^ O 


OS 

4J 

ns 

a 


00 ob 00 X OO 00 


o o o o o o 


CM CM CM CM 


CM CM CM CM 


o 



I 1 


cOco-^^-ct-i-^^^i- 

CMCMCMCMCMCMCMCM 

t I I I I I I I 

oooooooo 


r- X O' O -J CM 
sj >J m to m 



r'-XO'Oi-lf'JfO'sS’ 
in uS in nD so vo vo o 



CM 


- 102 - 



Performance 


1 . 


A summary of the full scale performance data is presented in 
Tables XII and XIII in English and S. I. units, respectively, and is plotted 
in Figures 51 through 53. The performance is listed in the following table 
for the four injectors at the nominal design point (206.7 N/cm^ (300 psia) 
chamber pressure, 0/F = 4, 7.62 cm (3 in.) chamber length). 

Injector Delivered Performance Comparison 


Injector 


Energy Release 
Efficiency - % 


Triplet 100 
Premix 98.4 
-2 Zoned 96.8 
-1 Zoned 96.2 


The energy release efficiencies for both the triplet and premix injectors are 
very high in a very short chamber. The achievement of this performance level 
underscores the value of conducting single element cold flow and combustion 
testing. 


Triplet Injector Performance 

The triplet performance was the highest recorded. In 
chambers of 7.6 cm (3.0 in.) length the performance of this injector was 100% 
ERE, as shown in Figures 51 and 52. The data groups very closely and with 
the exception of mixture ratio, no single parameter appears to affect the 
performance level of the injector. The highest performance occurs at the 
mixture ratio four, design point. This is consistent with the design 
parameter MP generated from the single element testing and discussed in 
Section VI. The full-scale triplet data is replotted in terms of this design 
parameter in Figure 53, which is a curve fit of the test data. This figure 
indicates an optimum performance in the range of MP = 0.75 a value consistent 
with what was predicted based on cold flow and single element testing. One 
deleterious aspect of this design point is that it leads to very high pressure 
drop fuel elements or very low pressure drop oxidizer elements. The observed 
drop in performance on both sides of the optimum can be explained on the basis 
of the cold flow work accomplished earlier in the program where it was apparent 
that at low values of MP the fuel momentum completely penetrated and split the 
oxidizer jet. At high MP values, the low fuel momentum jet failed to penetrate 
the oxidizer but merely surrounded the oxidizer jet which resulted in unmixed 
flow. At the design point, optimum penetration is achieved and a well mixed 
homogenous fan is formed. 

The triplet ERE as a function of chamber geometry is found 
in Figure 52. Performance is not a strong function of chamber length 
especially at 0/F = 4. This fact is not surprising considering the high 
overall performance in the short chamber. In fact, it is probable that good 
performance could be obtained in a 5.08 cm (2 in.) long chamber. Contraction 
ratio influences were insignificant. 
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TABLE XII 

aerojet liquid rocket company 
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TABLE XII (cont.) 
triplet Injector *♦**»♦•»*» 
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ENGLISH UNITS (SEE TABLE NOMENCLATURE LIST) 
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TABLE XIII (cont.) 
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f Irij ector Pe rf ormance 

The premix performance data is plotted as a function of 
mixture ratio in Figures 51 and 52. These figures indicate that ERE is 
relatively insensitive to 0/F changes. ERE was over 98% across a mixture 
ratio range extending from 1 to 8. This characteristic was observed in both 
the single element combustion and cold flow testing. Insensitivity to 
mixture ratio is an advantage from an engine system standpoint. Shown in 
Figure 51 are the influences of P and propellant temperature. Lowering P 
is in the direction of slightly i^icreased performance [approximately 0.1% 

ERE per 68 N/cm^ (100 psia)] and lowering temperature is in the direction of 
lower performance [approximately 1% ERE/55°K (100°R)]. The effect of chamber 
geometry is illustrated in Figure 52. As shown in this figure varying contrac- 
tion ratio did not influence the performance. Increasing the chamber length 
from 7.6 to 13.9 cm (3.0 to 5.5 in.) resulted in changes in performance that 
were small enough to be within the data scatter (+1?0- 

Zoned Injector Performance 

The ERE of the zoned triplet injectors is presented in 
Figure 51. These injectors deliver substantially less performance than either 
the triplet or premix injectors. The low performance of the zoned triplet 
injector when compared to the triplet injector is due to two effects: (1) the 
mixture ratio distribution loss due to zone cooling and (2) the off-optimum 
element design dictated by the barrier concept. As discussed earlier in this 
report, 0.8 and 0.5 percent I loss can be attributed to mixture ratio 
maldistribution for the -1 ana"-2 zoned injectors, respectively. If one 
accounts for this loss, the remaining difference between the zoned triplet 
and triplet injectors must be attributed to injector element design. 

Referring to Figure 53, where the triplet element ERE as a 
function of the mixing parameter is displayed, note the design point mixing 
parameter values for the core and barrier circuits of both zoned injectors. 

Both circiuts were composed of triplet elements removed from the optimum of 
a F-O-F triplet element design. This departure from an optimum design is 
forced because the barrier flow split is controlled by injection element area. 
Without a complicated injector manifold (i.e., separate barrier and core 
circuits) an overall optimum element can not be obtained. Therefore, in any 
engine that uses F-O-F triplet elements, performance losses will result from 
both the effects of barrier coolant (mixture ratio distribution) and non- 
optimum element design. Based on this concept and the work presented in 
Section III on film vs varrier cooling, it would appear that film rather than 
barrier cooling would be a better design approach. In general, any GO 2 /GH 2 
injector-chamber that is comprised of injection elements that are mixture 
ratio sensitive, should be cooled with a film rather than a barrier. 
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2. Hedt Transfer 


a. Chamber Heat Flux 

The full scale charber heat transfer was evaluated in 
a manner similar to the single element data. This approach is discussed in 
Appendix D and it relies on a ratio of measured experimental heat flux 
nominalized by the nominal heat flux based on conventional turbulent heat 
transfer correlations (Bartz). This nominal ization accounts for differences 
in Pc, 0/F, and chamber diameter. The direction of this ratio away from one 
therefore is a measure of the injector related influences on chamber heat flux. 

The data for the various injectors is plotted in 
Figures 54 through 63 and is discussed in the following paragraphs. All four 
injectors are summarized in Figure 54 for the nominal test conditions (0/F = 4, 
Pc = 205 N/cm2, (300 psia) -1 chamber). 

The experimental chamber heat fluxes for the triplet 
and premix injectors when compared to the predicted heat flux (as defined by 
the simplified Bartz relation) were high in the injector end of the chamber, 
decaying to one at the throat. This trend and the values of the heat flux 
ratio is nearly identical to the data from the single elements. The chamber 
heat flux using the zoned triplet injectors was a factor of three less than 
the triplet injector near the face. This decrease in heat flux is not entirely 
due to the effects of zoning. It is in part generated by the inefficient 
injection elements of the barrier circuit. These elements, because they are 
not well mixed, do not release heat equivalent to an optimum triplet element. 
Therefore, the measured flux levels reflect a depression generated by both 
the barrier cooling and decreased energy release • 

(1) Premix Injector Data 

Heat flux ratios calculated from the premix injector 
data are plotted in Figures 55, 56 and 57. Figure 55 shows the -1 chamber 
data for Pr = 68 N/cm2 (100 psia) (0/F = 2, 5, 7) and Pc = 206 N/cm2 (300 psia) 
((/F = 1, 2, 4, 8). Figure 56 shows the -1 chamber data for = 344 f'i/cm<^ 

(500 psia) (0/F ^ 2, 4, 6) and the -2 chamber data (Pc = N/cm^ (175 psia) 
and 0/F = 2, 4, 6). Figure 57 shows the premix data obtained with the L* 
section positioned between the injector and the -1 chamber. 

The premix data indicate that mixture ratio does 
not influence significantly the thrust chamber heat flux. Also chamber 
pressure does not appear to influence the heat flux ratio distributions. This 
is demonstrated on the bottom graph of Figure 55 where the band drawn for the 
68 (100 psia) and 344 N/cm2 (500 psia) chamber pressure data matches the 
206 N/cm^ (300 psia) data. There is a slight tendency for the -2 chamber 
(CR = 2.0) heat flux ratios to be less than the -1 chamber (CR = 3.4) values 
(bottom graph of Figure 56). 
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A distinct pattern of heat transfer characteristics 
is evident in the premix injector/short chamber data (Figures 55 and 56): 

(1) Heat fluxes are somewhat above nominal near the injector face (heat flux 
ratio = 1.1 to 2.0); (2) The heat flux ratio increases quite sharply downstream 
of the injector and a peak ratio of 2.0 to 2.9 occurs at 1.9 cm (0.75 in.) 
axial distance; (3) The heat flux ratio decreases from this point until it is 
1.0 or slightly less (heat flux about nominal) at the throat; (4) The heat 
flux remains at a nominal level in the divergent nozzle section to an area 
ratio of 2.3. 

The data from tests using the -1 chamber and the 
L* section is shown in Figure 57. These data demonstrate that the high heat 
flux ratios observed in the cylindrical section of the shorter chambers are 
a function of distance from the injector face. With the L* section in place, 
heat fluxes in the -1 chamber cylindrical are near nominal and the high heat 
flux ratios exist upstream in the L* section. Apparently the -1+L* cylin- 
drical length is sufficiently long that a near fully developed boundary layer 
exists at the downstream end, since the throat data are slightly lower than 
the short chamber data. 

(2) Triplet Injector Data 

Chamber heat flux ratios obtained from triplet 
injector tests are plotted in Figures 58, 59, and 60 in the same format used 
for the premix data. Figure 58: -1 chamber results, = 68 N/cm (100 psia) 

0/F = 2, 4, 5) and Pc = 206 N/cm2 (300 psia) (0/F =1,2,4, 6). Figure 59: 

-1 chamber results for Pq = 344 N/cm^ (500 psia) (0/F = 2, 4, 6), and -2 
chamber results for Pc - 119 N/cm^ (175 psia) (0/F = 2, 4, 6). Triplet data 
for the -1+L* chamber configuration are shown in Figure 60. 

Over most of the mixture ratio range tested, the 
0/F effects are similar to those noted for the premix injector; slight, 
irregular changes in heat flux ratio are produced by changes in 0/F. 

However, a fairly regular low 0/F effect is evident in the triplet data as 
the heat fluxes tended to be lower than nominal at the lower 0/F values. 

This occurred at 0/F = 1, Pc = 206 N/cm2 (300 Psia) (bottom graph. Figure 58) 
and at 0/F = 2, Pc = 344 N/cm^ (500 psia) (top graph. Figure 59). This 
reduction in heat flux may be a film cooling like effect produced by the 
outer impinging fuel elements when excess amounts of fuel and high fuel 
injection momentum create a fuel rich zone adjacent to the wall. 

Some Pr effects are evident in the triplet data. 

The heat flux ratio 0.634 cm (0.25 in.) from the injector face tends to 
decrease with increasing Pc (Figure 59). Heat fluxes in the -1 chamber were 
slightly higher than the -2 chamber values at comparable MR and Pc conditions. 

Heat flux ratios calculated from the triplet 
injector/-l+L* test data, shown in Figure 60, are similar to the premix 
injector. The somewhat lower heat flux ratio calculated from the -1+L data 
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nearest the injector face may be due to the use of the ID data reduction 
model for the L* section (see Appendix D). 

(3) Zone-1 Injector Data 

Data for the Zone-1 triplet injector are plotted 
on Figure 61. The top graph of Figure 61 shows that the effect of converting 
the triplet injector to the Zone-1 triplet design was to produce a pronounced 
decrease in heat flux ratio in the cylindrical chamber region (compare with 
top graph. Figure 58). The nominal heat flux decrease due to the decrease in 
O/F along the wall is only 10 - 25%, but the data show that the actual 
decrease ranged from 15% to 60% at 1.9 cm (0.75 in.) axial distance. The 
largest effect was at 0/F = 4 and the smallest effect was at 0/F = 2. Heat 
fluxes downstream of the cylindrical section were about the same as the 
triplet. Chamber pressure effects are small except that streaking at the 
throat occurred at 344 N/cm2 (500 psia) chamber pressure (bottom graph, 

Figure 61). This streaking was probably due to the non-optimum element 
design of this injector. 

(4) Zone-2 Injector Data 

The Zone-2 triplet injector heat flux ratios are 
plotted on Figure 62. Comparison of the P(~ = 206 N/cm^ (300 psia) data (top 
graph) to the triplet results (top graph. Figure 58) reveals that the Zone-2 
triplet injector also produced lower heat fluxes in the cylindrical region. 
However, very pronounced streaking occurred downstream of the cylindrical 
length, and the Row B heat fluxes were higher than the triplet injector heat 
fluxes in the convergent region; again the result of non-optimum element 
design. The low heat fluxes observed at 0/F = 1 with the triplet injector did 
not occur with the Zone-2 injector. No Pc effect is apparent in the Zone-2 
data since the Pc = 68 N/cm2 (100 psia) and 344 N/cm2 (500 psia) data compare 
reasonable well with the 206 N/cm^ (300 psia) data band (bottom graph, 

Figure 62). 


(5) Generalized Correlation of the Data 

The chamber throat heat transfer data are compared 
to traditional turbulent and laminar boundary layer correlations in the graph 
shown as Figure 63. The graph shows the St PrO-6 parameter as a function of 
throat Reynolds number. This is the approach utilized in the Reference 13 
study of rel ami nari zati on . The turbulent equation shown on Figure 63 is the 
equation used to calculate nominal heat flux values. Data from the ALRC APS 
thruster program (Reference 12) are also shown in Figure 63. For the purposes 
of Figure 63 only, the MR 4, 6, and 8 heat flux ratios were divided by 
correction factors which adjust the data to a nominal heat flux based on 
enthalpy driving potential rather than temperature driving potential. This 
was done so as to be consistent with the Reference 13 approach. The correction 
factors used were: 1.1 for 0/F =4, 1.4 for 0/F =6, 1.5 for 0/F = 1.8. 
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The data indicate a tendency toward laminar 
boundary layer behavior at Reynolds numbers below 5 x 105 (about the same 
transition point noted for the APS data). At higher Re values, the data 
generally agree with the turbulent equation although a slightly lower 
coefficient (0.023 instead of 0.026) is indicated. The trend toward 
laminarization is more gradual than indicated by the APS results. 

A few of the data appear unusually low as they 
compare well with the laminar correlation at Reynolds numbers in excess of 
106. These data are the Zone-2 row B throat values for Pc = 206 N/cm^ (300), 
0/F = 1 and 8; the Zone-2 row B value for Pg = 68 N/cm2 100, 0/F = 4, and 
the triplet value for Pg = 206 N/cm^ (300, 0/F =1). It is not clear why 
these four data points are so low, combustion streaking effects may be a 
factor in these data. 

b. Injector Face Heat Transfer 

The face temperatures of all the injectors are listed 
in Table XIV. Tests 101 - 109, 134 - 144, 154 and 155 are with the premix 
injector. The thermocouple locations can be referenced to the injector face 
by referring to Figure 45 of this report. This data indicates a dependence 
of face temperature on chamber pressure i.e., low chamber pressures result 
in lower face temperatures. It is also apparent that the face temperatures 
are hotter than would be desired with a long cycle life injector. They are 
as high as 844°K (1060°F). 

The triplet data are shown in Tests 111 - 133. The 
locations of these thermocouples can be determined by referring to Figure 46 
of this report. The data presented in Table XIV is not entirely representative 
of the actual face temperature as the locations of TJ-1 and TJ-2 were in a very 
cool portion of the face. In fact, the visual appearance of the injector face 
indicated temperatures in excess of 925°K (1200°F), a value too high for long 
cycle life injectors. 


The zoned injector data is included in Tests 147 - 152B 
and 157 - 164. This injector was designed and manufactured with increased 
cooling capability when compared to the triplet injector. The injector was 
constructed with a face bleed circuit designed to flow 13% of engine fuel 
flow. The triplet, by comparison had a design bleed flow of 6% of the total 
fuel. Inspection of this data indicated temperatures as hot as 1084°K, 

(1491 °F) at the 0/F = 2 run condition. Mixture ratios on either side of 0/F = 
2 were in the direction of lower temperature. Again face temperatures are too 
high for long cycle life applications. By visual comparison, however, the 
effects of increased coolant did reduce the overall face temperature, although 
local hot spots were still present. 
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INJECTOR FACE TEMPERATURE 
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INJECTOR FACE TEMPERATURE 
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3. Stability 


The stability data for this program are listed in Table XV 
by run number. The relevant run conditions and the stability conditions 
including frequencies and amplitudes of unstable tests are included. As it 
was beyond the scope of this program to attempt to formulate a generalized 
stability model for gas-gas combustion, the data presented are intended to 
be used for further development of existing analytical models.* To permit full 
utilization of this data, a complete description of the engine system is given 
in Figure 64 with the relevant engine flow parameters. 

A summary of the data reveals the following stability 

characteristics : 

(1) The triplet injector exhibited unstable combustion in 
the first radial and third tangential mode at mixture ratios near four. At 
mixture ratios two and six the injector was stable. 

(2) Using the premix injector, two tests with unstable 
combustion were observed at mixture ratio of four. AH the remaining tests 
were stable. This instability appeared to couple with the first tangential 
mode of the chamber. 

(3) The -1 zoned injector produced unstable combustion at 
a mixture ratio of six which coupled with the first tangential mode of the 
chamber. The remaining -1 and -2 tests were stable. 

Since at least one case of high frequency instability was 
noted with each injector, stability damping devices probably would be needed 
on a production engine. A designer could incorporate either a resonator or 
face mounted baffles in any injector design. 


Recently emerging gas/gas combustion stability models are underwork at NASA 
Lewis Research Center. Dr. R. J. Priem using this data was able to successfully 
correlate at least a part of this data. This model will be published by 
LeRC as a NASA TN in the near future. This document title is listed in 
Reference 20. 
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GAS/GAS INJECTOR STABILITY DATA 



Bomb 2 grain RDX 



TABLE XV (cont.) 

GAS/GAS INJECTOR STABILITY DATA 



157 3.38 5.32 13.51 107 73.8 0.948 0.430 0, 251 0.L14 







cm 7>iRU 









Mixture Ratio Distribution 

Specific Impulse Loss Barrier & Core Mixture Ratio 

7o of ODK 


ZONED INJECTOR 
FLOW CHARACTERISTICS 




-Barrier 


Lines -2 Injector 
Lines -1 Injector 



Overall Mixture Ratio 


PREDICTED SPECIFIC IMPULSE LOSS 
DUE TO ZONED COOLING 


-1 Zoned 
Injector 


-2 Zoned 
Injector 


- 124 - 


Barrier Flow 





Figure 50. Full Scale Test Hardware and Assembly 
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FIGURE 51 . EFFECTS OF MIXTURE RATIO AND CHAMBER PRESSURE ON FULL SCALE PERFORMANCE 
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NG PARAMETER ON TRIPLET ELEMENT PERFORMANCE 
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Mixture Ratio Effects 
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FIGURE 62. ZONE-2 TRIPLET CHAMBER HEAT FLUX 
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Vni. DATA EVALUATION 


A. INTRODUCTION - METHODS OF APPROACH 

The objective of this task was to synthesize analytical models that 
were based upon the experimental data of this program. The resulting design 
equations and curves were prepared to provide an injector designer with a tech- 
nique to estimate performance and chamber heat flux for prototype injectors. 

This work culminated in the procedures contained in the Design Handbook. The 
intent of the modeling work was to bring together all aspects of the performance 
and heat transfer problem and unify the data within the structure of a firmly 
based analytical formation. It became clear early in the program that these 
aims could best be served by two fundamentally different approaches to the 
mixing/combustion problem: (1) one based entirely on empiricism tempered by 

physical fundamentals and (2) one based on an analytical approach tuned with 
empirical data. 

The first approach has as its advantage ease of computation and 
accuracy. Its disadvantage is that it is tied exclusively to specific injec- 
tion element and test conditions and can not be extrapolated beyond the tested 
envelope with confidence. The second approach is more complex, but it has 
modeled the first order mixing/combustion mechanisms. Thus it is general in 
nature and can handle most design problems. However, approximation had to be 
incorporated to make the model practical which reduces the accuracy of the 
model as a prediction tool. 

The empirical approach is based on correlation of the cold flow 
data, leading to design equations that predict mixing efficiency (Em). These 
equations are modified (where data was available) by appropriate constants to 
account for combustion influences. The basic combustion influence was treated 
by an "equivalent-length" function. This equivalent cold flow length is the 
technique by which the cold flow based equations can be scaled to equivalent 
combustion conditions. Energy release efficiency is obtained by scaling Em 
by an empirically determined function that relates the two parameters. 

In addition to the performance aspects of the model, chamber wall heat flux 
data was correlated. This correlation estimates the local flux in the forward 
end of the chamber as a function of injection momentum ratio. 

The analytical approach results in much the same information except 
that it is not element-specific. The details of the analytical modeling are 
found in Appendix A. In this section a generalized approach is discussed that 
simplifies the rather tedious computational procedures implicit in the analyti- 
cal model for clarity in presentation of the concept. The discussion in this 
section is a result of correlations formed by using the analytical model. 

These correlations result in four non-dimensional parameters that characterize 
the mixing process for both axisymmetric and two dimensional jets. By 
appropriately selecting and/or calculating the various parameters contained 
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within this correlation virtudUy any element or test condition can be analyzed. 
Combustion is treated via an analytical technique that considers the thermo- 
chemical impact of combustion in conjunction with the mixing process. In addi- 
tion, chanter heat flux can be predicted based on computation of outer-edge 
combustion gas properties. 

The following paragraphs describe the details of both approaches. 

The empirical approach is discussed first followed by a synopsis of the 
analytical approach. 

B. DISCUSSION OF RESULTS - EMPIRICAL APPROACH 

The objective of this section is to document the work that resulted 
in the design equations and plots presented in the Design Handbook. The approach 
taken here was empirically based as opposed to analytical approacn described in 
the next section. Basic design equations were formulated by correlation of the 
cold flow data with the relevant data from the single element hot fire senes. 

The effects of combustion were treated by empirically determining a functional 
relationship between cold flow length and contiustion length. This functional 
relationship was unique for each element type, and allows the design equations 
based on cold flow data to be scaled to a combustion environment. In order to 
predict combustion performance, it was also necessary to determine a relation- 
ship between mixing efficiency and corribustion efficiency. In addition, the 
single element and full scale injector heat flux data were correlated as a 
function of momentum ratio; resulting in a technique that estimates the char- 
acteristic heat flux for various elements. The following paragraphs describe 
the details of the cold flow correlation; the functional relationship between 
E[j| and combustion efficiency; the determination of an equivalent cold flow 
length, and the correlation of the heat flux data. 

1 . Cold Flow Data and Correlations 

The cold flow data discussed earlier in Section IV, were used 
to form correlations between E^ and various operating and design variables. 

These correlations were intended to organize the data into an easily used for- 
mat that accounted for the effect of the primary variables on Em. In general, 
this aim was realized. Single element data was correlated with an equation 
that is of the following general form. 

= K f a/D^, 0/F, D„/Df, ApA^) (6) 

where K is a function of the particular variables plotted vs K in Figures 65 
through 69 for each element type. 

The technique used to correlate the data was similar for each 
element type. Each variable of the cold flow data was investigated to ascer- 
tain its particular influences on E^. The variables were then grouped and 
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plotted point by point versus experimental Em. If second order influences then 
become evident, a new variable prouping was synthetized that accounted for the 
observed influence and the resulting data were replotted again point by point 
versus experimental Em. This process was continued until the experimental data 
grouped within + or better. The experimental Em versus the variable-grouping 
then defined a unique relationship between Em and the variable-group. 
relationship was then put into equation form by defining a new variable (K) that 
was a function of the relevant physical parameters and resulted in Equation 6. 
The single element and full-scale combustion performance data was then investi- 
gated to see if the equations would predict the absolute values and trends of 
the combustion performance data. In general, the combustion did influence the 
mixing process resulting in modifications of the cold flow equations. These 
modifications were constants that were applied to the cold flow equations to 
account for combustion influences. The resulting equations, therefore, were 
formulated such that one constant was used for the cold flow correlation and 
another constant was used for the combustion correlation. The following 
paragraphs describe the results of this work for each element type. 


Coaxial Element 

The coaxial element and increased shear element data is cor- 
related in Figure 65. Referring to Figure 65, note that the equation: 


E = 100 - 26.2 In (rTfi^) 


(7) 


where K is a variable plotted in Figure 65, correlates the experimental 
data to^within + 5% over the entire ranges of design and operating variables 
tested. Inspection of the above equation and the relationship with 
V -V 

— 5L) reveals the following information: 


Em is related inversely to L/Dp. (Long chambers and small 
diameter elements are in the direction of high performance.) 

V.-V V.-V^ 

Em is related to -l^with high -J— ^ in the direction of 

^0 0 

high performance. At zero Vf-Vo is a maximum (minimum E^). 

becomes asymtotic to a value of 10 at about 10 -n . 

^ 0 

The maximum was not precisely defined by the experimental 
data of this program as there were no tests conducted with 

V --V 

-1 2. less than 0.27. However, the fundamental physics of 

0 
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turbulent mixing, dating back to Prandtl's mixing length 
(lypothesis. , is believed to be modeled by a criterion that 
mixing is minimum at Vf - Vg ~ 0. Therefore, the Kg curve 
was defined with a zero slope at Vf - Vo = 0, which was 

V.-V V.-V 

used to extrapolate K from — -- = 0.27 to -r — -°= 0. 

0 ^0 

The effects of an increased shear area oxidizer tube can 
be accounted for by use of the oxidizer hydraulic diameter 
Dh. Increased shear area (smaller hydraulic diameter) when 
compared to an equivalo'it area circular element (Deg) is in 
the direction of increased performance. This correlation 
treats increased sliear as a multiplier to the t\u. ic Vf-Vg/V 
relationship. Tfiat is. Mu' <uldition of added d i i iii-. i on 
surface area has the same offoit as increasing iiio ri'Ialive 
velocity. In practice however, increased shear designs an* 
difficult to fabricate and a better design choice would be 
to increase the relative velocity between the propellants. 

Swirl er Elements 


The relevant design parameters and curves necessary to esti- 
mate the Em of a swirler (oxidizer swirl) coaxial element are found in Figure 66. 
It must be noted that this is an equation valid for a swirl element with a tan- 
gential to axial velocity ratio of 1 only. The equation correlates the bulk 
of the swirler cold flow data to + 2%. 


E 

m 


100 


Cln (j^^p ) 
0 


( 8 ) 


Other values of swirl must be scaled using the empirical relationship defined 
in the upper right hand portion of Figure 66. Because the swirler element is 
such an efficient mixing element*, the E^ data generated during this program 
was all over 90%; a fact which limited the extent of the correlation. How- 
ever, the single element combustion data discussed in Section V was as useful 
as the cold flow data. By backing out E^ from ERE, the combustion data was 
used to supply some of the information necessary to complete this correlation 
The cold flow data was correlated with the constant "C" in the correlation 
equation equal to 2.7. However, using the 2.7 value to predict performance 
resulted in consistent overpredictions . A better combustion correlation was 
obtained using C = 5. Therefore, C = 2.7 should be used for cold flow mixing 
analysis and C = 5 should be used for the combustion problem. In addition to 
the above the following information is relevant. 

*T’rTco 1 d f Tow ~olTly~ - combustion influences the mixing of this element more 
adversely than any other element tested. 
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is a function of L/Dq. Again, long chanter lengths and 
small element sizes are in the direction of high performance. 

The degree of swirl does influence the resulting mixing per- 
formance. The curve shown in Figure 66 is a crossplot of 
the data shown in Section IV. Em for values of Vt/Va less than 
one were computed by applying this factor to the Em calculated 
through use of the swirl er equation. 

Inspection of the variable, Ks reveals a maximum { min Em) 9t 
a Vf-Vo/Vo equal to ten. It can be reasoned that values to 
the left of the maximum are controlled by swirl and to the 
right of the maximum by shear mechanisms. Since the total 
swirl component is directly related to the amount of oxi- 
dizer the mixing process is controlled by the oxidizer 
flow rate. Low values of Vf-Vo/Vo are generated by high 
oxidizer flow rates and thus swirl controls the mixing pro- 
cess. At high values of Vf-Vo/Vo less of the total flow 
is swirling, (less oxidizer), however, the relative shear 
is increased and becomes the controlling parameter. 

Premix Elements 


The relevant design correlating parameters and auxiliary 
curves for both the premix triplet and pentad elements are found in Figure 67. 
The design equation is 


100 - K In (t4tU) 

r P . fql 

'"m " 1 + 3 ( 0.07 W/Dq - 0.02)/(L/Djj) ^ ^ 

This equation correlates both the premix triplet and pentad cold flow data to 
within +4% Em using the cold flow Kp curve. Applying the cold flow Kp to the 
single element combustion data overpredicted the performance at low fuel flow 
rates (high mixture ratio). A Kp for combustion was found by applying the 
above equation to the single element hot fire data. Figure 67 shows the 

'^f ''f 

influence to be large (high Kp) at low - becomes identical to the cold 

”0*0 

flow at high values of momentum ratio. Because these elements were premixed 
i e., mixing occurs before combustion, it was expected that there would be 
very little corrtustion influence. However, a combustion influence can be reasoned to 
be a function of the bulk velocity in the mixing cup. As the H2 flow rate is reduced 
bulk velocity of the gases through the mixing cup is reduced increasing the 
likelihood that some combustion could take place inside the cup. Visual 
inspection of both the full scale and single element injectors did reveal that 
some combustion occurred locally inside the cup, which evidently modified the 
mixing process. Therefore the combustion Kp should be used for the combustion 
problem while the cold flow Kp should be used for cold flow predictions. 
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The design parameters and plots reveal the following addi- 
tional information. 

Em is sensitive to L/Dq. Long chamber lengths and small 
oxidizer (mixing cup) diameters are in the direction of 
higher Em- 

W V 

. An optimum design is achieved at a ^ — y— equal to three. 

”o ^0 

The width of the fuel injection orifice (w) is an important 
mixing process variable. When the fuel slot is narrow with 
respect to the oxidizer dia (0.33 W/Dq in Figure 67), the 
injected fuel penetrates the oxidizer and produces a well 
mixed product. When the fuel slot is wide with respect to 
the oxidizer diameter, the fuel jet is split by the oxidizer 
forming a striated stream that is fuel rich on the periphery 
and oxidizer rich in the core. 

F-O-F Triplet Element 

The relevant design parameters and curves are found in 
Figure 68 for the F-O-F triplet element. The design equation is as follows. 

E = 100 - K. In [9.5/(L/D )] (10) 

This equation, together with the appropriate Kt, correlates the data to within 
+ 5%. Referring to Figure 68 note that a combustion influence constant "C" is 
Applied to the momentum ratio to obtain Kt. The data plotted in Figure 68 was 
generated with "C" equal to one and is based on cold flow data only. However, 
as with the other elements discussed previously, combustion effects influenced 
the correlation. Based on the single element and full scale data, an optimum 
performance was obtained at Wq Vg/wf Vf sin u = 0.75; a value significantly below 
the value of two obtained in cold flow. A possible reason for this shift 
might be that the burning was occurring at the propellant impingement inter- 
face. The resulting change in temperature and density formed a hot gas barrier 
between the propellants analogous to the situation created by the liquid-liquid 
"blow apart" phenomenon found with hypergolic propellants. In order to pene- 
trate this barrier, relatively more fuel momentum was required resulting in an 
optimum at a fuel momentum higher than the equivalent cold flow situation. 

C = 2.3 correlates the combustion data, while C = 1 correlates the cold flow 
data. 


In addition to the combustion effects the design parameters 
and plots reveal the following additional information. 
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Performance increases with length and decreases with element 
size as with every other element. However, the length here 
should be based on distance from the impingement point, not 
from the face. This distance was used to correlate the data 

The F-O-F triplet element has a well defined optimum at a 
(Wq Vo)/wf Vf sin a = 2 for the cold flow situation and 
0.75 for the combustion case. Higher fuel momentum over 
penetrates the oxidizer jet resulting in a fan that is 
fuel -rich in the center and oxidizer-rich on the periphery. 
High oxidizer momentum results in the reverse situation - 
oxidizer in the center and fuel on the outside. 

Impingement angle affects mixing in that it defines the 
component of fuel momentum that is directed at the oxidizer 
jet. Increasing the impingement angle increases the fuel 
momentum available for mixing the propellants. 

Like Doublet 


The like doublet relevant design parameters and plots are 
found in Figure 69. The design equation is as follows. 


E = (K. - 


100 
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(L/D ) 
eq 


1.08 
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0.77 


0.07 

, L/D , 
L eq^ 


.36 I 


W V 
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(Ao/Af) 


1 .44 


] 


+ 0.99 (11) 


where: 


It should be noted that no combustion data was generated with this element and 
the correlation is based entirely on cold flow data. Considering the effects 
noted before it is likely that combustion would modify the mixing process. The 
following conclusions can be made. 

Lenqth and element size were the single most important 
variables affecting the mixing process. 

Inspecting the Kj curve, a weak optimum in mixing efficiency 
(high Kj implies high mixing efficiency) occurs at an element 
spacing (S/Do) of one. Element spaced either further or 
nearer than one oxidizer diameter, tend to result in lower 
mixing efficiencies. Large fuel to oxidizer spacing allows 
the propellants to dissipate before they mix. Small spacing 
does not produce better mixing because of interaction of the 
individual oxidizer and fuel fans. That is, the secondary 
impingement of the fans produces striated flow that effec- 
tively blocks further mixing. 
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Small momentum ratio and area ratio influences were noted. 

High values of fuel momentum and small fuel areas were in 
the direction of increased mixing efficiency. 

2 . T h e_ R elatio ns hip Between Combustion Efficiency and 

The work presented previously predicts for various elements 
as a function of the injection > flow, design and operating variables. Although 
Em Is a measure (and a very sensitive one) of an elements ability to produce a 
well mixed stream, it does not in itself predict performance, i.e., it is not 
contustion efficiency. In order to predict performance it was necessary to 
determine a functional relationship between E^^ and combustion efficiency. 

If the local properties of a gas flow field can be determined 
either experimentally or analytical ly, combustion efficiency can be calculated. 

A flow field can be subdivided into streamtubcs each of which has an associated 
mixture ratio and flow rate. A cotrtbustion efficiency for the total flow field 
can be computed by summing the mass weighted contribution of each stream tube 
and dividing by combustion efficiency based on overall properties. 


nn 

L 

i 


w.. c*., 

n n 




( 12 ) 


This calculation was performed for every cold flow data point tested in this 
program (see Table III). 

A calculation of nc is based on the local flow properties and 
Em as well is based on the same local properties (i.e., local mass flux and 
0/F). It follows therefore that a functional relationship exists between Em 
and nc- Plotting the values of Em against the corresponding values of nc 
(found in Table III) should define this relationship. This data isfound in 
Figure 70 and a functional relationship does exist. The two following items 
are restrictions to this relationship; 

a. This relationship is valid only for the H 2 /O 2 propellant 
combination. Other propellant combinations could result in a function much 
different than that represented in Figure 70. This difference stems from the 
relationship of C* with mixture ratio. H 2 /O 2 propellants have the character- 
istic of a flat C* vs 0/F curve while some propellants have a steeper func- 
tional relationship between these parameters. The insensitivity of H 2 /O 2 to 
mixture ratio distribution is demonstrated by the fact that an Em of 70% 
results in a combustion efficiency of 96.2%. Other propellant combinations 
at the same Em level may exhibit a lower combustion efficiency for an 
equivalent Em. 
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b. The distribution of mass and mixture ratio through the 
flow field is similar for each element and test condition; i.e., the same 
proportions of oxidizer rich and fuel rich streamtubes are found in each ele- 
ment streamtube. For the axisymmetric elements i.e., (the coaxial, swirl 
coaxial and increased shear coaxial) this relationship is quite valid. In 
general, the center of the flow field is oxidizer rich and progressi vely 
becomes more fuel rich towards the edge. As can he seen from Figure 70, the 
axisymmetric element data spread is less than the total spread of all the data 
and is more tightly grouped around the line identified as 0/F = 4 on Figure 70 
The external impinging elements, however, can not be said to have similar flow 
fields regardless of 0/F. For instance, depending on the relative fuel to 
oxidizer momentum ratio the centers of the flow field can vary between very 
oxidizer rich and very fuel rich. Fortunately this characteristic was experi- 
mentally determined as is shown on Figure 70 for the curves identified as 
0/F = 6 and 0/F = 2. This data indicates that the impinging element data 
tend to group above and below the 0/F = 4 curve, established from the coaxial 
data, for mixture ratios of 6 and 2 respectively. The curve in Figure 70 
therefore is the functional relationship that allows the designer to convert 
E|n to combustion efficiency. The data spread is such that this technique 
will yield a combustion efficiency that is usually accurate to within Z%. 

3 . Discussion of Combustion Influences 

The preceding discussion and data evaluation has been based 
to a large extent on element mixing data obtained during the cold flow experi- 
ments. It was noted in that discussion that combustion influenced the cold 
flow mixing, resulting in the need for modification of the correlations based 
on cold flow data alone. In order to explain this effect, the fundamental 
differences between cold flow and combustion must be addressed. This is done 
from an empirical standpoint in the following paragraphs and from a more 
rigorous analytical standpoint in the following section and Appendix A. 

The combustion process results in an overall flow field that 
appears to be "stretched" in length.* That is, when comparing identical tests 
0/F and mass profiles that have decayed to a uniform value in cold flow, still 
■ exhibit considerable nonuniformity under combustion conditions at the same 
length. Therefore, cold flow correlations based on a cold flow length when 
used to predict combustion performance would seriously overpredict that per- 
formance. It was necessary therefore, to develop a functional relationship 
between cold flow and combustion. 

Chamber length or combustion length is a parameter that can 
be scaled i.e., there exists a function that will transform combustion length 
to an equivalent cold flow length. Neglecting the influence of combustion on 

♦Se^Te^FiVn'TwhFr^^tTiFresults of identical combusting and noncombusting 
testing was discussed. 


- 149 - 



the basic mixing process and considering only the differences between cold 
flow and hot flow; results in the conclusion that the change in density 
between the two processes is the single biggest difference. The effect of 
combustion is to increase the temperature, decrease the density, and for con- 
fined flow increase the flow field velocity when compared to the equivalent 
cold flow case. 


The last effect is the key to generating a stretching func- 
tion. It implies that the common parameter between cold flow and combustion 
flow is stay time; i.e., a relationship between cold flow length and hot flow 
length can be equated on an equal time basis. 
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Equating time and noting that in cold flow the field velocity is a constant 
results in the following expression. 



cold 
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comb ^ 
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(15) 


Evaluation of the integral in Equation (15) requires a knowledge of how the hot 
flow velocity varies as a function of chamber length. The evaluated integral 
then defines the cold flow length as a function of the combustion length and 
the respective field velocities of the cold and hot gases. 

The key to evaluating the integral was knowledge of how hot 
gas velocity varies with length. A priori this information was not available. 
However, a number of reasonable relationships were postulated and various 
solutions of the above equation were found that related equivalent cold flow 
length to hot flow length. The validity of any one solution was determined 
by comparison of combustion efficiency calculated from the empirical equations 
at an equivalent cold flow length to the actual combustion efficiency at its 
nominal length. If any one solution provided a good correlation it was 
deemed a valid solution. 

Shown in the following sketch are 3 possible ways that velocity 
can vary with length. 
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The first solution Vconib = constant, assumes a step function 
change in density from cold reactants to combustion products at the injector 
fact. This solution is one that is a limit solution (i.e., it is the maximum 
deviation from the cold flow situation). The second solution assumes a linear 
distribution of temperature (velocity) with axial distance, and the third solu- 
tion assumes an exponential increase of temperature with distance. The 
equivalent cold flow length was determined by substituting for velocity in 
the right-hand side of Equation 15, each of the three possible velocity-length 
relationships described above and integrating the equation between zero and 
Lcomb- The constants in the various velocity relationships were evaluated by 
assuming that V = Vcomb at L = Lcomb and V = Vcold at L = 0. The resulting 
solutions expressed L(~oid 3^ 9 function of Lcomb, and the results are shown 
in Figure 71 for the 3 assumed velocity-length relationships. 

Inspection of these results reveal that a factor as large as 
3.7 in length exists between hot fire and cold flow conditions. This means 
that if a confjustion chamber is 37 cm long, the combustion efficiency pre- 
dicted from cold flow correlation should be based on a 10 cm length. 

In order to test this model, combustion efficiencies were 
computed using the correlations of the previous section. These computations 
were done at various cold flow lengths for the test conditions and element 
types of the single element testing. The cold flow length equation that best 
matched the experimental observed combustion efficiency was then a functional 
equivalent of the hot fire length. These relationships are illustrated in 
Figure 72 for the single element F-O-F triplet, swirler, premix and coaxial 
elements. Also included in this figure is the analytical solution of Vcomb = 
constant from Figure 71. 

The Vcomb = constant solution results in the maximum pre- 
dicted change of cold flow length with hot fire length. Referring to Figure 72, 
note that without exception experimentally measured combustion efficiency is 
lower than the values analytically predicted by the Vcomb ~ constant limit 
correlation. Thus the combustion chamber length must be significantly longer 
than that predicted by Vcomb " constant analytical model. 
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Two reasons for this deviation can be offered; and their 
relative importance is a function of the element type. The first is that 
coni>ustion does affect the mixing process and, as can be seen from Figure 72, 
affects it to varying degrees depending on the element type. The elements 
least affected by the combustion process lie closest to the analytically 
predicted line based on Vcomb = constant. This element is the premix which 
is not affected for obvious reasons; the bulk of the mixing is accomplished 
before combustion can take place. The swirler element and triplet elements, 
however, depend on a mixing mechanism that is highly concentrated near the 
injector face. If combustion interfers with these processes, then the mixing 
is dramatically reduced and as a result longer chambers will be required to 
effect mixing by turbulent diffusion. The data of Figure 72 indicates that 
this is the case. The coaxial element is only slightly affected for the 
second reason, i.e., the mixing process is so distributed over the length 
of the chamber that the combustion does not have much impact on the mixing 
rate. The above discussion indicates nothing about the absolute mixing 
efficiency in a given length. Even though the triplet element is sensitive 
to combustion effects it is still a higher performing element in a given 
length than the coaxial element. 

The fact that each element type results in a different com- 
bustion pattern is unfortunate for it precludes the existence of a universal 
stretching function. The designer must therefore use the empirical function 
of Figure 72 to estimate the cold flow length appropriate for use in the 
performance correlations. For element types different from those tested the 
designer must use judgement or refer to the generalized prediction techniques 
discussed in the following section (VIII, C, Analytical Model). 

4. Empirical Heat Transfer Correlations 

Heat flux the injector end of the chamber was found to be 
very sensitive to injector type and operating conditions during this program. 
In fact, chamber heat flux near the injector' was observed to vary an order of 
magnitude depending on which injection element v/as used. The following para- 
graphs describe how these effects were modeled. 

A relevant chamber heat transfer parameter for design pur- 
poses is the ratio of experimental heat flux at the injector end to nominal 
heat flux based on conventional turbulent heat transfer correlations. The 
nomi nal i zati on parameter accounts for all combustion parameters (Pc» 0/F, 
chaiTfcer dia) that are related to the chemistry and thermodynamics of the 
corrfcustion process. The deviation of this parameter away from one is a 
measure of how the injector affects the chamber heat transfer. Values of the 
heat flux ratio higher than one indicates local heat fluxes higher than would 
be predicted based on conventional turbulent heat transfer correlations alone. 
Values of the heat flux ratio lower than one indicate a benign environment due 
to a relatively low heat flux. 
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Th£; relationship between the heat flux ratio and momentum 
ratio at the head end of the chamber, 1.9 cm, (0.75 in.) is shown in Figure 73 
for coaxial, swirl coaxial, premix and triplet elements. The bands shown in 
Figure 73 includes all the data from both the full scale and single element 
testing. The values of the heat flux ratio illustrated in this figure should 
be used to estimate local values of heat flux at the injector end of the cham- 
ber. Heat flux ratios for axial positions between the injector and throat can 
be estimated by interpolation. As can be seen from the data of Sections VI 
and VII the heat flux ratio tends towards one at the throat independent of 
injector type and varies reasonably linearly with length between X = 1.9 cm 
and the throat plane. The heat flux ratio at any length can be estimated 
using the following equation. 





X 



(16) 


Inspection of Figure 73 reveals heat flux ratios over three 
for the triplet and premix elements at momentum ratios of 0.2 - 0.4. It is 
important to note that all the elements tend to have a functional relationship 
between heat flux and momentum ratio that is similar to that already dis- 
cussed for combustion performance. The triplet element for instance has a 
maximum heat flux at the identical momentum ratio that the correlation equa- 
tions and experimental data would predict optimum performance (see Figure 68). 
In addition the shape of the local heat flux curve is similar to the analogous 
performance curve. The premix and coaxial element curves indicate the same 
relationships; i.e., high local heat fluxes at injection conditions that favor 
high performance. The significance of the above is that an injector designer 
can expect high heat flux with high performing injectors, or conversely 
designing for low wall heat flux will be in the direction of low performance. 

5 . E mpirical Model ERE Prediction 

In order to test the validity of the empirical model, it was 
used to predict the performance of the single element injectors of this pro- 
gram and the full scale injector data of this and other relevant APS sized 
thrusters . 


The model was used as detailed in the Design Handbook and 
outlined below in the following procedure. Chamber lengths were converted 
to equivalent cold flow lengths via Figure 1 of the Design Handbook. Mixing 
efficiencies were computed using the specific equations and related curves 
for each element. ERE was determined from the empirical relationship between 
ERE and Ejj, (Figure 70). 
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Predicted vs actual ERE for the single element data is shown 
in Figure 74. This correlation groups all the data to within +5%. Reflected 
in this + 5% accuracy is the inherent inaccuracy of the single element data. 
Because of the small size and large heat losses, the measured energy release 
efficiency is probably only accurate to + 3%. It can therefore be reasoned 
that the accuracy of the empirical model may be better than + 5/». 

A better indication of the accuracy of the model is shov/n in 
Table XVI. This table compares the experimental performance of various H 2 /O 2 
engines with the model predictions. The predictions were based on the nominal 
design point of the particular engine With the exception of the TRW engine 
the nonnnal design point was 0/F = 4, Pr “ P07 N/cm^. The TRW engine predic- 
tion was hasc'd on an ovc'rall O/i of b.J. Ihe two, zoned injector predictions 
were based on a mass wc'igtit(Ml prediction o! the barrii'r aiul i.oce s tre.uiitiilH's . 
As can be seen by refering to the table the model predicts ttu' pc'rtormance to 
within +1.4% total scatter. It is also significant that the model predicts 
the ranlT of the injector performance accurately. That is, it predicts the 
highest performance for the injector that was highest performing and predicts 
lowest performance for the injector that was lowest performing. 

The fact is especially encouraging considering the different 
variety of injectors included in the list. There are four triplets (all at 
different element design points*; two coaxial injectors, one designed by ALRC 
and one designed by Rocketdyne; and two premix injectors both designed at ALRC 
but with different element design points. Using this model, the injector 
designer has a powerful tool for predicting GO 2 /GH 2 performance. 

C. DISCUSSION OF RESULTS - ANALYTICAL APPROACH 
1. Introduction 


The preceding section has described the empirical approach to 
data correlation. In addition, a parallel effort was carried on throughout the 
program utilizing a more general analytical approach. The advantage of the 
latter is that it identifies design and operating parameter interactions and 
their influences on mixing rate which can not be accounted for using an empiri- 
cal approach alone. The physical insight of the mixing and combustion pro- 
cesses provided by the analytical model permits extension and extrapolation 
of the data acquired on this program to untested future design requi rements . 


*The TRW O-F-0 triplet injector was not modeled by the empirical correlation. 
The model is based on F-O-F elements rather than O-F-0 elements. An accurate 
performance prediction for this element was obtained by using fuel properties 
in the design equations wher'’ oxidizer properties are indicated and vice versa 
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COMPARISON OF EMPIRICAL PERFORMANCE 
PREDICTION TO MEASURED DATA 
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A description of the analytical model formulation requires a lengthy discussion 
and is therefore documented separately in Appendix A. This section summarizes 
how the analytical model was used to correlate the experimental data. The 
methodology to use the analytical model to analyze or design future gaseous 
injectors is covered in the Design Handbook. 

Based on the preliminary data correlations from the cold flow 
portion of the program, the analytical model was initiated by postulating prob- 
able physical mechanisms that might explain the observed data behavior. The 
selection of pertinent process variables were ascertained by an approximate 
mathematical mixing model formulation described in Appendix A and was based 
on such considerations as the momentum and continuity equations, Newton's laws 
of motion, equation of state for idealised perfect gases, boundary layer equa- 
tions, potential flov/ equation, etc. This step of the analysis was the most 
critical and most important step in the analytical data correlation. The 
mathematical analysis model suggested probable groups of primary design vari- 
ables and their interactions. Absolute quantitative predictability was not 
initially considered a primai'y objective. The latter was achieved empiri- 
cally by appropriate curve fitting of experimental data vs predicted analytical 
trends . 


It was recognized that modeling of fundamental mixing pro- 
cesses was necessary to obtain meaningful scaling criteria that could be 
extrapolated to untested design variables. In addition, it was recognized 
that basic differences betw'een cold flow vs hot flow, single element hot fire 
and multiple element hot fire mixing processes exist. Only by recognizing 
and properly accounting for differences in the fundamental mixing processes 
was it possible to accurately extrapolate to and predict hot fire mixing char- 
acteristics. It was for this reason that the emphasis had been placed on 
modeling and correlating mixing processes, rather than correlating injection 
element design parameters as v/as done in the empirical approach. 

The concept behind the model is basically simple. The model 
consists of defininq the characteristic potential flow element dimension from 
the functional relationship of design parameter interactions. A boundary layer 
approximation is then defined and integrated to analytically predict an axial 
dynamic pressure decay rate which is used to predict mass diffusion profiles. 
The model yields considerable understanding of the mixing mechanism, sheds 
physical insight upon experimental data and is suggestive of design modifica- 
tions required for element optimization. However, it is hindered by the dif- 
ficulty that the design parame-ter interaction differs according to element 
type as summarized in subsection 2 of Appendix A. 

2. Descri_ption of Generalized Analytical Model 

Specific cold flow and hot fire mixing models were developed 
during the program for each element type (.ested which successfu^y correlated 
all types of test data v/hich v/as experimentally generated on this contract. 
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However, to attempt to disseminate this detailed knowledge contained in six 
different cold flow and four additional combustion mixing models would be 
impractical. Therefore, based on the objective of reducing this vast body of 
knowledge into a practical form usable to the average design engineer, a 
simpler generalized approach was taken. This should encourage greater usage 
because of its relative simplicity and yet maintain its generality, while 
sacrificing some absolute accuracy. This approach was only possible because 
certain groups of non-dimensional design parameters were common to all ele- 
ments and had a consistently predictable impact on the model solution. 

Cross correlations of the results from the specific element 
models led to the development of only two sets of non-dimensional mixing 
correlation parameters. Bg for axi -symmetric elements and Bt for two dimen- 
sional elements. This generalized model has successfully correlated chemical 
composition mixing rates, 0/F; mixing efficiency, performance efficiencies, 
% C* and % ERE; and chamber wall heat flux, Q/A, for both cold flow and com- 
bustion operation test data generated with every element type evaluated on 
this contract. The choice of correlation parameters was selected on the basis 
of maximizing its physical significance for the benefit of the designer. Those 
parameters having maximum significance on performance improvement are visibly 
identified. A conscious effort has been made to minimize the "black box" 
syndrome often encountered in analytical models. Even if the designer does 
not immediately satisfy all of his technical objectives on his very first 
design attempt, at least the direction of change required and the quantitative 
magnitude of change will be identified for the desired improvement by inspec- 
tion of the numerical values of his Ba or Bt parameters. 


Two versions of non-dimensional mixing correlation param- 
eters, Ba and Bt, have been defined to characterize axi -symmetric and two 
dimensional elements, respectively. The general form of these B parameters 
are shown below for both the fuel and oxidizer jets emanating from their 
respective source in any given injector element concept. 


^L/Dho^ 
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(L/D ) riiMii 


0.66 


f [Eletn. Type] 


to r- 


^do 


Vo 


10 ' 


Pd 


eq 


(L/D^^) f [Elem. Type] 


’tf 


0.8 + 0.2 


\ 


4 (0/F) 

MW /MW. 
0 f 


C . (1 + Da.) 

Of 


10^ 


( 19 ) 


Pd^ Vf 

Pd" 
eq 

( 20 ) 


Physically speaking, these parameters represent a generalized 
element L/D ratio. It has repeatedly been shown for all element ty pes that L/D 
is the single most important parameter for correlating mixin g efficiency; al 1 
elements tested show improved performance at longer L/D's. 

An identification of specific element design parameters (type, 
size, hydraulic characteristics), propellant physical properties (viscosity, 
density, molecular weight, temperature, combustion characteristics) and engine 
operating conditions (mixture ratio, injection aP's, chamber pressure) uni- 
quely define discrete numerical values of Bg or Bt for both fuel and oxidizer 
elements. This information is sufficient to completely define everything of 
interest about the gas mixing distribution. The generalized solution permits 
bypassing all of the intermediate computation steps required by the various 
specific element models. 

To calculate the local oxidizer composition evaluate either 
Equation (17) or (19). For that numerical value of the correlation parameter, 

B, on the abscissa of Figure 75 find the local oxidizer-rich mixing ratio ni 
on the ordinate, using the appropriate curve (axisymetrical or two dimensional). 
Then the local oxidizer rich stream tube mixture ratio is as follows: 


0/F. = ^ (21) 

Similarly the local fuel composition can be determined by using either 
Equation (18) or (20) with Figure 75 to find n. on the ordinate for the 
appropriate numerical value of the correlation parameter, B, on the abscissa. 
Then calculate the local fuel rich mixture ratio. 

0/F. = n.j X (0/F) (22) 
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To satisfy the overall continuity relationship, the mass 
fractions contained in the fuel rich and oxidizer rich stream tubes are uni- 
quely defined so that the overall mixture ratio is the design engine 0/F. 
This can be done by defining an intermediate parameter 


(1 - n^) 


(23) 


from which it is possible to define the oxidizer rich stream tube mass fraction 
and fuel rich stream tube mass fractions by Equations (24) and (25), respectively, 


X. 

J 


U 

1 + u 


"l + (0/F)j. 
1 + (O/F) 


(24) 


X. = 1 - X, 
1 J 


(25) 


The local streamtube mixture ratios and mass fractions then 
enable mixing efficiency and performance efficiency predictions as follows; 


= 100 < 1 - X . 
m 1 J 


■(0/F). - (0/F)' 
1 + (0/F), “ 


- X. 
1 


(0/F) - (0/F), 


(0/F) LI + (0/F).J 


nm = 100 [X. . n. + X. . n-] 


(26) 

(27) 


%C 


X, . C* (0/F), + X, . C* (0/F), 

★ = 

C* (0/F) 


(28) 


.pop - ^--^so (Q/^)j;^ -/sp 
/oERE - 1 T07T) 


(29) 


sp 


The term E^ as defined in Equation 26 is a two streamtube ver- 
sion of the Rupe mixing efficiency parameter (Ref. 1). Although Em is a widely 
used single value parameter which is useful for relative ranking of element 
mixing efficiencies, detailed information about the specific mass distribution 
is required to ascertain a quantitative performance efficiency prediction. 

For example, if low nozzle area ratio sea level combustion efficiency data 
(%C* or %ERE) were available, the corresponding mixing distribution and hot 
fire Em would still remain unknown; and it would be impossible to analytically 
extrapolate the sea level performance to high area ratio vacuum performance. 
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To circumvent the above shortcoming in E^, the mass weighted 
mixing parameter, nm, defined in Equation (27) is proposed. The term nm is 
Tikewisfrsingle vSlued function between 0 to 1.0 like Em, but it can be used 
directly to estimate performance efficiency. First of all, assume nj nn - 
Then using Equations (21) through (25), the oxidizer rich and fuel rich mixture 
ratios and respective stream tube mass fractions can be calculated, '^us, 
experimental hot fire %C* and %ERE are available, nm.can be estimated empirically 
using Equations (28) and (29) through parametric variation. Then using th s 
empirically inferred value of nm, the vacuum performance can be analytically 
extrapolated to any desired nozzle exit area ratio by using the appropriate C 
or Isp vs 0/F characteristics. 

The above discussion, encompassing Equations (17) through (29), 
is a summary of the generalized analytical procedure. A brief description of 
two of the terms in Equations (17) through (20) which are unique 


be 


described, i.e., D^f or and f [Elem. TypeJ 

(17) through (20) are standard parameters 


to this model 
The remaining 
defined in the 


will 

terms in Equations 
nomenclature list. 

Daf or Dao symbols for the combustion influence parameter. 
The combustion influence parameter was developed to relate the 
thermodynamics of the combustion process to parameters that controlled the 
mixing process. The detailed development of this parameter is discussed in 
Appendix A Fnr cold flow the combustion inf luence parameter has _a numeric al 
value of zero. For the G0?/6H? propellants tested on this program the com- 
bustion influence parameter varied from 0 to 10 as shown in Figu^ 76 If 

hot-fire performance and/or chamber heat flux values are required, then the 
coniiustion influence parameters (Da. and Da ) must be included in the calcu- 
lation For combustion performance problems several iterations may be 
required until the calculated combustion influence parameter converges upon 
the assumed value used in the equations for obtaining the parameter B as 
defined in Equations (17) through (20). This process is demonstrated in the 
sample problem in the Design Handbook. 

The second design correlation factor identified in the mixing 
correlation Equations (17 through 20) is dependent upon element type. This 
term designated "f [Elem. Type]" has a value of unity for shear mixing ele- 
ments without transverse injection momentum component such as ?r 

shear coaxial elements. The physical significance of this term is the ratio 
of the characteristic jet or fan dimension at the injector face plane to a 
downstream location. The mathematical derivation of th]S parameter is 
described in Subsections 2 and 4 of Appendix A. In addition, penetration 
mixing elements such as the premix or external impinging triplets are char- 
acteJized bj an initial composition at the injector face plfe which 
already pre-mixed to some degree as described in Equations (A-64) through (A-73) 
in Aopendix A. The initial premix composition significantly enhances mixing 
efficiency. For penetration mixing elements the adjustment for impingement 
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mixing is accounted for by determining an apparent L/D increase for an analogous 
shear mixing element. In addition, the swirl er co-axial element in a combustion 
environment requires special consideration as described in Subsection 8 of 
Appendix A. 


3. Analytical Data Correlation 

A correlation of analytically predicted cold flow Em's vs 
experimental test Em's are shown for various element families in Figures 77 
and 78. These elements were tested over the following ranges. 


Parameter 

F/E - newtons (Ibf) 
0/F 

L/D 

T - °K (°R) 

P 

a - Radians (degrees) 


Nominal Range 

67 (15) 13 - 222 (3 - 50) 

4 2-6 

1 0.5 - 2.0 

7 3-30 

300 (540) no - 300 (200 - 540) 

7t/3 (60) tt/4 - tt/2 (45 - 90) 


Of these tests, experimental data from the smallest 13 newton 
(3 Ibf) thrust element is least accurate. That is because all elements used 
the same measurement rake resulting in relatively coarser (fewer probe 
measurements) sampling with the smallest element. Furthermore, due to the^ 
large quantity of test data, both experimental data reduction and mathematical 
model analysis had to be highly automated. For example, element Cd’s were 
analytically estimated and input on the basis of orifice Reynold’s number, 

L/D, and assumed entry conditions. Only a few random tests were selected to 
compare predicted injection dynamic heads against experimental injector aP's 
to verify actual test Cd's. A more detailed input data selection could vary 
predicted Em's by a few percent. On some tests (see Table III) the integrated 
0/F differed from the injected O/F which may indicate some inaccuracy in the 
experimental Em. For some elements at short L/D's where the gradients are 
steep, the maximum fuel or oxidizer rich zones did not coincide with a probe 
position resulting in the loss of a high mass flux region between the probe 
positions. Thus, some disparity between the calculated and experimental data 
could not be avoided. 


The analytical predictions are in reasonably good agreement 
with the experimental data attesting to the validity of the generalized analyti- 
cal model. The analytical model is not only limited to its ability to predict 

absolute E . Even more important is its ability to identify significant design 
m 
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variables for a given type element and their interactions thus providing the 
designer an analytical means for optimizing his design or extrapolapng pre- 
vious test data to future design requirements. The physical insight which the 
mathematical model sheds on experimental data trends is likewise invaluable. 

By incorporating the combustion influence parameter into the 
generalized mixing parameters, Bg and Bt, the experimental single element h 
fire mixing performance data were predicted. As analytically predicted by 
the approximate mathematical model, combustion effects significantly impeded 
the mixing process and degraded performance compared to analogous cold flow 
conditions . 

The correlations for the single element hot fire test data 
using the generalized analytical model are shown in Figure 79. For the latter 
tests the cold flow GN 2 properties (MWq, po. uo) were replaced by hot fire GO 2 
properties. Othen^ise, the only difference between cold f ow and hot fire 
model calculations was the inclusion of the combustpn influence parameter, 
ba^ and Da to account for combustion effects. Satisfactory correlations 
befween predicted and experimental test energy release efficiencies were 
obtained for the shear coaxial, premix, and external impinging F-O-F triplets 
elements . 

The swirler coaxial element data were correctly predicted in 
their relative performance order (lowest to highest) but somewhat off on exact 
magnitude. The reason for this is as follows. As explained in ^ction 5 of 
Appendix A, the swirler coaxial element is characterized in a combustion 
environment with a collapsed oxidizer cone beyond which it is necessary to 
evaluate its mixing rate by axisymmetric shear ^^laracteristics. A 
swirler coaxial combustion model was developed consistent with actual physi 
wily Jc^Sm-ng mechanisms. The specific model predictions quantitatively 
agreed with experimental data in both magnitude and trend. Use of the modi- 
fications required to develop the generalized model, however, resulted in the 
(isULru? cr™ted. TlieVecomLnded procedure for utilizing the generalized 
cSrtustion correlation to describe the swirler coaxial element is described in 
Subsection 8 of Appendix A. 

In addition, the generalized hot fire swirler coaxial element 
predictions shown in Figure 79 were based upon using only ^1''^ 
stations (corresponding to the four measured thermocouple stations and the 
Szle th^aJ pUne) for the 5.1 cm (2;i"-).b''anfcer length tests e,g 
axial analvsis stations for the 14 cm (5.5-in.) chamber length tests. This 
coarse grid results in a bias which accounts for the consistent 
underorediction An alternative procedure would be to use a comparable coarse 
Vn^V^irler calculation technique but adjust the anticipated performance effi- 
ciency higher by using the empirical dependence between predicted and experi 
Sal performance shown in Figure 79 for all hot fire swirler elements. 
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4. 


Generalized Chamber Wall Heat Flux Correlation 


The analytical compatibility model described in the Design 
Handbook provides heat flux predictions in the absence of experimental 
data. Although the predicted absolute heat flux values are not as accurate 
as the performance mixing predictions, the model predicts relative trends 
and identifies the approximate heat flux magnitude within approximately + 30 
to 50%. 


The heat flux is calculated based on the properties of either 
the oxidizer or fuel rich streamtubes using the following relationship. 


0.026 




D 


T)T 


0.2 p 

P Lp 

n 0.6 


w. 


0.8 


(T 


gas 


\a 1 1 ^ 


(30) 


The designer must visually examine his injector pattern layout 
to identify the local wall mixture ratio. For example, if annular fuel coaxial 
elements are used the wall mixture ratio will be fuel rich. Hence the wall heat 
flux will be more closely approximated by using If external impinging 

F-O-F triplets are used, the wall 0/F will depend upon element orientation next 
to the wall and whether the fuel underpenetrates or over-penetrates the oxidizer 
jet. In the event of like doublets, predicted wall heat flux will depend upon 
whether a fuel doublet or oxidizer doublet is nearest the wall. Figures 80 and 
81 show correlations of predicted vs measured heat flux ratios for the various 
injectors tested on this program. 

The external impinging F-O-F triplet single element hot fire 
tests were underpenetrated above 0/F =2.0 resulting in a generally fuel rich 
boundary around two sides in line with the fuel orifices. However, two oxi- 
dizer rich lobes are squeezed out perpendicular to the fuel rich spray regions 
and somewhere in between, nominal 0/F zones occur. Analytically, one would 
predict widely varying heat flux ratios around the periphery of a single hot 
fire F-O-F triplet element. This thermal streaking tendency is verified by 
the measured heat fluxes experimentally measured around the single element 
F-O-F triplet as shown in Figure 81. The predicted minimum/maximum heat flux 
extremes predicted by the analytical heat flux model corresponding to fuel- 
rich, oxidizer-rich, and nominal 0/F stream tubes is in good agreement with 
experimentally measured thermal flux variations. 
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FlfiURF 65. COAXIAL AND TNCREASCD SHEAR COAxlA: DESION CORRF! ATTON:; 



Effect of Tangential Velocity or 
Swirler Performance 
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FIGURE 66. SWIRLER DESIGN CORRELATIONS 
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FIGURE 68. TRIPLET ELEMENT DESIGN CORRELATIONS 
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FIGURE 70. THE RELATIONSHIP BETWEEN H^/0 
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FIGURE 74. PREDICTED VS. MEASURED ENERGY RELEASE EFFICIENCY 
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IX. 


CONCLUSIONS 


The objective of this program was to develop an injector/chamber com- 
bustion technology base for gaseous hydrogen-gaseous oxygen rocket engines. 

The end product of this program was a design handbook that summarized the 
resulting information into models that characterize performance and compati- 
bility. These objectives were accomplished by: (1) empirically determining 

the mixing characteristics of a number of elements in a comprehensive cold 
flow test program. ( 2 ) empirically determining the effects of combustion by 
comparing the mixing characteristics of a swirl er element in a combustion and 
non-combustion environment. (3) conducting an extensive series of single 
element rocket tests to determine the performance characteristics and chairfcer 
heat flux of selected elements. (4) designing and testing four full scale 
injectors to obtain multielement performance and chamber heat flux information 
and (5) utilizing the data from (1), (2), (3) and (4) to synthesize models 
that relate performance and chamber heat flux to the design variables of the 
injector/chamber. 

The major conclusions from this program were: 

1 . Energy release efficiencies of 100% and 98.4% were demonstrated 
with F-O-F triplet and premix injectors, respectively, in chambers 7.62 cm 
(3 in.) long at nominal operating conditions. 

2. Chamber heat flux is highly dependent on injection element type 
with injector-end wall fluxes as much as three times that predicted by con- 
ventional turbulent heat transfer correlations. 

3. High performance 6 O 2 /GH 2 injectors can be unstable at high fre- 
quencies and provision for damping devices should be incorporated into the 
injector and/or chamber design. 

4. The zoned or barrier cooling concept results in a high performance 
penalty in comparison to GH 2 film cooling to effect a chamber wall temperature 
and heat flux consistent with long life thrusters. 

5. Mixing characteristics can not be modeled from cold flow data 
alone; combustion is a first order effect on the mixing process. 

6 . Two models were developed that allow an injector designer to 

a priori estimate the performance and chamber heat flux characteristics of 
potential injector concepts. 
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APPENDIX A 


ANALYTICAL MIXING MODEL DESCRIPTION 



This appendix describes the analytical mixing model which was developed 
to predict the influence of design and operating variables upon gaseous cold 
flow and combustion mixing rates. To facilitate cross reference, the appendix 
is subdivided into eight headings. The first section below describes the 
inter-relationships of the remainder of the model. 

1 . Physical Basis for Analytical Model 

One of the first significant discoveries which occurred while 
analyzing the cold flow data was the observation that those injection element 
types which rapidly dissipated their injection dynamic heads tended to be high 
performing (rapid mixers). From this it was apparent that if the axial rate 
of dynamic head dissipation could be analytically predicted, it might be 
possible to analytically predict mixing rates, ALRC had previously predicted 
dynamic head decay rates by integrating the element potential flow velocity 
profile modified for viscous boundary layer effects. The methods for calculating 
element potential flow characteristic dimensions are developed in subsection 2 
and the boundary layer approximation in subsection 3 of this appendix. Some 
element types are characterized by decreasing potential flow spray fan thick- 
nesses. These elements, however, are limited in their spreading rate at a 
critical length to diameter, (L/D)*, described in subsection 4. Combustion 
effects and chemical species change on mixing rate are analytically developed 
in subsection 5. Subsection 6 utilizes a shear coaxial element for illustra- 
tion and mathematically predicts its mixing rate based on the above calculated 
parameters. The mathematical mixing model is then extended to describe all 
other element types tested on this program in subsection 7. Finally, sub- 
section 8 describes the special treatment required for analysis of swirl 
coaxial elements. 

2. Potential Flow Characteristics of Several Conventional Gaseous 
Injection Elements 


A primary requirement of the analytical mixing model is a 
quantitative description of the characteristic element fan or jet dimension 
as a function of element design parameters and injector operating conditions. 
This subsection summarizes characteristic fans for several conventional 
injection elements which were tested on this and other program or which may 
be of interest for future gaseous injection element designs. ALRC had 
developed extensive injection models to parametrically predict element potential 
flow dimensions as a function of variable design and operating parameters on 
other programs. Liquid spray models had been developed at ALRC to analytically 
predict and correlate experimental atomized drop sizes with significant element 
design and operating parameters. These liquid element spray models were 
modified to describe gas injection element flow fields. The results of these 
modifications are described below. 
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Shear Coaxial Element 


The parallel stream shear coaxial element has been used exten- 
sively with the O 0 /H 2 propellant combination on numerous programs. It is 
mathematically characterized by a central oxidizer jet having a circular 


diameter. 


where: 



‘4 A 


1/2 


(A-1) 


Dq = oxidizer jet diameter 

Aq = orifice area of a single oxidizer element 


The oxidizer jet is surrounded by a concentric fuel annulus having 
inner diameter Df ] , and outer diameter, 0^2- The fuel jet is characterized by 
its hydraulic diameter- 


hf “ ^f2 ” ^fi (A-2) 

Both oxidizer and fuel potential flow characteristic dimensions 
remain constant as a function of increasing axial length. 

Increased Shear Coaxial Element 

A minor variation of the shear coaxial element is the increased 
shear element. The increased shear element is achieved by convoluting the 
oxidizer jet so that it is non-circular. Thus for the same oxidizer orifice 
area (indicative of constant thrust/element) its wetted perimeter is greater 
than the circular shear coaxial element, resulting in a smaller hydraulic 
diameter. Similarly, the fuel wetted perimeter across which turbulent mixing 
occurs is also increased resulting in reduced fuel hydraulic diameter. 

Both oxidizer and fuel hydraulic diameters are calculated using 
the standard equation. 

n = A (Area) 

^ h Wetted Perimeter ^ ‘ 

Swirl er Coaxial Element 

The swirl er coaxial element is characterized by a tangential 
oxidizer injection velocity component. The tangential velocity may either 
be imparted by a hydraulic swirler or a mechanical swirler. 

Three important characteristic design parameters of the swirler 
coaxial element are its oxidizer diameter. Dp-, discharge coefficient, Cj^q*, 
and oxidizer spray cone angle, Oq. Due to the tangential velocity component, 
the oxidizer stream adheres to the oxidizer element wall resulting in an 
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annular oxidizer jet at the injector face characterized by the oxidizer 
element discharge coefficient. The oxidizer jet hydraulic diameter at the 
injector face can be calculated from the continuity equation: 


^ho (L = 0) ~ ■ ^do^ ^o 


(A-4) 


The oxidizer sheet thickness at the injector face is one-half the above 
hydraulic diameter. 


^0 (L = 0) 



(A-5) 


The tangential velocity component results in a hollow oxidizer 
cone of angle, Gq, in the absence of fuel injection. Because of the annular 
axial fuel injection momentum impinging on the oxidizer cone, the resultant 
bi propellant cone angle is reduced. The resultant spray cone angle, ()> , can be 
calculated as follows from the momentum equation. 




tan 


w V sin 0 
0 0 0 


Un + w V cos e 
f f 0 0 0 


(A-6) 


It is readily evident that the bi propellant cone angle depends upon the element 
0/F and design A^/Aq ratios, as well as the oxidizer spray cone angle. 


Because of the divergence of the spray cone, momentum and conti- 
nuity considerations dictate that the oxidizer potential flow thickness 
diminishes with increasing displacement from the injector face as shown below. 


^0 _ ^do 

Dq " 4 [1 + 2 (L/D^) tan 


(A-7) 


Similarly, the fuel potential spray 
follows. 


t 


f 

D 


(L) 

f2 


^df - 
4 [1 + 2 


fan thickness can be shown to vary as 


11 . 

^ 2 ] 


'f2 


tan $] 


(A-8) 


Equations (A-4) through (A-8) above can then be used to provide the necessary 
input for the analytical swirl er mixing model. 


Parallel Sheet or Showerhead Element 


The showerhead injection element is the simplest to fabricate and 
analyze. It merely consists of axial injection orifices. In the case of 
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circular showerhead elements the circular orifice diameter is the character- 
istic potential flow dimension. 

D . JUSir (A-9) 

For non-circular parallel sheet elements the hydraulic diameter is the charac- 
teristic dimension and can be calculated from Equation (A-3) above. In any 
case the showerhead element has a constant characteristic potential flow 
dimension for any axial length. 

Like on Like Doublet 

The like on like doublet element incorporates both impingement 
and shear principles. The impingement momentum is mathematically characterized 
to describe a gaseous potential flow fan similar to a conventional liquid 
propellant spray fan. A like doublet is mathematically characterized by equal 
stream diameters, equal density, equal injection velocity, and equal momentums 
from the opposing orifices. Thus, the impingement fan is syrrmetrical . 

Consider a like-on-like doublet having stream diameter, D ; 
discharge coefficient, C^; injection velocity, Vq; impingement half-angle, 
propellant density, p; and viscosity, p. The axial component of stream 
injection velocity is (Vq cos a). The component of injection velocity normal 
to the plane of impingement is (Vq sin a). This component is locally stagnated 
at the point of impingement on the stream axis, but from energy conservation 
must result in a component (Vq sin a) flowing radially outward from the impinge- 
ment point in the impingement plane. The axial component of injection velocity 
(Vq cos a) is conserved at the outer fan boundary envelope after impingement, 
from axial momentum considerations. 

This impingement process corresponds to a constant radial 
(Vo sin a) velocity source traveling in a constant axial (Vq cos a) velocity 
coordinate system. Thus, solving for the resultant fan propagation angle 

(j) = sin"^ (tan a) (A-10) 

defines the shape of the fan envelope. 

From continuity relationships 

Wt = 2p Cj (J Vq = 2p (L tan <^) • t V^ cos a (A-11) 

wh6rG i 

wt = total flow rate from both orifices 

L = axial fan location from point of impingement 

t = fan thickness at axial location, L 
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Therefore, the like double fan thickness decreases with 
increasing displacement from the impingement point as 

t(L) _ ^^d 

4 (L/D) tan (t> cos a 


(A-12) 


Thus, in principle the like on like doublet is similar to the swirler coaxial 
element described previously from the standpoint that its characteristic 
potential flow spray dimension decreases with increasing L/D from the injector 
face. The local fan thickness can also be functionally related to the orifice 
diameter, discharge coefficient, and impingement angle as shown in 
Equation (A-12). 


External Impinging F-O-F Triplet 


The triplet element is characterized by both impinging fan 
potential flow and boundary layer shear principles. The fuel/oxidizer orifice 
diameter ratio is an important parameter for unlike impinging element designs. 
Except for the Af/Ao = 2.0 design which has equal fuel and oxidizer diameters, 
the other F-O-F triplet designs tested result in partial oxidizer mis- 
impingement because of the smaller fuel diameter. The fraction of oxidizer 
which does not impinge with the fuel is calculated below. 

Om = (A-13) 


where: _i 

3 = TT - 2 sin [Df/Dg] 


(A-14) 


The two misaligned oxidizer lobes are each approximated by an equivalent 
showerhead orifice of apparent diameter, Dq. 


D' = D 
0 


0 



3 - sin g 

2tt 


(A-15) 


The next step in the F-O-F triplet element analysis is to evaluate 
the relative fuel penetration distance into the oxidizer jet to determine 
whether the F-O-F triplet should be approximated with an oxidizer rich core 
between two fuel rich fans or a single fuel rich core with a split oxidizer 
fan. 


From the momentum equation considerations the F-O-F triplet fan 
propagation angle normal to the impingement plane can be mathematically shown 
to be _ . — 1 



w^ sin 

(w^ cos a^) + w^ (1 - Om)^ 


(A-16) 
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In the case of underpenetrated jets (Y < 0.5 D^) from continuity 
equation formulations the fuel potential fan thickness can be calculated as 
indicated below. 


tf (L) _ ^ ^df 

” 8 (L/D)^ cos tan <p 


(A-17) 


Similarly, 

equation. 


the oxidizer fan 


t„ (L) 


thickness can be calculated by the following 


•n (1 - Om) 
8 (L/Dq) tan <{> 


(A-18) 


If the fuel injection momentum over penetrates the oxidizer jet 
(Y > 0.5 Dq), the comparable fuel and oxidizer spray fan thicknesses are 
indicated By Equations (A-19) and A-20), respectively. 


t^ (L) 


t. (L) 


TT C 


df 


4 (L/D)^ cos tan 4 ) 
TT (1 " Om 1 

TlT^ 


(A-19) 


(A-20) 


D 16 (L/D ) tan <p 

Deep Cu^ Premix Elemenfe 
The premix elements evaluated on this program consisted of a 
circular axial oxidizer jet with either two (triplet) or four (pentad) 
rectangular fuel jets injected normally into the oxidizer stream. The pre- 
mixing cup consisted of a circular cross section identical in area to the 
oxidizer jet and a moderate (approximately 1 to 5) cup L/D ratio. The deep 
cup provides constraints for the premix gases and directs the gas jet in a 
nearly circular axial jet. For the premix injection element it was assumed 
that the mixing cup diameter is the characteristic potential flow dimension 
for both fuel and oxidizer jets. 

iU2 


D^ = D, 


4 A 


cup 


(A-21) 


Unlike the other impinging elements such as the like doublets, 
swirler coaxial, or F-O-F triplet, the characteristic potential flow dimension 
does not diminish with increasing L/D from the injector face. In the event 
that shallow cup (L/D < 1 or Acup » Aq) premix injectors are evaluated in the 
future, its characteristics could be expected to behave more like the external 
F-O-F triplet (or pentad) model described previously with = tt/2 radians. 


Other Element Types 

Based upon the above illustrated analysis concepts utilizing the 
injection momentum and continuity considerations, it is possible to synthesize 
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potential flow models for any conceivable injection element type based upon 
the physical insight of the injection process. For example the unlike 
doublet (which was not evaluated on this program) can be synthesized with 
only minor modification of the like doublet and F-O-F triplet element models. 

3. Boundary Layer Approximation to Predict Axial Velocity Head 

Decay Rate 

ALRC had previously developed a mathematical model describing 
liquid injection element hydraulics in reference 14. In this analytical 
hydraulic model, potential flow characteristics and approximate viscous 
boundary layer solutions had been combined to describe many diverse hydraulic 
phenomena including axial static pressure distributions. 

To reiterate the analytical objective, it was not the intent to 
develop a mathematically rigorous gas mixing model. It was only intended to 
utilize an approximate mathematical model to identify significant design 
variables and their interactions which could then be used to correlate the 
experimental data for numerical accuracy. Consistent with this limited 
objective. Figure 82 schematically illustrates the basic concept underlying 
the mathematical assumptions. 


The injection velocity profile at the injector face plane is 
assumed to be described by potential flow. It is readily recognized that a 
finite velocity profile exists at the injection plane depending upon feed 
system design parameters, but the uniform velocity approximation greatly 
simplifies the mathematical solution. During cold flow testing, the gases 
within the combustion chamber are relatively quiescent. Therefore, the high 
injection velocity gases from the injection element are decelerated by viscous 
drag. Strictly speaking the viscous boundary layer equations should be solved 
in cylindricdl co-ordinat6s for the jet with incredsing rddidi dnd dxidl 
static pressure gradients simultaneously with the momentum and continuity 
equations. Because of its mathematical complexity it was more practicle to 
solve the problem using boundary layer approximations to identify the 
general form of the solutions and use the empirical data to fine tune the 
model for numerical accuracy. For the gas injection elert«nts evaluated on 
this program the Reynolds numbers were on the order of 10i>; indicative ot 
turbulent boundary layer characteristics. For simplicity the velocity boundary 
layer growth rate as a function of axial distance was approximated from the 
flat plate solution. 


0.375 L 


Re 


0.2 


(A-22) 


For convenience the above equation was transformed to the following 
representation. 

6 0.75 (L/D)°-^ 

R o. 0.2 


(A-23) 
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To be exact, the injection velocity profile should have been 
integrated to simultaneously satisfy the momentum and continuity equations 
with appropriate boundary conditions. However, for mathematical simplicity, 
the assumption was made that the initial potential flow diameter is reduced 
radially inward at the same rate that the boundary layer grows radially out- 
ward. Then a reduced oxidizer potential flovj boundary coefficient, C^, can 
be defined by 

Similarly, an increased oxidizer boundary layer contraction coefficient, C^-, 
can be defined by 

Cc = 

Next, it was approximated that the local oxidizer velocity profile in the 
boundary layer varies linearly from VqJ st C[j to zero velocity at Cq as shown 
schematically in Figure 82. The momentum transfer (dynamic head dissipation) 
rate can now be calculated by integrating the dynamic head throughout the 
local velocity profile. 


Pd(L) = 


C ttR ‘ 
C 0 


^0 




Ro 

/ 

0 


(1/2 yj) 2.rdr 


1/2 % V 


1 - 


■ */^/J 


2Trrdr 


(A-26) 


Integration of Equation (A-26) results in the following solution. 



(A-27) 


The dynamic head dissipation rate in the near zone can be adequately approximated 
by neglecting the second term in the right handmember of Equation (A-27). This 
makes the mathematical solution computationally tractable for manual 
calculation. 
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(A-28) 


Pd(L) 


C, I /C . - 6 /R f- 

b _ do 0^ Q-* 

E + 6 /R 
do 0 0 


The above derivation has been limited to an axi -symmetri c element such as the 
oxidizer circuit of a shear coaxial element. 


For a two dimensional element such as the coax annular fuel jet, 
a parallel derviation is as follows. From subsection 2 it was shown that the 
fuel hydraulic diameter for a shear coax is 

“hf ' Of2 - “fl <''-29) 


where: 

"" outer fuel annulus diameter 
= inner fuel annulus diameter 

and the fuel gap thickness 


tf = 1/2 (0^2 " ^fl^ 


(A- 30) 


Equation (A-22) can then be transformed as follows. 
, 0.75 

to ~ 0.2 


Re 


hf 


(A-31) 


The comparable definitions of the two dimensional boundary and contraction 
coefficients are indicated in Equations (A-32) and (A-33), respectively. 


Cb = - 2 (6^/t^) (A-32) 

Cc - -t- 2 (6^/t^) (A-33) 

The dynamic head ratio associated with the two-dimensional annular fuel velocity 
profile described by Equations (A-32) and (A-33) was numerically integrated. 

The two dimensional solution is shown on Figures 83 and 84 along with the axi- 
symmetric solution. 

4. Minimum Potential Flow Fan Characteristic 


Of the element types evaluated on this program, the swirl er 
coaxial element, like-on-like doublet, and external impinging F-O-F triplet 
result in fan spreading characteristics. A mathematical description of the 
jet spreading rate is developed in subsection 2. 
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Swirler Coaxial Element 

The oxidizer spray fan thickness for the swirler is given as 


follows: 


do 0 


^ 0 ^^^ " 4 [1 + 2 [l/U) tan 


(A-34) 


The velocity boundary layer thickness can be approximated as shown below. 


0.75 (L/0„)“-® (D^/D^„) 


6(L) ■ 


vO.8 r 


ROr 


0.2 


'ho 


1 - /n^ 


do 


(A-35) 


The critical (L/Dq)* can be solved by equating twice Equation (A-35) to 
Equation (A-34) and approximating with the quadratic solution. 


Re 


0.2 


2 tan ML/D^)^-^ + (L/D^)°*^ - -f- 


-hfl- 


1 - 


/ I- 


(L/Do) - 


Re 


0.2 




'Ho 


1 - /T^ 


do 


Cjo tan i 


do 


1/2 

- 1 


= 0 
• 111 




4 tan (|i 


Similarly, solving for the fuel side, 


(A-36) 


[1 - (Pfi/0f;) ] °f2 
f2- 


tf(L) = — ; 

T 4 [1 + 2 (L/D..) tan <t>] 


0.8 


0.8 f2 


0.75 (L/D. 2 ) 0 


f2 ■ ^fl 
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^f2 ■ °fl 


Rer 
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hf 


(A-37) 


(A-38) 
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2 tan 0 ^ [1 + ^] [1 - ^] Re^ = 0 


"f2' 3 ‘-' ■ 0^2'* '■' ’ Df2'* '"'^ hf 


( r 


(L/D^2^ = 


i0.8 


-I 


1 + 


II 
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II 

V2 
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,1/2 n 


i.in 


Re. C .. tan -1 

°Hf J 


4 tan 4> 


(A-39) 


If the local L/D's exceed the critical values indicated by Equations 
(A-36) and (A-39) for oxidizer or fuel, respectively ; the oxidizer and fuel 
minimum fan thicknesses predicted by Equations (A-34) and (A-37), respectively, 
must be evaluated at the critical length. 


Like-on-Like Doublet 


A similar derivation can be developed for the like doublet. 


TrCd D 

" 4 (L/D) tan <j) cos a 


(A-40) 


6(L) = 
(L/D)* 


0.75 (L/D)°-® ,Dx 
0.2 ^ 2 ^ 


Re 


D 


= Re 


0.111 


D 


3 tan 4> cos a 


0.555 


(A-41) 


(A-42) 


External Impinging F-O-F Triplet 

★ 

Prior to calculating critical (L/D) 's for the F-O-F triplet, it 
must first be determined whether we are analyzing under-penetrated or over- 
penetrated conditions. 


If 


Y < 0.5 D 

0 


to ( l -) = 


'CjQ (1 - Om) P„ 
8 (L/D^) tan 4. 




The oxidizer velocity boundary layer can be calculated from Equation (A-41) 
above. Solving for the above set of equations. 
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6 tan <i> 


0.555 


(A-43) 


Th6 corr6spondi ng fusl side eguations are (A-41 ) and A-44) shown below. 

ttC 

t h) = f (A-44) 

8 (L/D^) cos tan <p 


0.555 


from which (L/D^) = Rep 


0.111 


^df 


, 6 cos tan 4> 


(A-45) 


On the other hand it can be shown that if Y > 0.5 D^, Equations (A-46) and 
(A-47) replaces (A-43) and (A-45), respectively. 

0.555 


(L/D^) = Re^ 


0.111 


”Cdo 


(L/Df)* = Re^^' 


0.111 


12 tan (() 

TlCdf 


3 cos a^: tan 


0.555 


(A-46) 

(A-47) 


5. Hot Fire Mixing Model Derivation 


Prior to this program 9 ALRC had recognized that combustion effects 
significantly impeded the mixing process. Elements which developed 90-100% 
cold flow Em within approximately 2.54 cm (1 in.) from the injector face were 
observed to deliver only about 70% Em or less under combustion environments 
with approximately 12.7 cm (5 in.) chambers. 

In recognition of this limitation a series of one atmosphere 
chamber pressure combustion tests designated herein as hot flow tests 
utilizing GHp/GO? (without gaseous helium) propellants was conducted with the 
low velocity swirl er coaxial element (see Section V). Other nominal element 
design parameters were F/E = 222 Newtons (50 Ibf), Af/A^ = 0.5, ~ 

T = 300°K (540°R). Mass spectrometer data were obtained at six radial 
locations at each of three axial stations to obtain local chemical composition 
data both with and without combustion. As expected on the basis of previous 
observations radical differences were noted in the mixture ratio profile with 
and without combustion. In support of this hot flow test data analysis, a 
general combustion model was formulated to describe the effects of combustion 
on the gas dynamics of the diffusion/mixing process. This portion of the 
analysis is general and applicable to any element type or propellant combina- 
tion tested. Results of this general model were then extended to apply to the 
GH 2 /GO 2 propellant combination and in particular to the swirl er coaxial 


- 193 - 



element which was hot flow tested. It was thereupon found that the 
mathematical cold flow model previously derived for non-chemical reacting 
elements must be modified as indicated herein to account for combustion 
effects on the mixing process. 

General Combustion Model Formulation 

The continuity equation can be written and differentiated as 

follows: 

w^ = PAV (A-48) 

(A-49) 

^t 

From the principle of conservation of total mass, dw. = 0. Therefore, from 
Equation (A-49), 


^ rM + 
p ■ " ^ A V-i 


(A-50) 


The equation of state for a mixture of assumed perfect gas approximation can 
be written and differentiated as follows: 


P = 



(A-51) 


^ ^ dP dMW dT 

P p ■ MW f 


(A-52) 


Equations (A-50) and (A-52) can be combined and transposed into the following 
form. 


P A V “ T ■ MW 


(A-53) 


The right hand member of Equation (A-53) characterizes the 
combustion effect of a gaseous propellant stream. The left hand member 
describes the gas dynamic result acting upon the reacting stream. Equation 
(A-53) is a basic derivation which is applicable to any element type and any 
gaseous propellant combination. (To be applicable to a liquid or solid 
propellant mixture, additional terms are required to account for the inter- 
mediate phase change.) 

The physical significance of the right hand member is as follows. 

If the combustion stoichiometry is such that the mole-weighted molecular weight 
of the product gas mixture is lighter than the mole weighted reactant gas mix- 
ture, its contribution is positive, and vice versa. This chemical effect is 
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a relatively small contribution in comparison with the thermodynamic effect 
reflected in the temperature term. All reactions of interest for rocket 
engines depend on propellant systems whose product gas (combustion flame) 
temperature is approximately one order of magnitude higher than the incoming 
reactant gas temperature. Therefore, the summation of the right hand member 
of Equation (A-53) is always positive and generally between 3 to 30 and 
depends only upon the propellant combination and 0/F, 

The physical significance of the left hand member is as follows. 

If combustion were to occur at constant volume, dA = dV = 0 and the static 
pressure would have to rise proportional to the chemical effect characterized 
by the right member. On the other hand, near the injector face plane the 
total propellant orifice area is considerably smaller than the chamber cross- 
sectional area. From the momentum equation dV = 0, and if dP ~ 0 then the gas 
flow area must increase; dA > 0. If low contraction ratio chambers are used 
it is probable that at some point the hot gas cross section may fill the 
entire combustor flow area. Beyond that point no recirculation exists, 
dA = 0 and dV > 0 due to flow acceleration dictated by the continuity equation. 
From the subsonic compressible flow equations, the increase in gas velocity 
must result in a reduction in static pressure, dP < 0. Even if chemical 
reaction is completed (right hand member of Equation (A-53) equal to zero), in 
the throat convergent section (dA ■ 0) the gas velocity must accelerate 
(dV - 0) and the static pressure must decrease (dP < 0). Incidentally, if the 
right hand member of Equation (A-53) is set equal to zero, the trivial 
solution which results, characterizes non-reacting cold flow conditions. If, 
in addition, dP = 0, the solution is characteristic of a free jet diffusion. 

In spite of the gross simpliciations made in the derivation of 
Equation (A-53), it is obvious that many additional real conditions of gaseous 
propellant systems can be described by Equation (A-53) depending upon the 
assumptions made concerning applicable boundary constraints. The inter- 
relationships of chemical reaction, heat release, static pressure, flow area, 
and velocity are linked through Equation (A-53). The necessary limiting 
assumptions and simplifications can only be readily obvious by considering 
the real physical mechanism under analytical consideration, i.e.. Equation 
(A-53) is a mechanistic equation. 

Application of Hot Flow Model to Swirler Coaxial Element Using 

GH 2 /GO 2 Propellants 

The following simplified stoichiometry is suggested to describe 
GH 2 /GO 2 combustion at an overall 0/F = 4.0. 

4 H 2 + O 2 - 2 H 2 O + 2 H 2 (A-54) 

A numerical example is presented that illustrates the combustion 
effect on the right hand member of Equation (A-53). In this example the 
reactants are characterized by the subscript 1 and the products by the sub- 
script 2. Assume initial reactant temperatures are at 300°K (540°R), 
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From which 

^ ^ (0 0/F = 4) 

n n n 

Similarly the right hand member of Equation (A-53) was evaluated 
in an identical manner at a multitude of other mixture ratio points ranging 
from 0/F = 0.1 to 1000. The numerical results are shown on Figure 76 of the 
main text. At 0/F = 0 (pure H 2 ) and 0/F = °° (pure O 2 ), the limits correspond 
to non-reacting cold flow and the right hand member vanishes. Figure 76 shows 
the right hand member varies from approximately 0 to 10. It is not surprising 
to find that the peak occurs at the stoichiometric 0/F. Thus incipient 
combustion effects will always tend to separate the unreacted fuel from the 
unreacted oxidizer to impede the mixing process. The method of incorporating 
the above calculated reaction data into the earlier cold flow mixing model is 
described in Section VIII of the main text. 

It has been shown that the non-reacting cold flow swirler coaxial 
model incorporates a diverging hollow cone oxidizer jet with an annular fuel 
stream. The initial injection momentums define the resultant hydraulic cone 
divergence angle. Incipient combustion must initially occur along this fuel- 
oxidizer interface. Viscous boundary layer drag initially decelerates a 


- 196 - 



small portion of the fuel and oxidizer causing initial species diffusion from 
the potential core. This causes minute concentrations of oxidizer and fuel 
to diffuse into each other forming a near stoichiometric boundary separating 
the two. Due to the momentum equation and in spite of the viscous boundary 
layer drag, dV is negligible very close the injector face plane. Empirically 
it is known that gases cannot sustain high shear forces or radial pressure 
gradients and dP/P] is on the order of 1% or less and is negligible. There- 
fore the energy release must be accounted for almost entirely by an increase 
in the injection gas stream area. 

As an initial approximation, the combustion effects on the Reynolds 
number were neglected since the calculated boundary layer thickness is only 
weakly dependent to the 2/10-th power of the Reynolds number. On the other 
hand, the heat release was assumed to have a first order effect on the 
potential flow gas stream characteristic dimension. Heat addition to and 
area increase of the hollow oxidizer cone relative to its resultant hydraulic 
momentum angle show that the rate of oxidizer flow area increase is greater 
than the rate of cone divergence causing the hollow oxidizer cone to collapse 
into a solid cylindrical jet within a short distance from the injector face. 

Thus at any given L/D from the injector face, the hot flow boundary layer 
thickness is comparable to the corresponding cold flow boundary layer thickness; 
but the hot flow potential flow stream characteristic dimension is considerably 
larger than the cold flow dimension. This results in a very slow rate of 
decrease in the hot flow boundary coefficient (Cb). Likewise, the local 
dynamic head to injection dynamic head ratio (pV /pV“) , decreases very slowly 
accounting for the slow rate at which the initial injection concentration 
approaches the equilibrium design value. In addition to the above features, 
the hollow oxidizer cone collapse reduces the wetted perimeter on which dif- 
fusion occurs from both inner and outer surfaces of the ho^ow oxidizer 
annulus to only the outer surface of tne collapsed cylindrical surface. This 
increases the hydraulic diameter, reduces effective L/D and further retards 
turbulent mixing with combustion. 

A high oxidizer swirl cone angle and low oxidizer discharge 
coefficient can delay the swirl cone collapse. This in turn improves perform- 
ance. Downstream of the collapse point the mixing rate is similar to the shear 
coaxial element. However, within the distance from the injector face plane to 
the initial swirl cone collapse point, the swirler coaxial element achieves a 
much higher mixing efficiency than a shear coaxial within the same length. 

Thus, it is necessary to calculate the longer equivalent chamber length which 
would have been required by a shear coaxial element to achieve the same 
mixing distribution as the swirler coaxial element at its collapse point. 

This length differenticl quantitatively represents the advantage of the swirler 
over the shear coaxial in a combustion environment. 


- 197 - 



6. Shear Coaxial Element Mixing Model 


Having established the above mathematical preliminaries, it is 
now possible to analytically predict and verify with experimental data the 
mixing rate for a shear coaxial element. 


Two fundamental types of cold flow data were measured at various 
transverse probe positions and variable axial stations on this program. A 
mass spectrometer measured the local mole fraction chemical compositions of 
H 2 (fuel), N 2 (oxygen simulant), and He (base bleed). From this data it was 
possible to determine the local mixture ratio profile. Pitot tube measure- 
ments of local dynamic pressure heads were also obtained to determine local 
mass flux at each probe position. 


For shear coaxial elements the oxidizer injection dynamic head 
at the injector face plane is essentially equal to the oxidizer injector 
manifold pressure to differential. At the injector face plane it is known 
that the initial cold flow oxidizer mole fraction must be 100% Ho and that 
Em = 0. Far downstream when the gas is assumed uniformly mixed (Em = 100%), 
and the local dynamic head is negligible, the equilibrium N 2 mole fraction 
can be calculated as follows. 

^N2 


^H^, eql 




w, 


MW, 


w. 


MW. 


(A-55) 


Figure 85 shows the experimental relationship between the measured 
dynamic head ratio vs the measured N 2 mole fraction on the element axis. It 
was noted from Figure 85 that the empirical rate of change of species mole 
fraction from its 100% pure initial value to its asymptotic equilibrium value 
was directly proportional to the rate of change in local element dynamic head, 
starting from its initial injection value to its asymptotic uniform negligible 
value. The above relationship indicated that the empirical mass diffusion 
rate is directly proportional to the momentum transfer rate. Based on the 
above data, it became clear turbulent diffusion is a force initiated mechanism 
which can be characterized better by Newton's equation of motion rather than 
by the molecular diffusion equation modified with a turbulent eddy viscosity. 
The driving force responsible for turbulent mixing is the radial pressure 
gradient potential. The dynamic pressure dissipation rate was analytically 
predicted by integrating the velocity head across the potential flow character- 
istic dimension described in subsection 2 of this appendix modified for the 
boundary layer calculated from subsection 3. Thus, it was postulated 
that if it were analytically possible to predict the axial rate of dynamic 
pressure dissipation that by inferrence it would also be possible to predict 
the mass diffusion and mixing rate. 
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The analytically predicted axial profile of oxidizer local dynamic 
head using the numerical solution from Figures 83 and 84 is plotted on Figure 
86 vs L/D for the 67 Newton (15 Ibf) shear coaxial element. The experimental 
oxidizer dynamic head measured on the element axis probe is shown for comparison. 
Figure 86 shows that, in spite of the many gross mathematical approximations 
which were made, the calculated dynamic nead decay rate based on a viscous 
boundary layer mechanism closely approximates the experimental trend. The 
predicted fuel dynamic head profile is also shown on Figure 86. It is obvious 
that the calculated fuel dynamic head dissipation rate occurs considerably 
faster than for the oxidizer. 


The experimental dynamic pressure downstream of the fuel annulus 
was not available for direct comparison with the analysis. This is because 
the midpoint of the fuel annulus occurs on a 0.092-in. radius and the nearest 
experimental probes are on a 0.125-in. radius. Near the injector face where 
the gradients are steep, the 0.033-in. radial mismatch is significant in 
comparison with the fuel annulus dimension (tf = 0.028 in.). Nevertheless, 
for lack of better data, the eight probes on the 0.125-in. radius were averaged 
as a measured lower limit of the fuel dynamic head. The correlations indicated 
in Figure 86 show that it is possible to analytically predict either axi- 
symmetric or two dimensional jet dynamic head dissipation rates based on a 
consideration of design parameters. 


The above dissipation rates have been treated independently as if 
each were acting alone. In reality, the two solutions must be coupled and the 
interaction of one must be considered on the other. The above equations are 
only valid upstream of the point where the oxidizer and fuel boundary layers 
initially come into contact. Downstream of this point, the two jets combine 
to act as one. The contribution of each to the resultant jet must be 
proportional to their respective mole fractions since mixing occurs on a 
molecular basis. Since the rapid mixing zone near the injector face is 
controlled by Newton's equations of motion, the primary normalizing parameter 
must be the injection dynamic head since this provides the mixing potential. 
Therefore, define an equivalent injection dynamic head, Pdg which is the 
overall result of the separate fuel and oxidizer contributions. 


Pd 



w„ Wx 

Q X Pd + — — 
MW 0 MW. 
0 f 


x Pd^ 


w„ w. 

■ Q _- + —f_ 

MW MU. 

0 f 


The above relationship can be 


simplified to 


(A-56) 
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(A-57) 
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It is then possible to analytically predict the axial mole 
fraction decay rate of both fuel and oxidizer as indicated below. 


" ^ 0 , eql Pd 
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eq 


^^0 initial ^o, equilibrium^ (A-58) 
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Pd. 


Pd. 


eq 


^^f initial ^ ^f» equilibrium^ (A-59) 


In general the cold flow interaction exponent (F) represents the 
ratio of the common fuel-oxidizer wetted perimeter available for inter- 
propellant mixing divided by the total fuel or oxidizer wetted perimeter 
available for diffusion and shear mixing. 


Bi Propellant Interaction Area 
Total Wetted Perimeter 


(A-60) 


For example, the shear coaxial element has Fq = 1.0 because the oxidizer jet 
is completely surrounded by the fuel annulus; but Ff = 0.5 because bipropellant 
interaction only occurs on the inside surface of the fuel annulus. 

The analytically predicted oxidizer mole fraction decay rate for 
a shear coaxial element using Equations (A-1), (A-2), (A-8) , (A-57), and 
(A-58) is plotted on Figure 87. The experimentally measured maximum oxidizer 
mole fraction which occurs on the coaxial element axis is plotted for compari- 
son. Again, the mathematically predicted N 2 composition profile provides an 
adequate first order approximation to the measured local compositions. Like- 
wise, the calculated fuel mole fraction profile is plotted on Figure 87, and 
the 0.125-in. radius probe measured compositions were averaged and plotted with 
the full realization that the actual maximum fuel mole fraction might be some- 
what higher than shown. 

It should be remembered that, near the injector face, the radial 
dynamic head gradient and species gradient are very steep. A deviation as 
small as one probe diameter from perfect concentricity between the element 
centerline and the probe measurement position could significantly influence 
the ifieasured data. Ihus, it is probable that the real (unmeasured) data may 
actually be in better agreement with the analytical prediction than the 
experimentally measured data shown in Figures 86 and 87. Therefore, it can 
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be concluded that the viscous boundary layer solutions provide an acceptable 
theoretical model of the shear coaxial mixing rate as a function of length. 


The experimental shear coaxial mixing data in Figures 12 and 13 
indicate high performance at low Af/Ag or low 0/F. Both cases result in low 
Pdo/Pdf ratio. From Equations (A-57) and (A-58) these influences are 
quantitatively predicted. Figure 88 shows a correlation of experimental 
vs the aforementioned oxidizer diffusion parameter prediction. In the range 
from 0 to 30 L/D at nominal (4.0) 0/F and 2.0 to 6.5 0/F at L/D = 7, a nearly 
linear correlation is obtained verifying the validity of the approximate 
mathematical mixing model. 

The analytical development described thus far has utilized the 
shear coaxial element as an example because many O2/H2 injectors have been 
developed using coaxial injectors. Another reason is that two basically 
different analytical solutions result for axi -symmetri c and two dimensional 
mixing characteristics and the shear coaxial incorporates both characteristics 
for oxidizer and fuel respectively. The fact that axi - symmetri c and two 
dimensional mixing characteristics differ has been experimentally verified 
by empirical turbulent eddy viscosity correlations which note that two different 
empirical correlation constants are required. 

7. Other Element Mixing Models 

The model derived to this point is valid for those element types 
whose potential flow characteristic dimensions remain constant with increasing 
axial length. For a summary of elements that satisfy this condition review 
subsection 2. 


Another series of element families exist whose transverse injection 
momentum components cause diverging fans whose characteristic potential flow 
thickness initially decrease inversely proportional to increasing L/D. 


t(L) !l_ 

dW l/d 


(A-61) 


or |/ 

ufel “ 1 + K3 (L/D) 

where: 

t(L) = characteristic potential flow dimension at length, L 

D(0) = characteristic potential flow dimension at the injector face 
plane 

K^. = appropriate design constants defined in subsection 2 
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Again, see subsection 2, to determine which elements satisfy this behavior. 
These latter element types dissipate their injection dynamic heads at a 
considerably faster rate with L/D because not only do their velocity boundary 
layer thickness increase proportional to (L/D)0<8, but their characteristic 
potential flow dimension from which 6 is subtracted to define Cb (Equations 
(A-4) and (A-12)) decreases . 

The characteristic potential flow dimension of these diverging 
fans do not decrease in thickness indefinitely with increasing (L/D), however. 
A critical (L/D)* exists for each element which characterizes its minimum 
allowable potential flow spray fan thickness. This critical (L/D) occurs 
when the increasing velocity boundary layer thicknesses from both sides of 
the spray fan just equals the local potential flow thickness. 

2 X r, [0 (L/D)*] = t (0 (L/D)*] (A-63) 


Specific application of this criterion to the various applicable elements is 
derived mathematically in subsection 4. The physical significance of (L/D) 
is as follows. The normal component of injection momentum causes the spray 
fan to spread perpendicular to the orifice plane. Lateral spreading of the 
spray fan causes the fan to become thinner from momentum jnd continuity 
considerations (see subsection 2). At the critical (L/D) when Equation (A-63) 
is satisfied the radial component of the injection momentum has been dissipated 
by viscous drag with the chamber recirculation gases and further radial fan 
sprGdding cGases. Thus (L/D)* rGprGSGnts the asymptotic limit beyond which the 
potential flow fan geometry essentially remains constant. 


The equations introduced thus far in conjunction with those from 
subsections 2 and 4 are sufficient to predict the axial mole fraction decay 
rate of any shear mixing family of elements. The shear mixing elements 
include shear coaxial, increased shear coaxial, swirler coaxial, like-on-like 
doublets, and parallel sheet (showerhead) elements. The shear mixing elements 
are characterized by having their initial fuel mole fraction start with IUU« 

H 2 and whose initial oxidizer mole fractions are 100% N 2 or O 2 . On the other 
hand, the penetration mixing elements (premix and external impinging E-O-r 
triplets) start with an initial concentration already premixed to some degree 
very near the injector face plane even with zero boundary layer growth and 
diffusion. The degree of premix effectiveness depends upon the design 
parameters as described below. 


In spite of the difference in numbers of fuel orifices per element, 
the premix triplets and pentads can be characterized similarly. The primary 
design parameter of the premix element was discovered to be its cup aP. In 
the designs tested on this program the mixing cup area was equal to the 
oxidizer area; and the fuel area varied between 0.5 to 2.0 times the oxidizer 
area. The measured cup aP was averaged between the fuel and oxidizer injector 
aP's less their respective injection dynamic heads entering the mixing cup. 
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(A-64) 
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The initial prettiix oxidizer and fuel concentrations at the mixing cup exit 
were calculated from Equations (A-65) and (A-66) below. 
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(A-65) 

(A-66) 


Notice the similarity between Equations (A-65) and (A-66) with Equations 
(A-58) and (A-59), respectively. Again, the similarities are based upon the 
premise that turbulent mixing is governed by Newton's equation of motion and 
pressure potentials. The initial concentrations at the injector face plane 
predicted by Equations (A-65) and (A-66) are then utilized in Equations (A-58) 
and (A-59) to determine local concentrations as a function of chamber length. 


In general, the cup aP calculated by Equation (A-64) is maximized 
at low 0/F which increased the total number of gas moles at constant weight 
flowrate because of lower fuel molecular weight; or at constant 0/F, cup aP 
increases at low A^/Aq ratios because higher fuel injection momentum normal 
to the axial oxidizer jet decreases cup discharge coefficient. 


The initial premix composition of the F-O-F triplet element 
depends upon the fuel penetration ratio into the oxidizer jet. For two fuel 
impinging diametrically opposed on a single oxidizer, a penetration ratio of 
one-half is optimum and results in highest performance. 


that 


In reference 15 it was found from empirical penetration studies 



Vf sin af 



(A-67) 


where: 

Y = fuel penetration distance into oxidizer stream 
Kp = empirical penetration coefficient 

The investigators in Reference 15 found Kp = 2.5 for gas/liquids and Kp - 3.2 
for GHe/GN 2 . 

Based upon the experimental test data from the current program, 
Equation (A-67) above was transformed to the following form for single element 
F-O-F triplet cold flow tests. 
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Xf(0) (0/F)p (A-73) 

^ MW /MW ~ 

0 f 

The above equations can then be used as applicable to predict the 
initial concentrations to be used in Equations (A-58) and (A-59) for external 
impinging F-O-F triplets. 

As part of the data correlation analysis, specific mixing models 
were developed for each family of elements which were cold flow tested. 

Element design parameters and test operating conditions were then input into 
the analytical mixing computer model to predict local dynamic heads, mole 
fractions, mixture ratio distributions, E^, nm. C* efficiency, energy release 
efficiency (ERE), and other hydraulic parameters for each test point summarized 
in Table III. 
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Although the details contained in each specific element model 
are too numerous to itemize, a cursory description of the basic equations 
will be summarized below to provide the reader with an understanding behind 
the model formulation concept. It was felt that to include all the details 
of each element model formulation would become so complex that it would 
discourage the reader from attempting to utilize the models for design or 
analysis purposes. In Section VIII of the main text, the solutions from 
these individual models are condensed into a generalized analysis technique 
which is applicable to any of these elements and in simple enough form to 
encourage its use as a universal gas/gas Design Calculation Procedure. 

For each element and test point, the following parameters were 
given or assumed as model inputs. For sake of brevity where both fuel and 
oxidizer parameters are applicable the propellant subscripts will be deleted 
and the equation will only be written once with the understanding that in the 
model, oxidizer and fuel were accounted for separately. The input parameters 
were propellant temperature, Tp; molecular weight, pressures, P; weight 
flow rates, w; base bleed mass fraction, WHe/(wQ + Wf + WHe)^ orifice areas, 
A; discharge coefficients, C<j; viscosities, p; interaction exponents, F, and 
axial length, L. 

The first step was to calculate propellant gas density. 

p . lip- (A-74) 

0 P 

Injection velocities and injection dynamic heads were calculated next. 


V = 



Pd = 1/2 P 


(A-75) 

(A-76) 


The hydraulic diameters at the injector face, L/D^,, and one dimensional 
Reynolds numbers based on diameter were calculated next. 


n = ^ (Area) 

^0 Wetted Perimeter 


and L/Dj^ 


(A-77) 




h 


(Cd °h 

P 


(A-78) 


Equilibrium mole fractions and equivalent dynamic heads were previously defined 
in Equations (A-55) and (A-57). Equations (A-23) and (A~31), based on injector 
face potential flow properties, were used to calculate the velocity boundary 
layer thicknesses for axi-symmetric or two dimensional element characteristics. 
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'V, div‘er( :nt fan two dimensional elements, if L/Dh £ (L/D)* use L/Di •, if 
/i:, • (L/l;)'' use (L/D)* to calculate t(L) at (L/Dj* (see subsections 2 

;i. ’ i ). Ffi' two dimensional characteristics 
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(A-79) 


C;, - 0 find axi -symmetric dynamic head ratio from Equation (A-28) or 
Lire 33 and two dimensional dynamic head ratio from Figure 83. If C 5 < 0 
-yKHiiii: head ratios from Figure 84. 


For penetration mixing elements, the initial premix mole concentra- 
, > defined in the front of this subsection. For all shear mixing 

t A; : j r t'Le initial mole concentration is 100%. Equations (A-58) and (A-59) 

L -tLii predict the oxidizer rich and fuel rich mole fractions for any 
ale rent and any cold flow test. 
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Further, define the following parameters 
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(A-81) 

(A-82) 

(A-83) 

(A-84) 


I rom which the mass fractions in stream tubes i and j can be calculated so 
that the overall 0/F is nominal. 


U [1 + (0/F)j] 

""j " (1 + UT [1 + (0/F)J 
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Based upon the above calculations the following performance parameters can 
be predicted. 


E = 100 
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8. Generalized Treatment of Swirler Coaxial Element with Combustion 
Effects 


Section VIII,C in the text described the synthesis of the specific 
element mixing models detailed herein into a simpler generalized mixing 
correlation procedure. In a combustion environment, the swirler coaxial 
element requires special treatment in the generalized model. The 
following procedure utilizing the generalized charts is recommended. Evaluate 
the hollow cone swirler for hot fire conditions using small increments of 
axial step sizes. Calculate the axial length at which the hollow cone 
collapses. 


"collapse 



1 - (1 + DaJ 


4 A 

0 

i 

2 tan <t> 


TT 


1/2 


(A-91) 


liotice that Da in Equation (A-91) is a function of L. To assure numerical 
accuracy it is recommended that at least ten points be evaluated between the 
injector face plane and the collapse point and evaluate the oxidizer composi- 
tion at the collapse point. Then parametrically analyze a similar shear coaxial 
design as a function of L and define a transition length, L-Prime, at which 
the shear coaxial oxidizer composition is identical to the swirler composition 
at the collapse point. Again, at least ten intermediate points from the 
injector face to L-Prime are recommended for numerical accuracy. Finally, use 
the hollow cone swirler combustion profile from the injector face to the 
collapse point and the shear coaxial solution evaluated at a pseudo-length 
defined by Equation (A-92). 


^Apparent Shear ^Physical Swirler ^ ^Collapse ^ 
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FIGURE 82. SCHEMATIC VELOCITY DECAY PROFILE FOR SHEAR MIXIN 
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FIGURE 83. THE EFFECT OF BOUNDARY LAYER POTENTIAL FLON 
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Local Dynamic Head Ratio, Pd^/Pd 
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FIGURE 84. LOCAL DYNAMIC HEAD RATIO FOR FULLY DEVELOPED 
BOUi’IDARY LAYER FLOW 
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FIGURE 85. EXPERIMENTAL RELATIONSHIP BETWEEN PRESSURE DECAY AND MASS DIFFUSION 
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FIGURE 86. DYNAMIC HEAD DECAY 
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FIGURE 87. MOLE FRACTION DECAY 
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FIGURE 88. ANALYTICAL CORRELATION OF EXPERIMENTAL DATA 




APPENDIX B 


COLD FLOW TEST AND DATA CORRELATION PROCEDURES 




1 . Cold Flow Facility and Test Procedures 

The basic means of generating mixing information is to measure dis- 
crete species and mass concentrations throughout the injector element gas flow 
field, A parameter termed, Em, relates these discrete measurements to a single 
number which indicates deviations from uniform mixture ratio. Any measuring 
technique must therefore be assessed by its ability to accurately measure local 
mixture ratio and mass flux. 

The cold flow test setup shown schematically and pictorial ly in 
Figure 89 consists of a twenty element probe and a gas supply and flow meter- 
ing system for accurate control of mass flow rate. The flow tests were con- 
ducted using hydrogen for fuel and nitrogen as an oxygen simulant. The need 
to scale different molecular weight propellants was avoided by using nitrogen 
as an oxygen simulant during the cold flow tests, and ignition of the exhaust 
gases was precluded. The cold flow hardware was used also for single element 
hot fire testing. 

Gas samples are taken with the twenty element probe which can be 
traversed axially and radially throughout the complete mixing zone. The rake 
shown in detail in Figure 10 has 20 total pressure and 3 static probes to yield 
mass flux data. The individual probes are arranged so that by rotating the 
rake, fine sampling of the flow field can be obtained at any radial and circum- 
ferential location. Composition samples were sequentially drawn through the 
total pressure probes and routed to the mass spectrometer via a 24 position 
valve. 


The individual gas samples were obtained with the sampling system 
that connected the total pressure tube of each probe to one of the 24 input 
ports of a 24-position single-throw selector valve (Scaniswitch, Mod W1261, 
qualified for 500 psia operation). The single output of this valve was con- 
nected to the input of the mass spectrometer. Operational experience with the 
mass spectrometer indicated that the time constant for the system is on the 
order of three seconds. Therefore, the sample selector switch was triggered 
every three seconds by a remotely controlled rotary stepping solenoid; and a 
complete scan of the injector flow field was accomplished in about 160 sec. 

The output of the mass spectrograph consisted of four continuous analog signals 
which were related to preset atomic mass numbers. Any four atomic mass numbers 
(AMU) may be selected, but for this testing AMU's of 2, 4, and 28 -- corres- 
ponding to hydrogen, helium, nitrogen — were preset. The hydrogen and nitrogen 
concentrations gave mixture ratio data. Helium was used to prevent gross face 
recirculation. 

Measurement of local pressure ratio and total pressure was required 
to determine local mass flux. To obtain the necessary precision for mass flux 
resolution, both total pressure and the difference between total and static 
pressures had to be measured. Therefore, in addition to the sample being con- 
nected to the mass spectrograph, the sample had to be monitored for both static 
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and total prossure. This was accomplished by connecting the output of the 
first Scaniswitch to the home channel of a Scanivalve. The Scanivalve, an 
integral 48-position valve and pressure transducer, permits very precise 
repeatable data by making up to 48 pressure measurements on the same trans- 
ducer. For this testing the Scanivalve used a +2 osi differential transducer, 
providing less than ^% full-scale error (+0.02 psi) and capable of withstand- 
ing 500 psi overload in either direction. 

The static pressure taps of the pressure rake are connected to the 
low side of the Scanivalve differential pressure transducer. To obtain a pres- 
sure reading, a solenoid valve was closed, thus putting total pressure on one 
side of the differential transducer and static on the other. A separate pres- 
sure transducer was used to obtain total pressure. As a system check, a third 
transducer measured the absolute static pressure. 

A high-accuracy digital data acquisition system was used to record 
data during these tests. It provided an on-line digital printout of test data 
in engineering units as well as a digital magnetic tape for detailed computer 
data reduction. The system has the required accuracy (j^O.01% four significant 
figures) at reasonable sampling rates. The accuracy requirement was very 
important, considering the relatively low velocity flow field and thus small 
differences in total and static pressure. The data acquisition system operated 
at up to 40 channels per second, with 30 channels of input. 

During flow test the propellants were set at a desired flow rate 
and mixture ratio. When the flow rates stabilized, chamber pressure was set 
using a remotely operated throat pintle. When chamber pressure was set, the 
automatic rake scan was started. At the end of the complete rake scan — 
which required about 120 sec -- the gas flow was stopped, the system purged, 
and the rake repositioned. Marker voltages were recorded that indicated 
Scanivalve, Scaniswitch, and isolation valve position. Sufficient scans of the 
data system were made to assure that steady-state pressure and concentration 
data are obtained. 

The same setup and procedures as described above were used in both 
the single element and multiple element cold flow testing. To permit use of 
the same rake for the multiple element testing, it was mounted off the center- 
line of the chamber, centered on one of the elements. Rotation of the aft 
closure relative to the injector permitted mapping the complete injector face. 

2. Cold Flow Data Acquisition and Reduction 

The initial step in reduction of a typical cold flow test was to 
create an Engineering Unit File (EUF) from the raw signal received by the 
digital acquisition system. The full data acquisition process is outlined in 
Figure 90. The EUF computed and created a master tape that contained the 
following information. 
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a. Local sample pressure 

b. Local sample composition 

c. Rake sample position 

d. System flow parameters 


This program adjusted the data for tare values of pressure, compo- 
sition and flow rate. Pressure data sampling was adjusted for transducer zero 
shift which was monitored continuously during the test. Composition data as 
recorded on the EUF was corrected by two factors: (1) the raw mass spectro- 

meter peak voltage reading (for a particular species) was reduced according to 
the average voltage residual recorded during the time in which pressure data 
were being sampled. That is, the mass spectrometer transmitted a small voltage 
reading (on the order of ten percent of a typical steady-state peak height 
sample) when it was not actually sampling gas, which was subtracted from 
regular composition voltage values; (2) the second correction involved adjust- 
ing the mole fraction data to maintain consistency with a known composition 
calibration gas, which was continuously sampled by the mass spectrometer. 
Essentially, the calibration gas was utilized to calculate the constants Ki , 

Kp, and K3 in the following equation. The equation is derived from the fact 
that the partial pressure of a particular species is proportional to its mole 
fraction in a gaseous mixture. 


X, = 


PrI ^ ^2 


_H1 

Ph2 + K 3 • P^3 


(B-1) 


where: 

x-j = mole fraction of species one 

P^1 = mass spectrometer peak height voltage for species one 

Pj^2 = "*3ss spectrometer peak height voltage for species two 

P^2 = >Ti9ss spectrometer peak height voltage for species three 

When both of the adjustments have been made, the composition data for each 
rake position were entered in the EUF summary table. 

System parameters such as weight flows, pressures and temperatures 
were averaged over the time period. All data was then stored on tape where it 
served as input to the analysis program Datnal . Following tape storage of the 
engineering unit file summary block, the data analysis program acquired the 
EUF data and computed the key test analysis parameters. It first calculated 
parametric test data based on system information and computed 24 non- 
dimensional parameters such as injection velocity momentum ratio, velocity 
ratio and delta velocity. These variables were based on the run conditions 
and injector geometry of the particular test. 
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Measured aP and static pressure matrices were computed. A complete 
flow and comoosition matrix was constructed by assuming flow field symmetry. 
Thus the data in the nonmeasured quadrants was constructed from 
measured quadrant. For some elements (doublets and parallel sheets )syi^try 
was based on a 180° section of the flow field rather than a 90 section 
Fiqure 10 of the main report shows the rake element locations. Probe 1 lies 
on the zero-degree boundary, while probe 8 lies on the ninety-degree boundary. 
The measured total pressure matrix was calculated by adding corresponding terms 
of the delta pressure and static pressure matrices. 

The local mixture ratio matrix was calculated with the equation 
shown below, from the mole fraction data obtained from the EUF program. 






MW 


(B-2) 


H2 


Local total mass fluxes were calculated based on subsonic compres- 
sible flow equations, where the gaseous mixture properties (ratio of specific 
heats, molecular weight) used to solve the equation were calculated based on 
the mole fraction data input from the EUF program. In essence, local proper- 
ties were obtained for the flowing mixture based on mass or mole weighting 
local nitrogen, hydrogen, and helium properties. 

The analysis of the engineering unit data was accomplished in a 
program termed DATNAL. After the local mass flux and mixture ratio matrices have 
been established from the EUF program, DATNAL was used to characterize the total 
charriber flow field. The total, fuel, and oxidizer local mass flux matrices 
were integrated over the entire chamber. The integrated mixture ratio through 
any radial location was calculated by dividing the integrated oxidizer flow 
rate by the integrated fuel flow rate. Typically, integrated flow rates do 
not equal measured flow rates when the entire chamber flow field is examined. 

Most tests were characterized by integrated flow rates that were higher than 
the measured value. The reason was suspected to be chamber recirculation 
zones. To characterize all cold flow tests on a common basis, key analysis 
parameters were typically compared at the radial location, R. R is the J'adial 
location at which the integrated mass flow rate (oxidizer and fuel) equals the 
measured flow rate. If the flow field integrates low, R is equated to the 
chamber wall radius (0,75 in.). 

The analysis parameters E^, and percent c* were computed. 
is the mixing coefficient, defined by the equation below. 
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is the central parameter in cold flow element modeling. It was calculated 
for the six flow regions formed by considering each rake probe_radial location 
to be the outside flow streamline. was also calculated at R by interpo- 
lating between the Ef^, values at the bounding radial positions. The R Em was 
the value which was usually utilized in intra- and inter-element correlative 
analysis. 


is a mass-weighted mixture ratio distribution term defined by 
the following equation: 
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(B-4) 


Effectively, nm indicates the mass-weighted deviation from nominal mixture 
ratio, nm is calculated through the same flow field boundaries as Em. Per- 
cent nominal c* was the final analysis parameter calculated in DATNAL. It was 
determined by calculating a local c* matrix based on local values of 0/F. The 
local c* distribution was mass -weigh ted, integrated and divided by the total 
integrated mass flow rate to produce a c* value characteristic of a particular 
test. Percent nominal c* was found by dividing test c* by nominal c* at the 
flow field integ^rated mixture ratio. Percent c* was calculated for the same 
regions and at R in similar manner to Em and nm- 
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APPENDIX C 


PERFORMANCE EVALUATION TECHNIQUES 





Summaries of the single element and multi-element hot fire test series 
performance are presented in Tables VII, VIII, XII and XIII of the main report. 
Included in these tables for each test are the test conditions (Pc, 0/F, and 
propellant temperatures), the one-dimensional equilibrium specific impulse and 
characteristic velocity for the test conditions, the delivered vacuum specific 
impulse, the delivered characteristic velocity based on Pc-1 , Pc-2, and thrust, 
and the delivered Isp, C* (Pc-1)> afid energy release efficiencies. The tech- 
niques used to calculate the various performance parameters are described 
herein. 


The chamber pressure values reported in the performance have been cor- 
rected to stagnation conditions, i.e. 


Pc-1 = (Pc-l)meas 


(C-1) 


P 

C 


-2 = (P^-2) 


meas 


(P^/P) 


(C-2) 


(C-3) 

(C-4) 
(C-5) 

The chamber pressure calculated in the above manner is equivalent to the throat 
stagnation pressure providing that no additional momentum pressure loss occurs 
within the combustion chamber between the Pc tap location and the nozzle throat 
and that parallel flow conditions exist at the pressure tap location. The 
Pc-1 pressure tap is located on the combustion chamber wall immediately upstream 
of the contraction section, while the Pc-2 tap is located on the chamber wall 
an inch or more upstream of the Pq- 1 tap. In reality, the P^-l static pressure 
is probably increased by the subsonic deceleration caused by the flow turning 
at the inlet to the contraction section. In addition, both the Pc-1 and Pc-2 
static pressure measurements should be corrected (lowered) for losses due to 
heat addition and friction downstream of the pressure tap location. The exact 
pressure correction, however, is a function of the gas properties, chamber 
contour, tap location, and the energy release profile established by the injec- 
tion and mixing process. The derivation and calculation of these correction 
factors was beyond the scope of this effort. However, the thrust based C* 
and energy release efficiency (ERE) are not affected by these considerations. 

The theoretical (ODE) specific impulse and characteristic velocity are 
computed using the ODE option of the JANNAF TDK computer programvio). The 
calculations are made on the basis of one-dimensional flow assuming complete 
chemical equilibrium. For these calculations, the oxidizer was assumed to be 


where: 

1 n y/(y-1) 
(P^/P) = (1 + ^M^) 

Y = f (0/F) 

M = f (y. eJ 
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comprised of 99.398% O 2 , 0.549% Ar, and 0.053% N 2 on a mass basis. The fuel 
was assumed to be comprised entirely of H 2 . The effect of mixture ratio, 
chamber pressure, propellant temperature, and expansion area ratio were con- 
sidered in evaluating the theoretical performance for each data point. 

The delivered specific impulse (Ispd) is calculated from the ratio of 
the measured thrust corrected to vacuum conditions (F + Pa Ae) to the total 
propellant flow rate. The delivered characteristic velocity is calculated 
by two methods. In the first method, c* is calculated by the conventional 
chamber pressure basis, i .e. 

c* = 

where Pr is the measured value (Pc“l 01 ^ corrected to stagnation condi- 

tions and At is based on the measured cold throat diameter. No corrections 
are applied for the nozzle throat discharge (two-dimensional) or boundary 
layer (displacement thickness) effects on c*. In the second method, the 
characteristic velocity is calculated on the basis of measured thrust and 
an analytical estimate of the delivered thrust coefficient (Cp). Thus, the 
thrust based c* is calculated as follows: 

c* (F) = g^/Cp^ 

where: 

^spd ~ ^vac^'^t 

^Fd "" (ODK)* ^"^DIV “ ^^BL/*"vac^ 

'^DIV ' Divergence Efficiency 

AFp, = Boundary layer thrust decrement from the 
nozzle throat to the exit 

The nozzle divergence efficiency (nniv) i^ evaluated using the charts from 
Appendix A of CPIA No. 178 (Ref. 17) while the boundary layer thrust decrement 

is evaluated using the turbulent boundary layer charts of Appendix B of CPIA 

No. 178. The thrust base c* is equivalent to the chanter pressure c* since no 
correction has been made for the chamber boundary layer and heat loss effects 
and the throat discharge coefficient effect. 

The specific impulse efficiency is the percentage ratio of the delivered 
Icn to the ODE Icn. The c* efficiency is the percentage ratio of the delivered 
c* based on Pc-1 to the ODE c*. Finally, the energy release efficiency is com- 
puted as follows: 
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ERE . (1.0 - ERL/l^p ((,[,£,) * 100 

where: 


ERL 


sp (ODE) ■ ^spd " KL ^^sp DL ^^sp BL^ 


Al = Kinetic performance loss 

sp KL 


Al 


sp DL 


= Divergence performance loss 


Aisp g|_ = Boundary layer performance loss 

The performance losses used to calculate the energy release loss (ERL) 
from the measured specific impulse were determined as follows: 


The kinetic performance loss accounts for the performance degradation 
because of chemical recombination lag during the gas expansion process. It 
is calculated using the ODK option of the JANNAF Two Dimensional Kinetic 
Computer program. Reaction rates used for the kinetic calculations were 
those recommended by the JANNAF Liquid Rocket Performance Subcommittee 
(Ref. 18). The kinetic loss was computed as the difference between the ODE 
and ODK specific impulse at the test conditions under consideration, i.e.. 


^^sp KL " ^sp (ODE) “ ^sp (ODK) 

Nozzle divergence loss is a measure of the performance which is lost 
due to non-axially directed momentum at the nozzle exit. The divergence 
efficiency used to calculate the I$p performance loss was taken from 
Figure A-1 of CPIA No. 178. The performance loss is then computed as follows: 

^^sp DL " ^^spd^"'DIV^ ■ ^spd 

The boundary layer loss (BLL), which accounts for the degradation in 
performance due to shear drag and heat transfer, was evaluated using the 
simplified design charts presented in Appendix B of CPIA No. 178. In the 
case of the single element test data where the heat loss was often substantial 
(0.5 to 12% of Isp). an additional heat loss term was included in the boundary 
layer loss calculation. The shear drag effect was included using the simpli- 
fied boundary charts assuming an adiabatic chamber wall. Next, the heat loss 
effect on performance was estimated by calculating the total heat transfer from 
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the chamber based on the heat flux data calculated from the thrust chamber 
thermocouple measurements*. The performance loss was calculated as follows: 


'^^sp ^sp (ODK) 

where : 

at = EQ/Cp 

Tjj = f (0/F) 

Cp = f (0/F) 



The total boundary layer loss is taken as the sum of the shear drag and heat 
flux components. The above procedure was used to calculate the boundary layer 
performance loss because it accounted for the test-to-test differences in heat 
transfer effects without expensive test-to-test TBL computer analysis. The 
procedure was checked for a few selected injector/chamber combinations and it 
was found that the boundary layer loss was within 5 % of the boundary layer 
loss value calculated using the JANNAF TBL computer program analysis procedure. 


* eQ = / q dA, where q (heat flux - Btu/in.2_sec) is evaluated from a curve 
fit of the 4 experimental values. Within the chamber contraction section 
the heat flux was assumed to be inversely proportional to the chamber 
diameter to the 1.8 power (q « D-1-8). For the -2 chamber, in which the 
last thermocouple measurement is 3.75 inches from the nozzle throat, the 
heat flux in the chamber downstream of the last thermocouple was assumed 
to increase to the theoretical value i.e., qtheo = DB * (Tn - 500)Dl*8. 
Since nearly 2/3's of the -2 chamber contained no thermal instrumentation, 
the heat flux and thus the resulting boundary layer performance loss for 
this configuration is much less certain than the B.L. losses calculated 
for the shorter -1 and -3 chanters. 
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APPENDIX D 


HEAT TRANSFER EVALUATION TECHNIQUES 




Sinnle and multiple element thermal data was evaluated using the same 
basic methods. Local heat fluy values were obtained and the ratio of local 
heat flux to nominal heat flux was computed. 


Values of local heat flux in the case of the single element testing 
were obtained using a data reduction computer program written specifically 
for this test program. Its essence this computer program performed a heat 
balance on an interior element of the chamoer, accounting for heat trans- 
ferred through the wall and heat conducted to or from adjacent elements. An 
internal surface heat flux transient was inferred based on the wall acting as 
a heat sink, with correlations for axial and circumferential conduction. A 
correction was also made for heat loss lg a constant temperature sink or to an 
external mass such as a flange. Thus. 

Heat Flux = (Internal Surface Area) ^ (Heat Storage Rate (D-1) 

+ Axial Heat Loss + Circumferential Heat Loss 
+ Heat Loss to Sink) 


In the case of the full scale 6660 N, (1500 Ibf) hardware, a different 
heat flux technique was used. The local values of heat flux were obtained 
from the wall temperature measurements by inputting them as boundary condi- 
tions to transient conduction models for the chamber walls. These transient 
models were constructed within the framework of the SINDA computer program 
(Ref. 19). Two dimensional and one dimensional models were utilized. Heat 
flux values were calculated by assuming a linear temperature gradient between 
the gas-side surface and the first interior node (a distance of about 0.060 in.). 
The equation for experimental heat flux is given as Equation D-2. 


exp 


exp 


6 


rad 

at 


R 


Al. 

Ar 


(G ,) aT 
rad' 

s 

2 

experimental heat flux (Btu/sec in. ) 

7 - p ^! y 

= conductance betv;een surface node 
and first interior node 

Temperature difference between surface node 
and first interior node 

Surface node area 

Average radius between surface node and first 
interior node 

Axial length of nodes 

Radial distance between surface node and first 
interior node 


(D-2) 
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This approach for calculating heat flux has been employed in previous programs 
at ALRC; however, only one-dimensional conduction models were used previously. 
Two-dimensional conduction models were used as much as possible to analyze 
data from the -1 and -2 full scale chambers. It was found that heat fluxes 
calculated with the 2-D model were about 15% lower in the throat region than 
the 1-D values. Smaller differences were noted upstream and downstream of 
the throat. 

The two-dimensional node network were determined for both the -1 and -2 
chani)ers, while a one-dimensional model was used for the L* section data. 

While processing the data, the situation sometimes occurred where one of the 
surface temperatures was not measured. In these cases, the radial string of 
nodes attached to the missing surface node was eliminated from the analysis 
by assuming it was insulated from the adjacent node strings. This assumption 
probably produced less accurate heat flux values and in some cases produced 
a one-dimensional analysis. 

Axial heat flux distributions and the effect of chatpber geometry, 0/F, 
Pc and propellant temperature on heat flux were evaluated in terms of the heat 
flux ratio defined below: 


Heat flux ratio = 


Measured heat flux 
Nominal heat flux 


The nominal heat flux values were calculated from a turbulent flow correlation 
for heat transfer coefficient (Equation D-3) which was then used in 
Equation (D-4): 


h = 0.026 T Pr°*'^ 

g a 


where hg is evaluated based on nominal combustion gas properties. 
Measured heat flux, = h (T^ - T^) 


(D-3) 


(D-4) 


where: 

T = theoretical combustion temperature 
0 

T = gas -side wall temperature 523K (960 R) (assumed to 
be representative) 

Equation (D-3) is a correlation for heat transfer coefficient in fully 
developed turbulent pipe flow. This equation is commonly known to ALRC as the 
DB equation and can also be stated in the more familar form shown as Equation 

(D-5). 
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(D-5) 



DB = the DB factor given as output by the THERMOCAL computer 
program (includes all physical properties and the 
appropriate constants) 

Experience at ALRC has shown that Equation (D-5) also yields reasonable 
predictions for heat flux within rocket thrust charrbers. In actual design 
practice, the 0.026 correlation coefficient is often adjusted to account for 
contraction ratio and injector streaking effects. The 0.026 value was used 
in this study because of convenience and because it yields a representative 
heat flux value which is adequate for comparing the relative heat transfer 
characteristics of injector elements. Use of a heat flux ratio for comparing 
the element heat transfer characteristics is preferred to the use of heat flux 
alone because the ratio accounts for the known effects of mass flux and 
physical properties variations (0/F and variations) in fully developed flow. 

The heat flux ratio parameter, groups some of the data reasonably well, 
but does not yield a completely generalized correlation for injector element 
heat transfer characteristics. The ratio does characterize the relative 
severity of the heat transfer environment produced by each element and shows 
the effects of 0/F and Pc variations over the above the effect normally 
expected due to overall changes in mass flux and gas physical properties. 

Axial distributions of the heat flux ratio were obtained for each injector 
element tested and are discussed in the main section of the report. 
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APPENDIX E 


LIST OF ABBREVIATIONS AND SYMBOLS 





A 

APS 

Ar 

a 

B 

C 

C* 

CR 

D' 

D 

Da 

d 

E 

ERE 

ERL 

F 

E-O-F 

G 

g 

H 

h 

I 

J 

K 

k 


Area, Alumel 

Auxiliary Propulsion System 
Argon 

constant in Rupe equation 

Mixing correlation parameter 

coefficient, chromel , constant, specific heat 

characteristic velocity 

contraction ratio 

apparent showerhead diameter 

Diameter 

combustion influence parameter 
differential 

Efficiency, activation energy, element 

energy release (combustion) efficiency 

energy release loss 

Fuel, thrust, interaction exponent 

Two fuels impinging on ox (Triplet Element) 

gaseous, conductance 

gravitational constant 

hydrogen 

heat transfer coefficient 
Impulse 

Injector related 
Empirical factor 
thermal conductivity 
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L 

L* 

In 

M 

ME 

MP 

MW 

N 

n 

Om 

0 

0/F 

O-F-0 

ODE 

ODK 

P 

Pd 

Pr 

Q 

R 

r 

Re 

R 

S 

SE 

St 


length, liquid 
characteristic length 
natural logarithm 
mass flux, Mach number 
multiple element 
mixing parameters 
molecular weight 
nitrogen 
number of moles 

Oxidizer niisimpingement fraction 
Oxygen 

mixture ratio 

Two oxidizers impinging on fuel (Triplet Element) 
One-dimensional equilibrium 
One-dimensional kinetic 
Pressure 

Dynamic pressure, = 

Prandtl number 
Energy Release Rate 
gas constant, radius 
radi us 

Reynolds number 

Bounding radius in Em calculation 
element spacing 
single element 
Stanton number 
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T 

TBL 

TJ 

TDK 

t 

U 

Y 

w 

X 

X 

Y 

ID 

a 

n 

Y 
A 
6 
e 

n 

n 

u 

0 

It 

p 


Temperature 

turbulent boundary layer 
Injector face temperature 
two-dimensional kinetic 
time, thickness 
mass distribution parameter 
Veloci ty 

premix fuel slot width 
flow rate 

axial location, mass fraction 
mole fraction 
penetration parameter 
one-dimensional 
impingement half angle 
non-dimensional parameter 
ratio of specific heats 
differential 

boundary layer thickness 
area ratio 

efficiency, mixing ratio 
ERE at area ratio one 
viscosity 

oxidizer spray cone angle 

3.1416 

density 

sutmiations 

heat flux, fan propagation angle, resultant spray cone angle 
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Subscripts 


a 

c 

cold 

comb 

d 

e 

eq 

eql 

exp 

f 

g 

h 

i 

ii 
in 
int 

j 

1 

m 

n 

n 

nn 

nom 


axial, average, axisymetric, ambient 
chamber, coaxial, combustion 
refers to cold flow length 
refers to combustion chamber length 
doublet, discharge, diameter 
exi t 

equivalent 

equilibrium 

experimental 

fuel 

gas or gas-side 
hot* hydraulic 

refers to fuel rich streamtube 
streamtubes 

injector related, injected 
integrated 

refers to ox rich streamtube 

local 

mi xing 

number of fuel rich streamtubes, nominal 
number of oxidizer streamtubes, nominal 
number of streamtubes 
nominal 
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0 

ODK 

P 

s 

sp 

T 

tcv 

t 

V 

wall 

X 


oxidizer, overall 

one-dimensional kinetic 

propellant, premix, at constant pressure 

swirl er 

speci fic 

reci rculation 

thrust chamber valve 

tangential, total, throat, triplet, tank, two-dimensional 
venturi 

chamber wall properties 
axial location 
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I. 


INTRODUCTION 


The purpose of this guide is to present succinct design procedures for 
gaseous propellant injectors. Methods are presented for estimation of combustion 
efficiency, chamber heat flux and stability characteristics. Except for 
thermochemical data, which are readily available, the design guide contains all 
required procedural and empirical data for making preliminary design calculations. 
The design procedures are derived from a data base consisting of 405 cold flow and 
fire tests with six basic injector element types. Thus, the designer is assured 
reliable injector and chamber preliminary designs and associated performance 
characteristics. 

The design procedures are based primarily on data from NASA Lewis Research 
Center Contract NAS 3-14379, Investigation of GH 2 / 6 O 2 Combustion. However, 
industry wide data sources for gaseous propellant data were consulted and 
utilized in derivation of the design procedures and correlations. 

The experimental portion of the subject contract was restricted in scope 
to tests using gaseous hydrogen/gaseous oxygen propellants under test conditions 
representative of auxiliary propulsion systems, APS, for the Space Shuttle. 

However, from the start, the motivation for this program was a recognition 
of the fact that no significant design procedures were available for injectors 
using propellants which are both in the gaseous state. Therefore, the intent 
during the development of design approaches on this contract has been to make 
them as applicable as possible to all gaseous propellant injectors. 

It became apparent that this could be done most effectively by pursuing 
two approaches to gaseous injector modeling. The first approach had the objective 
of using the test data and correlating it directly with injector/chamber design 
parameters which are recognized from both theoretical and empirical standpoints 
as the controlling variables. This approach was defined as the Empirical Model 
and is located in Sectionll of this document. It has the advantages of 
( 1 ) inherently being the most accurate procedure for gaseous injectors which are to 
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Introduction (cont.) 


be designed within the operating envelopes and using the same propellants as this 
program, and (2) simplicity in the calculation procedure itself. Both of these 
advantages accrue from the process used to develop the model, i.e., it was de- 
rived directly from the test data and only correlated with the sensitive design 
and operating parameters. However the Empirical Model has a corresponding dis- 
advantage; it lacks generality since it stressed utilization of the test data 
and did not concentrate on quantifying the mechanistic causal relationships 
of the mixing combustion process itself. Therefore it should be used with 
caution when the desired propellant combination or operating envelopes differ 
significantly from those used on this program. 

The second approach had the objective of understanding the mixing/com- 
bustion process to the maximum extent possible, using both available theoretical 
knowledge and new techniques suggested and developed from close observation of 
the test data. This approach was defined as the Analytical Model and is located 
in Section III of this document. It is somewhat more complex than the Empirical 
Model, but it has quantitatively characterized the mixing/combustion process for 
gaseous propellants so that it is general in nature and can handle all gaseous 
propellants and operating conditions. To make this procedure practical many 
second-order physical and mathematical influences had to be ignored. However, 
the procedure is based on sound principles and accurate modeling of the first 
order influences. It places at the disposal of the designer a preliminary design 
tool that is broadly applicable and should be particularly valuable in trade-off 
studies that need to consider both the influence of design parameters and propellan 
combinations. Further it provides the foundation that can later be calibrated, 
if required, in any particular region of interest encompassing gaseous propellant 

thrusters. 

Both the Empirical Model and Analytical Model are described in detail in 
the main volume of the final report. This design handbook provides only the 
procedural steps, with most of the required information displayed in charts, 
graphs, and tables for clarity of presentation. The figures and tables applicable 



I Introduction (cont.) 

to each methodology are located at the back, of their respective sections. 

An example problem is included for each procedure, so that the reader can be 
assured that he is using the model correctly. The nomenclature is listed at 
the end of the handbook and comprises Section IV of this document. 
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II 


EMPIRICAL MODEL 


A. Introduction 


This portion of the Design Handbook is intended to provide 
a design approach for estimating combustion performance, chamber heat flux 
and stability of GH 2 -GO 2 rocket engines based on empirical data correlations. 
The sources of data for the calculation procedures were from test series 
conducted over the following ranges of operating conditions: 


Injector mixture ratio, 0/F 
Injector orifice area ratio, A^/A 
Chamber pressure, N/cm‘^ (psia) 
Propellant temperature, °K (°R) 


1 to 8 
0.5 to 2 

68 to 340 (100 to 500) 
167 to 333 (300 to 600) 


The injector element types used to obtain these correlations 
and their basic characteristics are identified in Table 1. The accuracy 
of the combustion performance and chamber heat flux correlations are generally 
within the measurement accuracies of the test data used for the correlations 
(See Final Report). The stability characteristics are correlated with injector 
pressure drop for low frequency data, and the high frequency data is tabulated 
with corresponding design and operating conditions. 

The performance correlation of this section of the Design Handbook 
should be used to determine combustion performance predictions only i.e., energy 
release efficiency, ERE. The mixing efficiency, E^, correlations are different 
in some cases if the end product desired is a predicted cold flow E^ and not 
ERE. The cold flow E^ correlations and the reasons for their deviation from 
the corresponding hot fire values are contained in Reference 1. Section III, 
Analytical Model, includes cal cul ational procedures which are applicable to 
obtaining both cold flow and hot fire E^ predictions. 


B. Combustion Performance Prediction 


The calculation procedure is structured such that it analyzes an 
engine point design and predicts a specific combustion performance efficiency. 
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II. B. Combustion Performance Prediction (cont.) 


In general, a desired performance level must be found by iterations of the 
design variables, as shown in the following flow chart. 


Problem Specification 

1 . Injector design and 
operating variables 


Iterate Design 
Variables 


Chamber Length, 



STEP NO. 

1. To begin, a complete specification of the problem must be given. 

The following items must be either known, assumed or calculated: 

(1) Chamber length 

(2) Injection element type 

(3) Mixture ratio 

(4) Injection orifice area ratio 

(5) O 2 and H 2 injection velocity 

(6) Element size 

2. Using Figure 1 and the known chamber length, determine an equivalent 
cold flow length. (For elements not shown on this curve, an engineering judge- 
ment will have to be used to select the most applicable generic curve depicted 
on Figure 1; e.g., for unlike impinging elements such as doublets, pentads, or 
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II. 8. Combustion Performance Prediction (cont.) 

STEP NO. 

quadlets, the triplet curve would be the most applicable). 

3. Compute the mixing efficiency, E^, using the equations and figures 

listed in Table II. E^ for elements that are not represented can 
be approximated by selecting the closest generic element type. For 
swirler elements with less than V^/V^ = 1.0, combustion performance 
can be estimated by using the correlation displayed in Figure 6 for 
the effect of V^/V^ on ERE. 

4- Determine the energy release efficiency (ERE)by using Figure 7 which 

is the correlation between ERE and E^^. For the external impinging 
elements (F-O-F triplet), the mixture ratio dependent curves should 
be used. For all other elements, use the 0/F = 4 curve. 

An example problem, illustrating the use of this technique to calculate 
ERE is enclosed as the last subsection (II. E.) of the Empirical Model. 

C • Heat Transfer Predictions : 

This portion of the design guide provides estimates for the 
combustion chamber heat flux. Use Figures where the estimated to nominal heat 
flux ratio,0 displayed as a function of oxidizer to fuel injection 

momentum ratio. This ratio is the estimated wall heat flux to the heat flux 
based on nominal properties of the combustion products. 

The nominal heat flux, is calculated from the following 

riOril ^ 

equation: 

= 0-026 k/D Re'^ Pr'"^ (T ^ - T ) 

nom '■comb w^ 
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II. B. Combustion Performance Prediction (cont.) 


where: 

D = Local diameter 


Re 

( /^U)g D 

= Reynolds number based on diameter, over- 

all mass flux, and gas properties at 
nominal 0/F and P^ and adiabatic wall 
temperature. 

Pr 

k 

= Prandtl number of combustion gases baseci 

on nominal 0/F and P and adiabatic wal 
temperature. ^ 

(fu)3 = 

Nozzle mass 

flux (one dimensional flow value) 


Thermal conductivity, viscosity, and specific heat of 
combustion gases at nominal 0/F and P^ and adiabatic 
wall temperature. 

^comb 

Theoretical 

combustion temperature 


T = Nominal design gas side wall temperature 

w 

Calculate a nominal heat flux using the above equation or a similar 
technique. Then determine the chamber heat flux ratio, Pgst^^nom 
3.8 cm (1.5 in) from the injector face by referring to Figure 8 for a particular 
element and injection momentum ratio. The estimated value is the product of the 
heat flux ratio and the nominal heat flux: 

,3 = / 0 

*^est I 0 / ^0^ 

\ ’^nom / 

Estimated heat flux values at other chamber stations can be estimated using 
the experimental observation that the heat flux ratio tended to a value of one a^ 
throat plane independent of element type. The 0est^0nom station is 

determined by linearly scaling the head end 0est^0nom as a function of injector 
to throat distance. The following equation can be used to calculate a heat flux 

ratio at any chamber length. 
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D. 


Stability 


The injector designer should be aware that both low and high 
frequency instabilities can occur using GH 2 -GO 2 propellants under the conditions 
specified at the start of this section. Low frequency stability criteria for 
momentum mixing and shear mixing elements are shown in Figure 9A and B. This 
data should be used as an approximate guide to stiffness requirements. The 
high frequency stability data listed in Table III in conjunction with the 
engine system schematic shown in Figure 10, constitutes a consistent set of data 
which can be used by the injector designer m conjunction with existing 
analytical models such as Reference 5. 


E. Sample Problem for Hot Fire Performance (Empirical Model) 

TASK: Predict the following: 

Injector/chamber design required to achieve a ERE = 98% 


GIVEN DESIGN REQUIREMENTS : 


Oxidizer: GG 2 

Fuel : GH 2 


F 


vac 


P 

c 


0/F 


c 


1.000 Ib^ 

300 psia 

4.0 
40:1 


Tp = 540°R 

Injector Allowable 

Oxid: 60 psi (max) 

Fuel: 60 psi (max) 


ASSUMEU PRELIMINARY UESIGN PARAMETERS: 


Element Type: 

*'comb‘ 

F/E: 


Shear Co-Axial 
5 in. 

50 lb 
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II. E. Sample Problem for Hot Fire Performance (Empirical Model) (cont.) 


CALCULATIONS: 


Number of elements N = 


1000 


20 


and 


I 


sp 


vac ode 


453.5 


Ibf-sec 
Ibm ■ 


(Source: JANNAF One Din)ensional 

Kinetic Program, Ref. 2) 


Fixed I Losses = 15 [assumed using JANNAF 

sp lt)m procedure Ref. 3) 

Perfect Injector I_„ = 453.5 - 15 = 438.5 Ibf-sec 

Ibm 


The delivered Specific Impulse ~ ERE x 438.5 = .98 x 438.5 - 430 Ibf-^ 


Solving for weight flows 

Total flow rate (Wj) 


F 

^ spd 


1000 


2.33 


Ibm 

sec 


w? 

Fuel flow rate (W^) = ^T'oTF 


2.33 

5 


.467 


Ibm 

sec 


Oxidizer flow rate (W^) 



2.33 - .467 = 1 .86 


Ibm 

sec 


Total pressure at injection inlet (P^) 


p = P^ + ^P 
0 c 

On oxidizer side use .10 P^ as a AP criteria. 
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II. E. Sample Problem for Hot Fire Performance (Empirical Model) (cont.) 

= .1 P^ = .1 X 300 = 30 psia 

On fuel side use all the available drop i.e., 60 psia. 

AP^ = 60 psia 

One should note that this value of AP . is in the range of expected stability 
problems as shown in Figure 9 of this 'handbook. 

and 



From isentropic relationships the Mach numbers are; 


Mach 

o’ 

O 

X 

II 

.37 

Mach 

"°-fuer 

.515 


and solving for injection velocity: 



The AV ratio is: 

'^f ■ ^0 2220 - 405 

V " 405 

0 

and determining orifice injection area 
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lI.E. Sample Problem for Hot Fire Performance (Empirical Model) (cont.) 


0 *^o ^0 1.86 X 48.3 x 540 _ fv?n 

P V N " 300 X 405 X 20 

C 0 


and diameter (D^) 




From Figure 1 @ of is: 


.160 in. 


and 


•-cold 


.88 in. 


^cold^^o 


.88 _ 

7T60~ " 


5.5 


and from Equation 1 of Table lithe mixing efficiency E^^ is calculated from the 
following expression: 

r 1 

From Figure 2 determine at '^o ~ 

'^0 

K = 14.6 

c 

E = 100 - 26.2 

m 

and from Figure 7 the ERE is 97.6%. If this first try had not agreed with 
the ERE goal, the above procedure would have to be repeated, varying 
preliminary design parameters such as, F/E, or element type until 

agreement had been achieved. 


fff] = 74.4* 
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Equivalent Cold Flow Mixing Length, 

FIGURE 1. RELATIONSHIP BETWEEN CHAMBER LENGTH AND COLD FLOW 
MIXING LENGTH 
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FIGURE 2. COAXIAL AND INCREASED SHEAR COAXIAL ELEMENT CORRELATION 
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FIGURE 5. PKEhlX ELEMEl^T CORRELATION 
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Injector Pressure Drop - Oxidizer Circuit Injector Pressure Drop - Oxidizer Circuit 

(AP/PJ 


SHEAR MIXING ELEMENTS 



Injector Pressure Drop - Fuel Circuit, (AP/P^)p 


PENETRATION MIXING ELEMENTS 



Injector Pressure Drop - Fuel Circuit (AP/P^)^ 


FIGURE 9. LOW FREQUENCY STABILITY CHARACTERISTICS OF GASEOUS INJECTOR 
ELEMENT 






TABLE I GASEOUS INJECTION ELEMENT Cl-iARACTERISTICS 
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GAS/ GAS INJECTOR STABILITY DATA 



Bomb 2 grain RDX 




TABLE III (cont.) 

Triplet Injector ^ 
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ANALYTICAL MODEL 


A. Introduction and Goneral Procedure 

The analytical model described in this section can be used 
by the designer to select the element type and independently vary injector 
and chamber design parameters to determine performance and chamber heat flux 
predictions. The effects of allowable ranges of feed system pressure drops, 
element geometry and size, and hydraulic characteristics also can be evaluated 
by the model. The chamber length parameter can be varied parametrically 
to analyze the effect of length on mixing rate. The chamber contraction ratio 
and design (or experimental test) wall temperature need to be specified only 
if chamber heat flux estimates are required. 


The flow chart shown in Figure 11 illustrates how the general- 
ized analytical model is utilized to predict either cold flow or combustion 
performance and chamber wall heat flux. The upper left block represents 
required model input data which must either be specified by system design require- 
ments or assumed as preliminary design criteria. The Injector, Chamber, and 
Element Selection blocks represent the designer's discretion in selecting the 
design combination which best satisfies all of his requirements. The upper 
central blocks may be bypassed if chamber heat flux estimates or combustion 
performance effects are not required. The bottom group of blocks convert the 
non-dimensional mixing parameters B or B. into oxidizer-rich and fuel-rich 

a U 

stream tube mixture ratios. The two upper right blocks list the tables which 
contain generalized equations used to predict performance and chamber heat 
flux, respectively, for all element types. 


B. Description of Calculational Procedure 

The above outlined procedure can be used as a generalized 
model to calculate local mixing distributions. For the elements specifically men- 
tioned, both oxidizer and fuel calculational procedures are defined step-by-step 
in Tables IV thru IX. For elements not listed, the appropriate generalized 
equations can be developed by following the modeling instructions given in the 
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III.B. Description of Calculational Procedure (cont.) 

final report. Appropriate equations to calculate hydraulic diameters at the 

injector face plane and spray angles are listed on the tables if applicable. 

The tables specify recommended values of interaction exponents and specify 

whether the axi-symmetric correlation or the two dimensional correlation 

B labeled on Figure 12 should be utilized. As a final output from these 

tables, values for mixing correlation parameter, B, are determined. Values for 

B or B and B . or are then used to enter Figure 12 to determine the local n 

ao to af tT , . j. • £ 

ratios,^ j and respectively, which are required in the combustion performance 

and chamber heat flux predictions to define the fuel rich and oxidizer rich 

stream tube mixture ratios. 

If hot-fire performance and/or chamber heat flux values are 
required, then combustion influence parameters (Da) must be included 
in the calculation. For combustion performance problems several iterations may 
be required until the calculated combustion influence parameter converges upon the 
assumed value used in the equations for obtaining the parameter B as defined in 
Tables IV through IX. This process is demonstrated in the sample problem in- 
cluded at the end of this section. 

The combustion influence parameter for ambient temperature 
6O2/6H2 is plotted vs. 0/F in Figure 13, and for 82.6% FLOX/CH^ in Figure 14. 

For other propellant combinations, the following procedure must be followed 
to determine the combustion influence parameter. 

1. One dimensional equilibrium thermochemical data must 
be generated for the particular propellant combination 
and combustion chamber conditions (inlet temperature 
and pressure) over a range of mixture ratios from about 
10”^ to 10^ times stoichiometric proportions. A program 
such as JANNAF ODE is recommended (See Reference 2). 

2. From the data obtained in Step 1, the combustion in- 
fluence parameter can be calculated at each mixture 
ratio from the terms in the right hand side of the 
following equation: 
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III.B. Description of Calculational Procedure (cont.) 

Combustion influence parameter = ^ dA dV dT dMW 

IT = ~T ■ “Mir 

3. The combustion influence parameter is then plotted 

vs mixture ratio as shown for the two propellant com- 
binations used in this section (Figures 13 and 14). 

Using the ^ parameters obtained from Figure 12 the 
performance and chamber heat flux predictions can be determined according to the 
procedures outlined in Ta des X and XI, respectively. The only additional 
data required for the performance prediction procedure outlined in Table X are 
values for theoretical C* and 1^^ vs 0/F for the desired propellant combination. 

It is recommended that this theoretical data be obtained from the JANNAF ODE 
computer program (Reference 2). A sample problem is included in the next section 
to assure the reader that he is following correctly the intended calculation 
procedure. The only additional data required for the chamber wall heat flux 
prediction procedure outlined in Table XI are values for theoretical gas temp- 
erature and gas properties vs 0/F for the desired propellant combination. The 
theoretical gas temperature can be obtained from the Reference 2 JANNAF program and 
the gas properties can be obtained from other sources available to the user such 
as Reference 6 for conventional rocket propellants or Reference 4 for air- 
hydrocarbon reactants. 

After calculating heat fluxes for the oxidizer rich and fuel- 
rich stream tubes as outlined in Table XI and comparing them with the nominal 
uniform mixture ratio value, the designer must visually examine his injector 
pattern layout to identify the local wall mixture ratio. For example, if annular 
fuel co-axial elements are used the wall mixture ratio will be fuel rich. Hence 
the wall heat flux will be more closely approximated by using . If external 

impinging F-O-F triplets are used, the wall 0/F will depend upon element orienta- 
tion next to the wall and whether the fuel underpenetrates (Y < .5 D ) or over- 
penetrates (Y ^.5 D^) the oxidizer jet. In the event of like doublets, pre- 
dicted wall heat flux will depend upon whether a fuel doublet or oxidizer doublet 
is nearest the wall. For a detailed description of how these factors influenced 
chamber heat flux on this program, the reader is referred to the Final Report, 
Reference 1 . 
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c 


SAMPLE PROBLEM FOR COLD FLOW AND HOT-FIRE PERFORMANCE 


TASK: Predict the following performance parameters: 


1. Cold Flow E^ @ L = 1 & 5 in. 

2. Hot Fire E^ and % C* @ L = 5 in. 


GIVEN DESIGN REQUIREMENTS: 

MW 

T 

P 


Oxidizer: 82.6% FLOX 

36.8 

400°R 

0.86 

Ibm/ft 

• 95 X It 
Ibm/f t- 

Fuel: CH^ 

16. 

800°R 

0.19 ^ 

Ibra/ft 

9.0 X 1 
Ibm/ft* 

F = 1,000 lb , P 

vac 1 ^ 

= 100 psla. 

O/F -= 5.75, 

e = 40:1 



ASSUMED PERFORMANCE GOALS: 96^ C*, 927. 

C* (5.75) = 7100 fps; I (e = 40, 5.75) = 420 

s p 


ASSUMED PRELIMINARY DESIGN PARAMETERS: 


Element Type: Shear Co-Ax 



AP 

"d 

N = 20 

T 

W 

Oxid: 

15 psi 

• 

00 

Cn 

L = 5 in. 


Fuel : 

25 psi 

.75 

CR = 3.0 



= 1800°R 


CALCULATIONS : 

Iniecdor Flow Rates 

F 


Wt = 


vac 


« »sp, theo> 


( 1 , 000 ) 

(.92) (420) 


2,59 Ibm/sec 


W 

o 


0/F *, _ ^ — <'2.59') = 2.206 Ibm/sec 

1 + 0/F * t " 1 + 5.75 
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Injector Flow Rates (cont.) 


W 




1 + 0/F 


2.59 

1 -f 5.75 


•384 Ibm/sec 


Chamber Design 




t 


(%C*) (C* ^ ) w 

___ theo t 


P 

c 


g 


•96 (7100) (2.59) 
(100) (32.2) 


5.48 


, 2 
in 


. CR — 5.48 (3.0) = 16.5 in^ 

Ag = . e = 5.48 (40) = 219 in^ 


Injector Desig n 


A P 


2 g 


'I ' ^ 


400 fps 


= 


12 


2 (32.2) 


15 

.19 


1100 fps 


W 


(144 — ) (2.206) 

ft 

20 (.86) (400) (.85) 


N p V 

o o do 


= .0543 in 


144 (.384) 


20 (.19) (1100) (.75) 


= .0176 in 


4 (.0543 in ) 

TT 


= 0.263 in 


Ho 


Let Tube wall thickness (t ) = .010 in. 

w 


t).. - D + 2 (t ) 

I. L O W 


.263 4 2 (.010) = .283 in. 
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“£2 = \ 


2 

°fl ^ — 


= \1 


4 (.0176) 


.283 + 


= ,320 in. 


Hf 


Df 2 - = .320 - .283 = .037 in. 


Re 


P o ^Ho 


d o 


(.86) [.85 (400)] (.263) 
(.95 X 10'^) (12 i^) 


= 6.75 X 


Re 


df 


.19 [.75 (1100)] (.037) 
~ ^ 5T2 in V 

(9.0 X 10 ) ( — ^) 


= 0.54 X 10 


Pd 

o 


Pdf = 


Pd = 
eq 


.86 


2 P o " 2 (32.2) 


400 


n2 


12 in/ft 


= 14#8 psi 


.19 


W 


2 (32 


MW /MWp 

“ MW /MW, 
1 + o t 


= 24,8 psi 


0/F 


Pd 


eq 


(36.8/16) 

14.8 + 5.75 [24.8] 

1 + .40 


= 17.6 psi 


4 (0/F) ^ ^ (5.75) . ^ IQ.O 

MW /MW, 36.8/16. 

O L 

The above parameter is not within 0.5 to 2.0 range. Therefore we must use 
Figure instead of using mole fraction terms in or parameters 

g (lOl = 2.43 

for B B @ (1) 

ao ^ 

U») = .43 

tot i> , B e <i) 

tf t 
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Cold Flow Distribution 



Use the above value of to find from Figure 12. 

= 0.55 j (0/F) = 0/F/ n . = 5. 75/. 55 = 10.4 

J 3 




Use the above value of to find from Figure 12. 

= 0.66; (0/F) = .66 (5.75) = 3.79 = 0/F^ 


(1 - n ^) (1 - . 66 ) 

U = n = ,55 = .415 

(1 - Hj) (1 - ,55) 


- 32 - 



u 


1 + u 


1 + 0/F \ 

1 j = 

1 + 0/F / 


.415 

1 + .415 


1 + 10.4 
1 + 5.75 


= .495 


= I - .495 = .505 


E = 100 1 - X, 


m 



(0/F\ - (0/F) 
1 


(0/F) - (0/F)^ 


1 + (0/F)^ 


_0/F [1 + (0/F)^] 


B 


ao 


B 


tf 

U 




1.0.4 - 5.75‘ 


5.75 

- 3.79 

= 100 [1.0 - . 

495 

1 + 10.4 

-.505 

5.75 

(1 + 3.79^ 

in. (Cold Flow) 






= 37.7 


0/F 

j 

6.04 

= 145 

n. 

= .940 

0/F 

i 

5.40 

= 1.16 

= .560 


= 

.440 

= 100 [^l - .56 

0.29 > 

1 AA 1 _ 

0.35 

Vl _ 

97.2% 

[ 6.75 i 

1 - \ 5. 

75 (6.40:jfJ 


= 76.2 


For Hot Fire Mixing at L = 5 in. 

From above cold flow solution at L = 5-in,, (O/F)^ = 6,04 and (0/F)^ 
5.40. So from Figure 14 


Da = 13.9 at 0/F . ^ 6.04 

o J 

Da^ = 14.0 at 0/F^ - 5.40 

However, the above combustion influence parameters were for ambient (540 R) 
temperature propellants. Thus, they must be corrected for design inlet temp- 
eratures as shown below. 


-: 3 - 



Assume Da (T ) = [ Da (T J] 

o o I- Q ref^ 


Da (400 R) = 
o 


13 9 { 540 °R ] 
V 400OR J 


Da^ (800 "r) 


540°R 

800^R 


Go back and re-calculate hot fire B & B at L = 5 in 

ao tf 




i85 (1 + 18.8) 


.2517, 0.2 

fjTilj 


"j = -313 


0/Fj = 18.4 


Da (540 R) = 8.5 

o 


Da^ (400°R) = 8.5 = 11.5 


Iterate on 0/F and Da 
j o 


Assume Da 


/ 1 + 18.8 y 

V 1 + 14 / 


(5.07) = 6.03 


n. = .415 0/F. = 13.8 

J 2 


Da^ (540^R) = 10.3 


- 13.9 


O.K. (Close enough) 


Similarly for hot fire fuel 


5.0 r„ , 

.Ci'M [O. 43 J 


|T75 (1 + 9.4_)P*^^ [.054j‘^^||^j‘ 


rij, - .32 


0/F. = 1.84 


Da^ (540 R) = 7 g 


D“f woo”') = 8.6 
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6.0 


On next iteration assume Da^ 

1 + 9.4 


'tf 


1 + 6.0 


1.25 


(9.25) = 15.2 




Da. 


= .45 o/F^ = 2.59 


Da^ (540°R) = 9.5 


(800°R) 


_ 9 5 = 6.4 O.K. (Close enough) 

“ 800 


U 


= .415 


( 1 - -45 ) 

(1 - .415) 


= .39 


.39 

Xj = 1.39 


1 + 13.8 
1 + 5.75 


X. 

J 


= .615 


X. 

1 


E = 100 

m 


= .385 

13.8 - 5.75 


( i!#T-r59) 


= 60.7% 


Notice the difference from 97.2% cold flow to 60,7% hot fire with same 
propellants and chamber length. Combustion E^ at 5 in. is less than 76.2% 


cold flow E^ at 1 in. 


%C* = 


( X^ ) l^c* (0/F)j + ( x^ ) ^C* (0/F)j|^ 


C* (0/F) 


.615 [5540] + .385 [6600] 


%C* = 


7100 


5950 


7100 


= 83.8% 


Note that for this initial preliminary injector design that the predicted 
combustion performance is significantly lower than the 967. design goal assumed 
At this point 3 options are available to the designer: 

1. Increase numbers of injection element, N >-20, to decrease 

hydraulic diameters to improve efficiency. 


2 


Increase chamber length to improve mixing 


3. 


Select another element type which has higher mixing charac- 


teristics# 


Option 1 may be limited by already narrow fuel annulus gap, t^ = 1/2 
= #0185 in. Option 2 for present injector requires factor of 2.5 
to 3.0 increase in present chamber length and may be impractical. Option 
3 appears most attractive. 


Go back and re—iterate until a satisfactory design is achieved. 

In any event this one hour hand calculation using the generalized charts 
will have prevented the designer from selecting a dismal performing injector in 
spite of its attractive 97,2% which could be expected from experimental cold 
flow tests. 
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Flow Chart for Analytical Model 


























FIGURE 13. COMBUSTION INFLUENCE PARAMETER VS. MIXTURE RATIO GO,/GH„ 






FIGURE 14. FLOX/CH. THERMOCHEMICAL DATA 



o 
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TABLE IV 


SHEAR COAXIAL ELEMENT EQUATIONS 


OXIDIZER 


FUEL 


Diffusion 
Characteristic : 

Hydraulic 
Diameter : 


Ax i- Syrame t r ic 


Ho 



Two-Dimensional 


Hf 


f2 


D 


fl 


COLD FLOW 


Da = 0.0 

o 


Da^ = 0.0 


F = 1.0 

o 


0.5 


HOT FIRE 


Da^ = f (0/F), 

F^ =0.2 

o 


Da = f (0/F), 

i 1 

= 0.0 


B 

ao 



1.25 725~ F 

K1 

K 


NOTES: 


1. If 

use mole fraction correction factor in Figure 15 to evaluate 


4 (0/F) 

MW /MW^ 
o f 


is not ^0.5 and < 2.0 in 


and 


tf 


equations, 


tf 


.8 4- 




4 (0/F) ^ 

MW /MW., 
o f> 


(1 4 


Da^) 






TABLE V 


SWIRLER COAXIAL ELSMENT EQUATIONS 


Diffusion 
Characteristics : 

Hydraulic 

Diameter: 


OXIDIZER PUEL 

Two-Dimensional Two-Dimensional 



COLD FLOW 


Da 

o 

= 0.0 

Daf 

= 0,0 

F 

0 

= 0,5 


= 0.5 


HOT FIRE 



Da 

o 

o 

u-t 

11 

Da^ 

- f (0/F) 

F 

o 

= 0.0 

■'f 

= 0.0 



use mole fraction correction factor in Figure 15 to evaluate 



2* If design L/D < CL/D)* calculated in above equation, use design L/D instead 

of (L/D)* in B and B - equations# 

to tf 


- 43 - 




TABLE VI 




2 . 


If design L/D < (L/D)* calculated in above equation, use design L/D 
instead of (L/D)* in and equations. 
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TABLE VII 


SHOVrERHEAD ELEMENT EQUATIONS 


Diffusion 
Characteristics s 

Hydraulic 

Diameter; 


OXIDIZER 

Axi-Symmetric 


Circular: D, 


Ho 


Non-Circular: D, 


Ho 


2 (a + b) 


FUEL 


Ax i- Symmetric 


r 

1/2 


4 A 

o 

n 

4Af -j 

TT 

*^Hf 




^ J 

4 ab 


4 cd 




Hf 


2 (c + d) 


COLD FLOW 


Da 


0.0 

0.1 


Da, 


= 0.0 

= 0.1 


Ba 


HOT FIRE 


Da == f (0/^^ 


= 0.0 


‘■■'“ho' 


Da, 


Jj. 

MW 


f 

(0/F)l °- 


= f (0/F) 


= 0.0 


45 


^do 


-1 .625 .25 F 

] M [t-] ■ 


L/D. 


Hf 



- ^ 

r.i3 


Ba^ 

4 (0/F) 



MW /MW^ 


l^df 


0 f 




..625 . 

.25 


bfi 


Pdf 

J 1-10^ J 


Pd 

- eq_ 


NOTE; 


1. If 


r±jo£Fi.i 

I^MWo/MWf J 


is not ^ 0.5 and ^ 2,0 in and Ba^ equations^ 


usa mole fraction correction factor in Figure 15 to evaluate 




TABLE VIII 


PREMIX (TRIPLET & PENTAD) ELEMENT EQUATIONS 


Diffusion 
Characteristics : 

Hydraulic 

Diameter: 


OXIDIZER 


Ax i - Syimne t r ic 



FUEL 


Axi -Symmetric 



Ba 

o 


COLD FLOW 

= 0.0 

Fo =0.1 

HOT FIRE 

= f (0/F)j 

Fo =0.0 


Da^ = 0.0 

F^ =0.1 

Da^ = f (0/F)^ 

F^ = 0.0 
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TABLE VIII (cont.) 

PREMIX (TRIPLET & PENTAD) ELEMENT EQUATIONS 




.65 


P-if + “^eup 


(L/Dh£) 

1 

1 



Baf 

4 (0/F) 

.13 

<=df ' 

.625 

'^%f~ 

.25 

r^'^f ” 

1 

MW /MW, 
o t 



_io' _ 


Pd 

eq_ 


If @ f [ Ba^ ] ^ Hpi . 


NOTES ; 

1. If 


US£ 


^ (Q/^) is not 7 0.5 and ^ 2.0 in Ba 

MW /MW. 
o t 

mole fraction correction factor in Figure 15 


and Ba^ 
to evaluate 



Equations , 
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TABLE IX 


Diffusion 
Characteristics : 

Hydraulic 

Diameter: 


EXTERNAL IMPINGING TRIPLET ELEMENT EQUATIONS 
OXIDIZER 


Two-Dimensional 


Ho 


FUEL 


Two-Dimensional 


1 

o 

C 

1 

1/2 

pAf] 

L - J 

*^Hf 



1/2 


COLD FLOW 


Da = 0.0 

o 

Fo =0.1 


Da^ = 0.0 


= 0.62 




= 0.1 
0.1 


HOT FIRE 


Da = f (0/F) 

° J 

Fo =0.0 


0.62 


1 + Da 


Ip — sin 


Da^ = f (0/F). 



= o.O 


/V \ 

sin 


V°Ho J 

Pd 

o 


> 

1 

tan ”| 

V 

o 

'^f 

, 0/F 

V 

COS 

o 



0.5 


Om 


Om 


where 


HIO 


= 0.0 if 

°Hf ^ ®Ho 



6 - sin 6 



= TT 

“Hf ^ “„o 



B = 7r - 

- 2 Sin-1 



= exp 

Om (6*0 - 2*0 In 

e 

4 (0/F) 

1)' 

MW /MW, 
o f 

!/ 

I 
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iniitt, lA (.cone.; 





If t [ ^ H2 
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TABLE IX (cont.) 






TABLE X 


generalized performance equations 


n. = f 
J 

f»ao 

" "tol 

ni = f 

l«af 

or 

(0/F) 

J 

= 

(o/F)/n. 

(0/F)^ 

= 

nj*(o/F) 

U 

= 

(1-n 



J • - 

(1 -n 

X *” j 


u 

J 

1 + 

u 

X. = 

1 

1 - 

X. 

J 


1 + (0/F) 


1 + (0/F) 


E = 100 

m 


X 

1 

(0/F)j - (0/F)" 

■ ^1 

(0/F) - (0/F)^ 

__ 1 + (0/F)j 

(0/F) [1 + (0/F)^] 

[ • 

"J* \ • 

n± ] 



%C* = 


• C* (0/F)| + \ • C* (0/F)^ 

C* (0/F) 


%ERE = 




I^p (0/F) 
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TABLE XI 


GENERALIZED CHAMBER 

HEAT FLUX EQUATIONS 


Stream Tube 

Nominal 

Oxid. -Rich 

Fuel-Rich 

Mixture Ratio 

(0/F) 

(0/F)j 

(O/F)^ 

Gas Temperature 

T 

gas 

T 

T 

gas^ 

Specific Heat, C 

P 

C 

P 

G . 
PJ 

^Pi 

Viscosity, y 

M 



Film Coefficient 

h 

h 

h 


g 

gj 

gi 

Thermal Conductivity, k 

k 

k. 

J 


Heat Flux 

<P 




Pr 



T 

gas 


f [O/F] 


<!> 
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IV. 


<;ymbqls and abbreviations 


A 

= 

injector element area 

B 

= 

non-dimensional mixing parameter 

a,b,c,d 

= 

rectangular orifice dimension 

C 

= 

specific heat coefficient 

CR 


contraction ratio 

C* 

= 

characteristic velocity 

D 

= 

diameter 

Da 

= 

combustion influence parameter 

d 

= 

differential 

E 

= 

element, efficiency 

ERE 

= 

energy release (combustion) efficiency 

F 


thrust, interaction exponent 

g 

= 

gravitational constant 

H 

= 

hydrogen 

I 

= 

impulse 

JANNAF 

= 

Joint Army NASA Navy Air Force 

K 

= 

empirical factor 

k 

= 

thermal conductivity 

^col d 

= 

equivalent cold flow length 

1 VI 

^comb 

In 

= 

combustion chamber length 

= 

natural logarithm 

N 

= 

number of elements 

0 

= 

oxygen 

ODE 

= 

one-dimensional equilibrium 

Om 

= 

oxidizer misimpingement fraction 

0/F 

= 

oxidizer to fuel mixture ratio 

MW 

= 

molecular weight 

P 

= 

pressure 

Pd 

=: 

dynamic pressure 

Pr 

= 

Prandtl number 

R 

= 

gas constant 

Re 

= 

Reynolds number 

T 

= 

temperature 

t 

= 

wall thickness 
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Symbols and Abbreviations (cont.) 


u 

= 

non-dimensional parameter 

V 

= 

injection velocity 

w 

= 

premix slot width (see Table I) 

w 

= 

flow rate 

WP 

= 

wetted perimeter 

X 

= 

length from throat, mass fraction 

Y 

= 

penetration 

ot 

= 

impingement half angle 

B 

= 

non-dimensional parameter 

Y 

= 

ratio of specific heats 

A 


differential 

e 

= 

nozzle exit area ratio 

n 

= 

mixing ratio 

y 

= 

viscosity 

4> 

= 

heat flux, spray fan propagation angle 

7T 

= 

3.1416 

P 


density 

Subscripts 

a 

= 

axial, axisymmetric chamber, coaxial 

comb 

= 

combustion 

cup 

= 

mixing cup 

d 

= 

discharge, diameter , delivered 

e 

= 

exit 

est 

= 

estimated 

eq 


equivalent (see Table I, Comments column), 
equivalent dynamic pressure 

f 

= 

fuel 

H 

= 

hydraulic 

h 

= 

head end 

i 

= 

refers to fuel -rich streamtube 

J 

= 

injector 

j 

= 

refers to oxidizer-rich streamtube 
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Symbols and Abbreviations (cont.) 


m 

nom 

P 

w 

0 

ode 

s 

sp 

T 

t 

theo 

vac 

X 


mixing 

nominal 

at constant pressure, propellant, premix 
wall 

oxidizer 

one-dimensional equilibrium 

swirler 

specific 

total, triplet 

tangential, two-dimensional, thrust 

theoretical 

vacuum 

length 
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